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Abstract 

The ANTARES detector, completed in 2008, is the largest neutrino telescope in the Northern hemisphere. It is 
located at a depth of 2.5 km in the Mediterranean Sea, 40 km off the Toulon shore. The scientific scope of the 
experiment is very broad, being the search for astrophysical neutrinos the main goal. In this paper we collect the 22 
contributions of the ANTARES collaboration to the 32nd International Cosmic Ray Conference (ICRC 2011). At this 
stage of the experiment the scientific output is very rich and the contributions included in these proceedings cover the 
main physics results (steady point sources, correlations with GRBs, diffuse fluxes, target of opportunity programs, 
dark matter, exotic physics, oscillatinos etc.) and some relevant detector studies (water optical properties, energy 

reconstruction, moon shadow, accoustic detection, etc.) 
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Recent Results from the ANTARES Deep-sea Neutrino Telescope 



Aart Heijboer\ ANTARES collaboration 

^Nikhef, Amsterdam 

aart.heijhoer@nikhef.nl 

Abstract. The ANTARES observatory is currently the 
largest neutrino telescope in the Northern Hemisphere. 
Located at a depth of 2.5 km in the Mediterranean 
Sea, it aims to detect high energy neutrinos that are 
expected from cosmic ray acceleration sites. The status 
of the experiment will be discussed, including a broad 
target-of-opportunity program. The latest results will 
be presented, including searches for a diffuse high- 
energy cosmic neutrino flux, neutrinos from Gamma 
Ray Bursts, and for (galactic) point-like sources. 

I. The ANTARES Neutrino Detector 

Cosmic Rays are thought to originate in Galactic and 
extra- Galactic sources that accelerate protons and other 
nuclei up to high energies. Identification of the responsible 
objects could be achieved by detecting the distinct signa- 
tures of these cosmic accelerators, which are high energy 
neutrinos and gamma rays produced through hadronic 
interactions with ambient gas or photoproduction on in- 
tense photon fields near the source. While gamma rays 
can be produced also by directly accelerated electrons, 
the detection of high-energy neutrinos from these objects 
would provide unambiguous and unique information on 
the sites of the cosmic accelerators and hadronic nature 
of the accelerated particles. 

The ANTARES Collaboration has constructed a neutrino 
telescope 1 1 1 at a depth of about 2475 meters, offshore 
Toulon, France. Neutrinos are detected by photomultiplier 
tubes (PMTs), housed in pressure resistant glass spheres, 
which are regularly arranged on 12 detection lines. Each 
line accommodates up to 25 triplets of PMTs, located 
between 100 and 450 m above the sea bed. The lines are 
connected to the shore via a junction box and a single, 
40 km electro-optical cable, which provides both power 
and an optical data link. On shore, a computer farm runs 
a set of trigger algorithms to identify events containing 
Cherenkov light from high energy muons within the 
data stream, which otherwise consists mostly of signals 
from radioactive decay and bioluminescence. The selected 
events are stored for offline reconstruction. In 2007, the 
first 5 detector lines became operational, followed, in May 
2008, by the completion of the full 12-line detector. 



The reconstruction of muon tracks is based on the arrival 
time of the Cherenkov photons on the PMTs. For high 
energy neutrinos, the angular resolution is determined 
by the timing accuracy, which is limited by the transit 
time spread of the PMTs (1.3ns). Time calibration is 
performed by a number of independent systems, including 
LED and laser beacons Q located throughout the detector. 
The relative inter-line timing has been calibrated using 
the time residuals measured in a large number of down- 
going reconstructed muon events, in addition to the optical 
beacon systems. The positions of the PMTs vary with 
time because of the sea currents. Using an acoustic posi- 
tioning system, combined with information from internal 
compasses and tiltmeters, the positions of the PMTs are 
determined every 2 minutes with an accuracy of ^ 10 cm. 

Most of the analyses described here use a muon track 
reconstruction algorithm (based on |3J) that consists of 
multiple fitting steps. The final step is based on a full 
likelihood description of the arrival times of the detected 
Cherenkov photons, which also accounts for background 
light. The achieved angular resolution is, by necessity, 
determined from simulations. However, several aspects 
of the simulations were confronted with data in order 
to constrain the possible systematic effects in the timing 
resolution that would result in a deteriorated angular 
resolution. The angular resolution (i.e. the median angle 
between the neutrino and the reconstructed muon) was 
found to be 0.4 ± 0.1 (sys) degrees for the detector with 
all 12 lines operational. Studies of the detector and the 
optical water properties |4| are ongoing and may help to 
further improve and constrain the angular resolution in the 
near future. Moreover, a study to observe the shadow of 
the moon using down-going muons might in the future 
provide additional information on the (absolute) pointing 
accuracy ||5J. 

In the following, a selection of results recently obtained 
by the ANTARES experiment will be summarized; many 
of them are discussed in more detail in dedicated contri- 
butions to this conference. 
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Fig. 1: Upper limit on the diffuse neutrino flux of HE 
neutrinos obtained from the 2007-2009 ANTARES data, 
compared to theoretical predictions 6 and to limits set by 
other neutrino telescopes, (see [6} for references). 

II. Searches for high energy cosmic sources 

Searches for cosmic neutrinos and their sources com- 
prise a main goal of the ANTARES experiment. Various 
searches for high energy cosmic neutrinos have been 
performed using the first years of data. 

A. Search for a diffuse neutrino flux 

A search for a diffuse cosmic neutrino flux has been 
conducted using 334 live-days of data collected in 2008 
and 2009 [6]. Such a flux would result in an excess 
of high energy events over the irreducible background 
of atmospheric neutrinos. A measure of the energy is 
provided by an observable R, which measures the number 
of PMTs that detect multiple photons separated in time. 
The distribution of the R variable agrees well with the 
background-only simulations and shows no evidence for a 
contribution from a cosmic diffuse flux, which would 
result in an excess of high-i? events. Consequently, a 90% 
C.L. limit on such a flux is obtained in the energy range 
20 TeV - 2.5 PeV. The limit is shown in figure [T] together 
with previously published limits from other experiments. 

B. Point source search 

Cosmic point-like source of neutrinos have been searched 
for using 813 live-days of data from 2007 up to and 
including 2010 (7). An earlier version of the analysis is 
described in |8 1. Event selection criteria have been applied 
which optimize both the sensitivity and the discovery po- 
tential. Events are required to be reconstructed as upward- 
going and to have a good reconstruction quality, quantified 
by a variable based on the reduced log-likelihood of 
the track fit, and an angular error estimate better than 
1°. The resulting event sample consists of 3058 neutrino 
candidates, of which ^ 84(16)% is expected to be at- 
mospheric neutrinos (muons misreconstructed as upward- 



going). To search for point sources, the analysis uses an 
unbinned maximum likelihood method, which exploits the 
knowledge on the angular resolutiorQ of 0.5° and the rate 
of background events as a function of the declination. 

Two different versions of the search were conducted: 
in the 'full-sky' search, the full visible sky is searched 
for point sources. In the 'candidate search', neutrinos 
are searched for only in the direction of 24 a-priori se- 
lected candidate source-locations, corresponding to known 
gamma ray objects of interest. Neither search yields a 
significant excess of events over the background: the 
post-trial p-values are 2.5% (for a cluster of events at 
a,S = (-46.5°, -65.0°) for the full sky search and 41% 
for the most signal-like source in the candidate source 
list (HESS J1023-575). Limits have been extracted on 
the intensit}]^ of an assumed E^^ neutrino flux from the 
candidate sources. They are shown in Figure [2] The limit 
computation is based a large number of generated pseudo 
experiments in which systematic uncertainties on the 
angular resolution and acceptance are taken into account. 

These limits are more stringent than those from previous 
experiments in the Northern hemisphere (also indicated in 
the figure) and competitive with those set by the IceCube 
observatory [9| for declinations < —30°. The various 
experiments are sensitive in different energy ranges, even 
though they all set limits on E^'^ spectra. For this spec- 
trum, ANTARES detects most events at energies in a 
broad range around 10 TeV, which is a relevant energy 
range for several galactic source candidates. 

An independent point source analysis was performed 
using a different search method based on the 'EM- 
algorithm' |10|. This cross-check yielded similar results 
as the likelihood based analysis described above. 

The sample of neutrino candidates from the previous 
search |8| has been used for additional studies, which are 
also reported on at this conference: 

• An auto-correlation analysis fTT) was performed in 
order to test for unexpected (larger scale) structures 
in the neutrino candidate sample. No such structures 
were found. 

• An analysis has been performed to search for di- 
rectional correlations between the neutrino candidate 
events and 69 published ultra-high energy cosmic 
rays events detector by the Pierre Auger Observa- 
tory lfT2l . The data were found to contain no such 
correlations. 

III. Multi-messenger Astronomy 

Several analyses are performed in ANTARES, which 
focus on coincident measurement of neutrinos with a 

1. Since part of the data in this analysis was taken by a 5-line 
detector, the resolution is sUghtly worse than the 0.4° mentioned 
earlier for the full detector. 

2. The limits are on 0, which is defined by the the fol- 
lowing expression for the neutrino flux: dN/dE = (j> x 
(£/GeV)-2 GeV-icm-^s-i. 
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references). 
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Fig. 3: The upper limits (solid lines) for 37 GRBs obtained 
by the muon track search for the specified neutrino flux 
models (dashed lines) of gamma-ray bursts. The Umits 
were placed using data taken during 2007, when the 
telescope consisted of 5 detector lines. 



variety of external measurements. A selection is described 
below. 

A. Neutrinos from flaring blazars 

In addition to the time-integrated searches described 
above, a time-dependent point source search has been con- 
ducted to look for neutrinos in correlation to the variable 
gamma-ray emission from blazars measured by the LAT 
instrument on-board the Fermi satellite. By restricting the 
search to the 'high state' (typically 1-20 days) of the 
gamma emission, the background is reduced compared to 
the time-integrated point source search. An analysis using 
60 days of live time collected during 2008 is presented in 
I13J ; no significant excess above the expected background 
was observed. 

B. Neutrinos from GRBs 

Various models predict high energy neutrinos to be emit- 
ted by Gamma Ray Burst events. Restricting the neutrino 
search to the duration (i.e. Tgo) of the GRB virtually 
eliminates all background events. Hence, the detection of 
only a few events could already constitute a discovery. 
Two of such searches 1 14 1 have been performed. The first 
one uses the muon-neutrino channel, exploiting the good 
angular resolution of the detector to demand directional 
correlation in addition to the time. This search has so far 
been performed using 37 GRBs and ANTARES data from 
2007 (5 detector lines). No neutrino events were found in 
the a-priori-defined search cone and limits on the neutrino 
flux were obtained; see figure [3] 

The second search is ongoing and searches for 'shower' 
events, which are the result of a localized energy depo- 
sition in the detector. These events can be produced by 



e.g. electron neutrinos which produce an electromagnetic 
shower, or by neutral current interactions of all neutrino 
flavours. A reconstruction algorithm for these events has 
been developed. The sensitivity of this analysis to GRB 
neutrinos of all flavours is presented in |14J. 

C. Optical follow-up of ANTARES events 

To search for transient sources of neutrinos with an optical 
counterpart, a system has been setup to enable fast optical 
observations in the direction of detected neutrino events. 
A reconstruction algorithm that does not require full align- 
ment information is run online and alerts are produced 
for network of small automatic optical telescopes. Such 
alerts are produced for very high energy neutrinos or for 
multiple neutrinos that coincide in time and direction. 
Since February 2009, ANTARES has sent 37 alert triggers 
to the TAROT and ROTSE telescope networks, 27 of 
them have been followed. First results on the analysis of 
the resulting optical images to search for GRB and core- 
collapse SNe will be shown at the conference lfT3l . 

Another combination of measurements consists of corre- 
lating neutrino events with the signals from the gravita- 
tional wave detectors LIGO and VIRGO. A joint analysis 
is being performed that searches for a gravitational wave 
signal in coincidence with a sample of neutrino candidate 
events detected by ANTARES in 2007 |161. 

IV. Searches for Exotic physics 

ANTARES is also searching for signatures of physics 
beyond the Standard Model. An analysis is performed 
that looks for neutrinos produced by Dark Matter particles 
annihilating in the Sun and the Galactic center Iil7l . 
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Fig. 4: Flux upper limit (90 % C.L.) on up-going magnetic 
monopole for relativistic velocities (0.55 < (3 < 0.995) 
is compared to the published upper limits set by other 
experiments are also shown (see |18| for references) as 
well as the theoretical Parker bound. 



Magnetic Monopoles with masses between 10 and 10 
traversing the detector volume would be detected as a 
very bright track. A search for this signature has been 
conducted [ l8 1 and a limit on the flux of monopoles with 
(3 > 0.55 has been obtained; see Figure |4] This limit is 
more stringent than those from previous experiments. 

Another hypothetical form of matter is formed by Nu- 
clearites: particles composed of strange quark matter. 
The signature is a slow-moving (e.g. lO^'^c) bright point 
traversing the detector. The searches are described in llT9l . 

V. Conclusion 

The first deep-sea neutrino telescope, ANTARES, has 
been taking data for four and a half years now. A large 
number of analysis are being performed, looking for 
astrophysical signals of neutrinos, either stand-alone or by 
looking for coincident observations with a variety of other 
experiments. The geographical position, combined with 
the good angular resolution allow ANTARES to explore, 
in particular. Galactic neutrino sources in the relevant 
energy range. In addition, several analyses are aimed at 
detecting signals from non-standard model particles. 

The successful operation of ANTARES, and analysis of 
its data, is an important step towards KM3NET |20| a 
future km3-scale high-energy neutrino observatory and 
marine sciences infrastructure planned for construction in 
the Mediterranean Sea. 
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Abstract. ANTARES is currently the largest neutrino 
detector on the Northern Hemisphere. It consists 
of a tri-dimensional array of 885 photomultipliers 
arranged on 12 vertical lines, placed at a depth of 
2475 meters in the Mediterranean Sea near Toulon, 
France. The telescope, completed in 2008, detects 
the Cherenkov radiation of muons produced by high 
energy neutrinos interacting in or around the detector. 
Muon tracks are then reconstructed using a likelihood- 
based algorithm. One of the main goals of the exper- 
iment is the search for high-energy neutrinos from 
astrophysical point-like sources. Due to its location, 
ANTARES is sensitive to up-going neutrinos from 
many potential galactic sources in the TeV to PeV 
energy regime. New results from an unbinned method 
as well as the sensitivity of the detector are presented. 



I. Introduction 

One of the main questions in astroparticle physics is 
the origin of high energy Cosmic Rays (CRs). In the 
last decade progress has been made related to energy 
spectrum and composition [1|. However, the origin of 
CRs remains unknown. Many acceleration sites have been 
suggested, such as supernova remnants, microquasars and 
active galactic nuclei Q. The final signature of these 
cosmic accelerators are gamma rays and high energy 
neutrinos produced through hadronic interactions. The 
observation of a point-like source of neutrinos would then 
offer a unique occasion to study the mechanism of CRs 
acceleration. 

A. The ANTARES detector 

The ANTARES detector is located at a depth of 2475 
m in the Mediterranean Sea, 42 km from Toulon in the 
south of France (42°48A^, 6°10i;). It consists of a tri- 
dimensional array of 885 optical sensors arranged on 12 
vertical lines. Each line comprises up to 25 detection 
storeys each equipped with 3 downward-looking 10-inch 
photo-multipliers (PMTs), oriented 45° to the line axis. 
The spacing between storeys is 14.5 m while the lines are 



spaced by 60-70 m. A buoy at the top of the line keep 
them to stay vertical. 

The telescope operates by detecting Cherenkov light 
emitted by charged particles that result from neutrino 
interactions in or around the detector. The arrival time 
and amplitude of the Cherenkov light on the PMTs are 
digitized into 'hits' lU and transmitted to shore. 

II. Data Selection 

The data analysed in this work were collected between 
January 31st 2007 and December 30th 2010. The total 
livetime of the analysis is 813 days of which 183 days 
were with 5 lines, while for the remaining 630 days the 
detector consisted of 9, 10 and 12 lines 

The reconstruction of the muon track is achieved using the 
time and position information of the hits. The algorithm is 
based on a maximum likelihood method |4| where a multi- 
stage fitting procedure is applied in order to maximise 
the likelihood of the observed hit times as a function 
of the muon direction and position. The quality of the 
reconstruction is defined by the variable A, which is based 
on the maximisation of the log-likelihood |4 |. Figure [T] 
shows the cumulative distribution of A for upward-going 
events with the simulated contributions of atmospheric 
muons and neutrinos. Atmospheric muons are simulated 
with the MUFAGE package |4|; neutrinos are instead 
generated with the GENNEU 15J package and the Bartol 
model Q. 

Neutrino candidates events are selected requiring an up- 
ward going track, i.e. zenith angle < 90°, and a value 
for the lambda variable A > —5.2. The latter is obtained 
by optimizing the background reduction and the signal 
efficiency, in terms of the discovery potential. Another 
cut is then applied in order to reject mis-reconstructed 
atmospheric downward going muons using the informa- 
tion of the uncertainty on the reconstructed muon track 
direction obtained from the fit. This value is required to 
be < 1°. The final sample consists of 3058 events; from 
the simulations 84% of them are estimated to be neutrinos, 
while the rest are mis-reconstructed atmospheric muons. 
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Fig. 1: Cumulative distribution of the reconstructed quality 
variable A. The dashed line is for simulated atmospheric 
muons while the solid line corresponds to simulated 
atmospheric neutrinos. The bottom plot shows the ratio 
between data and Monte Carlo. 
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Fig. 2: Cumulative distribution of the angle between 
the reconstructed muon direction and the true neutrino 
direction for simulated upward going neutrinos that pass 
the cuts described in Section |ll] assuming a E:j^ neutrino 
spectrum. 



III. Detector Performance 

The angular resolution and the acceptance of the detec- 
tor have been obtained from simulation. The systematic 
uncertainty on the angular resolution has been computed 
following the procedure described in jg) by smearing the 
hit times according to a Gaussian with a width of at in 
order to artificially deteriorate the simulated timing accu- 
racy. The study leads to exclude an additional smearing 
of 3 ns which was found to be incompatible with data at 
the 2(7 level where a is the uncertainty on the flux model 
ISl . The best agreement between data and Monte Carlo is 
obtained for at — 2 ± 0.5 ns. This value is indeed used 
for the simulations presented in this paper. 

A. Angular resolution 

Figure [2] shows the cumulative distribution of the angle 
between the reconstructed muon direction and the gen- 
erated neutrino direction for neutrino events where we 
assumed an energy spectrum proportional to £^^7^ with 
the neutrino energy. The median of this angular error 
is 0.46 degrees. 

B. Acceptance 

The acceptance for signal neutrinos is also estimated using 
simulations. In the search, we deal with fluxes of the form 
of 

where is the flux normalisation. The acceptance. A, is 
defined as the constant of proportionality between (p and 
the number of selected events. Figure [3] shows exactly 
this proportionality: for a source at a declination of -90 
(0) degrees, A = 8.8(4.8) x lO''' GeV cm^ s. Systematic 
uncertainties on the acceptance are constrained by the 
agreement between the simulated atmospheric neutrino 
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Fig. 3: Acceptance, i.e. the constant of proportionality 
between the normalisation factor for an E^"^ flux and the 
selected number of events. 



sample and data. For the computation of the flux limits 
an uncertainty of 15% is assumed. 

IV. Search method 

The algorithm used on the analysis is based on the 
likelihood of the observed events which is defined : 

i 

where the sum is over the events, is the point spread 
function, i.e. the probability density of reconstructing an 
event i at a distance /3j from the true source declination 
and right ascension Ss,cts', B is a parametrization of the 
background rate, obtained from the observed declination 
distribution of the events; fisig is the mean number of de- 
tected events that the source produces and ^tot represents 
the total number of expected events and M{Nl^^^^) is the 
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Fig. 4: Probability for a 3(7 (dashed and solid lines) and 5a 
(dotted and dashed-dotted lines) discovery as a function 
of the mean number of signal events for the case where 
we use the number of hits information in the likelihood 
(dashed and dotted lines) and for the case where we do 
not use it (solid and dashed-dotted lines) for the fuU sky 
search. In this case the signal was added at a declination 
of -70°. 



probability for an event i to be reconstructed with Nhus 
number of hits (this probability was not included in the 
analysis with 2007 and 2008 data 0). 

In order to compute the test statistic the free parameters of 
the likelihood are maximized. We have now to distinguish 
between the two different analysis presented in this paper: 
in the candidate list search only the /.igig are fitted while 
in the full sky search we have in addition the source 
coordinates (Ss,as) to fit. In both cases the results of 
the fit are the maximum likelihood value C"'^^ and the 
estimates of the free parameters. The test statistic is then 
defined as: 

Q = - (2) 

where £b is the likelihood computed for the background 
only case. The higher Q the more the data are compatible 
with signal. 

Just using the number of hits information in the likelihood 
let us to gain a 25% (22%) factor for the 3 (5) a discovery 
probability as shown in Figure [4] for the full sky search. 

V. Results 

As mentioned above two different analyses have been 
done. The first one is a full sky survey with no assump- 
tions about the source position. In the second analysis, we 
made a search for a signal excess in an a priori defined 
spot in the sky corresponding to the position of some 
interesting astrophysical objects. 
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Fig. 5: Galactic skymap showing the 3058 data events. The 
position of the most signal-like cluster is indicated by the 
circle. The stars denote the position of the 51 candidate 
sources. 



A. Full sky search 

In the full sky search, no significant clusters of neutrino 
candidates were found. The most signal-like cluster is 
located at 0^,6^ = (-46.5°, -65.0°). The fit assigns 
5 events above the background. The value of the test 
statistic for this cluster is 13.0 which yields to a p- value 
of 2.6%. Figure |4] shows a sky map of the selected events 
in galactic coordinates with the location of the most signal 
like cluster. 

B. Candidate list search 

The results of the search in the direction of 51 pre- 
defined candidate sources are shown in Table U None 
of the sources have a significant excess of events. The 
most signal-like candidate source is HESS J1023-575 
where the post-trial p-value is of 41%. Figure [6] shows 
the 90% confidence level limits on ip using the Feldman- 
Cousins prescription flOl and assuming an E^"^ neutrino 
spectrum for each of the source candidates as a function 
of the source declination. The sensitivity of this analysis 
is also presented, defined as the median expected limit 
and resulting in a factor 2.7 better than the one obtained 
with data collected during 2007 and 2008 only ^ . Limits 
from other experiments are also shown. However, it should 
be noted that for this spectrum, ANTARES detects most 
events at energies around 10 TeV while the limits in the 
Southern Hemisphere published by the IceCube Collabo- 
ration 111 IJ apply to the PeV region. 

VI. Conclusions 

A search of high energy cosmic neutrinos has been 
performed. Data were taken during the first four years of 
operation when ANTARES consists of 5 line for most of 
the first year considered and 9, 10 and 12 for the rest. A 
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TABLE I: Results of the candidate source search. The 
source coordinates and the p-values (p) are shown as well 
as the limits on the flux intensity (fpo%CL.^ latter has 
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Fig. 6: Limits set on the normalisation (p of an E^^ 
spectrum of high energy neutrinos from selected candi- 
dates (see Table Also shown is the sensitivity, which 
is defined as the median expected limit. In addition to 
the present result, several previously published limits on 
sources in both the Southern and Northern sky are also 
shown. 



measurement of the angular resolution using MonteCarlo 
data yields to 0.46 degrees. Neither the full sky search 
nor the candidate list search show a significant excess of 
events therefore limits have been obtained on the neutrino 
flux. 
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Abstract. The ANTARES detector, currently the largest 
deep-sea neutrino telescope in the Northern Hemi- 
sphere, consists of a three-dimensional array of 885 
optical modules arranged over 12 detection Unes an- 
chored at a depth of 2475 m in the Mediterranean 
Sea, 40 km offshore from Toulon (France). The pho- 
tomultiplier tubes detect the Cherenkov light induced 
by the charged particles produced in the interaction 
of cosmic neutrinos with the matter surrounding 
the detector. The trajectories of the resulting muons 
are reconstructed with high precision, revealing the 
direction of the incoming neutrinos. The main scientific 
goal of ANTARES is the search for high energy 
neutrinos coming from astrophysical sources. This 
contribution describes a point source search using a 
dedicated clustering algorithm, based on the analytical 
maximization of the likeUhood. The results of shuch 
analysis using four years of data will be presented. 

I. Introduction 

The ANTARES neutrino telescope |2| started data taking 
in 2007 and is fully operational since 2008. Located at 40 
km off the coast of Toulon it consists of 12 detection lines 
anchored to the seabed at a depth of 2475 m and sustained 
vertically by means of buoys. Each line has 25 floors (or 
storeys) composed by a triplet of photomultiplier tubes 
(PMTs) housed in pressure resistant glass spheres called 
optical modules (OMs). The OMs are facing downward at 
45° from the vertical for an increased detection- efficiency 
for up-going neutrinos. 

This three-dimensional photo-detector array detects the 
Cherenkov light emitted by the charged leptons originated 
in the interaction of high energy neutrinos with the matter 
surrounding the detector. The tracks of the produced 
muons can be reconstructed using the position and timing 
information of the hits arriving to the PMTs. An accurate 
timing and position calibration [2] of the detector OMs is 
necessary in order to achieve the best attainable angular 
resolution. 



The main goal of the experiment is the detection of high 
energy neutrinos from extraterrestial origin, and one of 
the most promising ways of establishing their existence 
is the search for point sources. Here, we present such 
a search using data collected between years 2007 and 
2010 for a total of 813 days of live time. In section 2 
the track reconstruction method and data selection criteria 
are described. The detector performance is reviewed in 
section 3. The clustering algorithm applied in this analysis 
is explained in section 3. Finally the search results are 
presented in section 4. 

II. Data selection and track reconstruction 

Data runs used in this analysis were recorded in the first 
four years of detector operation. Taking into account the 
time spent on sea operations (like the deployment of new 
lines) and sporadic data taking problems of the detector, 
the total livetime of the analysis is 813 days; about 77% 
of this data were collected using 9, 10 and 12 detection 
lines, while the remaining 183 days correspond to data 
gathered with the initial 5-lines configuration. 

The reconstruction method |3| is based on the maximiza- 
tion of the likelihood function describing the probability 
density function (PDF) for the residuals, defined as the 
difference between the measured hit time and the ex- 
pected arrival time of the hits. The goodness of the track 
reconstruction is described by the A parameter, which 
is basically the log-likelihood of the fitted track. This 
parameter can be used to eliminate badly reconstructed 
tracks by selecting an appropriate cut on the A value. The 
corresponding cumulative A distribution for events recon- 
structed as upgoing is shown on Figure [T] The contribution 
from the different components of the expected background 
is also included. The simulation reproduces well the data. 

For this analysis atmospheric muons were simulated us- 
ing the MUPAGE package (TU], while the atmospheric 
neutrinos were generated with the GENNEU [5| package 
according to the Bartol model |6|. 

A cut on the quality of the reconstruction at A > -5.2 
was found to be the optimal for the search discovery 
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Fig. 1: Cumulative distribution of the quality of the 
reconstruction parameter for data and MC upgoing events. 



potential using muon tracks recostructed as upward going 
(6 < 90°). Additionally, the uncertainty on the muon 
direction estimated from the fit is required to be < 1°. 
The selected sample contains 3058 events, out of which 
it is estimated from MC simulations that about 84% are 
neutrino events and only 16% downgoing atmospheric 
muons mis-reconstructed as upgoing. 

III. Detector performance 

The two main parameters describing the performance 
of a neutrino telescope are the angular resolution and 
the acceptance. Both parameters are estimated from 
simulations. Figure |2] shows the cumulative angular 
resolution for upgoing neutrino events following an E^'^ 
spectrum and complying the selection criteria described 
in Section 2. The plot shows that roughly 80% of the 
signal events are reconstructed with an angular error less 
than 1°, being the median value of the reconstruction 
error equal to 0.46 ± 0.1°. The uncertainty on this value 
has been computed considering all the effects leading 
to a deterioration of the detector timing resolution |7|. 
In addition, the uncertainty on the absolute orientation 
of the detector, which is estimated to be of the order of 
0.1°, is also taken into account in the limits computation 
(see section 4). 

The acceptance allows us to relate the detected event- 
rate with the neutrino flux of the source at the Earth, 
and it is shown on Figure |3] as a function of the si- 
nus of the declination considering a flux normalization 
(f> = 10^^ GeV^^ cm~^ s^^ . Based on the agreement 
between data and simulations a 15% systematic error on 
the detection efficiency has been considered for the limits 
calculation. 

IV. Clustering method 

Point source clustering techniques try to identify and sep- 
arate events coming from real sources from background 
events. The Expectation-Maximization (EM) algorithm |8| 




Fig. 2: Cumulative angular resolution for E ^ upgoing 
neutrinos selected for this analysis. 
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Fig. 3: Detector acceptance as a function of the sinus of 
the declination. 



used in this analysis is a pattern recognition algorithm 
that analytically maximizes the likelihood in finite mixture 
problems. These mixture models are different groups of 
data described by different density components. In the 
case of a search for neutrino point sources [|9J the mixture 
problem can be expressed as: 



p{x) = TTtgPhgiS) + 7r,gP,,(x;/i; E) • (P"'"*7^; 



ihits 



bg ) 
(3) 

where 7r{,g and Tr^g are the so-called mixing proportions, 
X — (a, (5) is the position of the signal event in equatorial 
coordinates, /i = {fia,IJ-s) and S = {(Ta,as) are, respec- 
tivly, the mean and the covariance vector of the Gaussian 
distribution, and p"'"*'' is the probability for an event to 
be reconstructed using nhits number of hits. 

In this analysis the expected density distributions of 
background and signal events are parametrized. The pdf 
describing the background is obtained from the declination 
distribution of data events, while the source signals are 
supposed to follow a two-dimensional Gaussian distribu- 
tion. 

The EM algorithm works in two steps. In the first step 
called "Expectation" the expected value of the complete 
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Fig. 4: Skymap with the equatorial coordinates of the 3058 
events selected. The position of the most significant cluster 
is denoted by the green circle and the coordinates of the 
candidate list sources are indicated with the stars. 
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Fig. 5: Limits on the E^"^ flux for the 51 sources in the 
candidate list search. Upper limits previously reported by 
other neutrino experiments for both Northern and South- 
em sky are shown. The ANTARES sensitivity computed 
in this analysis is also included. 



data log-likelihood is computed for a given set of param- 
eters. In the "Maximization" step a new set of parameters 
that maximizes the Ukelihood is found. In our case, the 
parameters to be maximized are the two components of the 
Gaussian width, the expected number of signal events and, 
in the full sky search (see next section), the coordinates 
of the signal source. 

After likelihood maximization the so called test-statistic, 
defined as the likelihood ratio of the two mixture models, 
is computed. Lower values of this quantity indicate that 
data is more likely to be produced by the background, 
while larger values are more likely to be produced by the 
presence of the searched signal. 

V. Results 

Two different searches for point sources have been consid- 
ered in this analysis. In the first approach a blind survey is 
done looking everywhere in the whole ANTARES visible 
sky. The second search used a catalog of candidate sources 
to look for presence of signal at particular locations 
in the sky. The candidate list of sources includes both 
galactic and extra-galactic sources known to be gamma- 
ray emitters. The detector visibility and PSF was taking 
into account when defining the list. 

No significant excess of events was found neither in the 
full sky search, nor in the candidate list search. The 
most signal-like cluster was found at {5 = —64.87°, 
a = —46.49°) in the all sky survey. For this cluster a 
2.6% excess p-value is found. 

The locations (in equatorial coordinates) of the most 
significant cluster, the 3058 events selected and the 51 
candidate sources are shown on figure|4] Upper limit s H on 
the E^^ neutrino flux spectrum are reported in table [ijand 
in figure |5] as a function of the declination for the sources 



in the candidate list. The ANTARES sensitivity (defined 
as the median value of the expected limit) it is shown, 
as well as limits reported by other neutrino experiments 
included for comparison. 

VI. Conclusions 

This contribution presented the analysis of 813 days of 
livetime using data collected in the first four years of the 
ANTARES neutrino telescope operation. No statistically 
significant excess of events has been found neither in the 
search using a candidate list of interesting sources, nor in 
the full sky search. The most significant cluster, with a 
post-trial probability of 2.6 % was found at coordinates 
S = -64.87°, a = -46.49°. Some of the most stringent 
limits to i?j7^ flux were obtained for sources located in the 
ANTARES field of view. Using a different search method 
the results presented here are consistent with the main 
analysis |10| reporting upper limits for the 51 candidate 
sources using the Feldman-Cousins prescription. 
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TABLE I: Results for the 5 1 sources used in the candidate 
list search. The equatorial coordinates of the sources, Q 
values, number of fitted signal events, pre-trial p-values 
and upper limits on the E~'^ flux (10~^GeV~^cm~^s~^) 
are shown. 
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Abstract. Clustering of neutrino arrival directions 
would provide hints for their astrophysical origin. 
The two-point autocorrelation method is sensitive to 
a large variety of cluster morphologies and, due to its 
independence from Monte Carlo simulations, provides 
complementary information to searches for the astro- 
physical sources of high energy muon neutrinos. We 
present the analysis of the autocorrelation function as 
a function of the angular scale of data collected during 
2007-08 with the ANTARES neutrino telescope. 

I. Introduction 

The key question to resolve the long standing mystery of 
the origin of cosmic rays is to locate the sources and study 
the acceleration mechanisms able to produce fundamental 
particles with energies orders of magnitude above man- 
made accelerators. Over the last years it has become 
more and more obvious that multiple messengers will be 
needed to achieve this task. Fundamental particle physics 
processes like the production and subsequent decay of 
pions in interactions of high energy particles predict that 
the acceleration sites of high energy cosmic rays are also 
sources of high energy gamma rays and neutrinos. The 
detection of astrophysical neutrinos and the identification 
of their sources is one of the main aims of large neutrino 
telescopes operated at the South Pole (IceCube), in Lake 
Baikal and in the Meditemnean Sea (ANTARES). 

A. The ANTARES neutrino telescope 

Whereas physics data taking started already during the 
deployment phase, the ANTARES detector 11] became 
fully equipped and operational in 2008. The detector 
is composed of 12 detection lines placed at a depth of 
2475m off the French coast near Toulon. The detector 
lines are about 450m long and hold a total of 885 optical 
modules (OMs), 17" glass spheres housing each a 10" 
photomultiplier tube. The OMs look downward at 45° in 
order to optimize the detection of upgoing, i.e. neutrino 
induced, tracks. The geometry and size of the detector 
makes it sensitive to neutrinos in the TeV-PeV energy 
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Fig. 1: Schematic view of the ANTARES detector 



range. A schematic layout is shown in Fig. [T| 



B. Neutrino detection 

The neutrino detection relies on the emission of 
Cherenkov light by high energy muons originating from 
charged current neutrino interactions inside or near the 
instrumented volume. All detected signals are transmitted 
via an optical cable to a shore station, where a farm 
of CPUs filters the data for coincident signals or hits 
in several adjacent OMs. The muon direction is then 
determined by maximising a likelihood which compares 
the times of the hits with the expectation from the 
Cherenkov signal of a muon track. 



C. Astrophysical neutrinos 

Two main backgrounds for the search for astrophysical 
neutrinos can be identified: downgoing atmospheric 
muons which have been mis-reconstructed as upgoing 
and atmospheric neutrinos originating in cosmic ray 
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reconstruction quality A cos(zenith angle) 



Fig. 2: Distribution of the quality parameter A (left plot) and the cosine of the reconstructed zenith angle (right plot) 
of the selected data events (markers) compared to Monte Carlo simulations of atmospheric muons (red, dotted line) 
and atmospheric neutrinos (blue, dashed hne). 



induced air showers at the opposite side of the Earth. 
Depending on the requirements of the analysis both 
backgrounds can at least partially be discriminated 
using various parameters like the quality of the event 
reconstruction or an estimator for the deposited energy ||8l . 

In addition, analysing the reconstructed arrival directions 
of the events allows to search for an excess over the 
uniform atmospheric backgrounds. Despite important ef- 
forts, no clear signature for point sources of astrophysical 
neutrinos has been found so far fS], ||4l, jS], ||6|. Both the 
distribution and morphologies of sources potentially emit- 
ting neutrinos in the TeV energy range are yet unknown 
but are possibly very inhomogeneous with most of them 
being located in the Galactic disk and spatially extended 
(e.g. shell type supernova remnants like RXJ1713 |7|). It 
seems therefore interesting to study the intrinsic clustering 
of the arrival directions of neutrinos. Possible analysis 
biases are naturally avoided as no prior information about 
the potential sources is required. Covering a large angular 
range, i.e. neutrino emission regions of very different 
sizes, this study complements the searches for point like 
sources and, if successful, would provide hints for un- 
derlying, yet unresolved, source morphologies and source 
distributions. 

II. Autocorrelation analysis 

The most commonly used method to detect intrinsic 
clusters within a set of N events is the standard 2-point 
autocorrelation distribution. It is defined as the differential 
distribution of the number of observed event pairs iVp in 
the data set as a function their mutual angular distance 
AfJ. To suppress statistical fluctuations that would reduce 
the sensitivity of the method, we analyse here the cumu- 
lative autocorrelation distribution defined as 

N N 

A/'p(Af7) = ^ ^ i7(A17y - (4) 



where H is the Heaviside step function. 

A. Data set 

The analysed data set has been recorded by the ANTARES 
neutrino telescope in 2007 and 2008. During this period 
the detector was in its construction phase and has been 
operated in various setups ranging from 5 active lines 
at the beginning of 2007 to a fully operational detec- 
tor of 12 lines since mid 2008. After applying a run 
selection removing for example periods without precise 
detector alignment information, the dataset corresponds to 
about 300 effective days. Comprising mainly atmospheric 
muons, about 100 million events were reconstructed with 
the standard ANTARES reconstruction algorithm. Basic 
selection criteria include a cut on the reconstructed zenith 
angle 6 to ensure that only upgoing muon tracks are 
selected (cos (6') > 0) and a cut on the angular uncertainty 
/3 given by the covariance matrix of the final likelihood 
fit (/3 < 1°). The final selection criteria is a cut on the fit 
quality parameter A, which is derived from the value of 
the maximal likelihood itself Before unblinding the data, 
this cut has been optimized by means of MC simulations 
to yield the best average upper limit on the neutrino flux 
in the search for point like sources Q, ||6|. 2190 events 
pass the final criterion A > —5.4. The A and zenith angle 
distributions of events passing all quality criteria (except 
the ones shown in the plot) are shown in figure |2] 

Following eq. |4] the cumulative autocorrelation distribu- 
tion of the selected events has been determined. It is 
shown in the left plot of figure [3] 

B. Reference autocorrelation distribution 

To detect structures in the sky distribution of the analysed 
events we need a reference autocorrelation distribution 
to compare with. This reference has been determined by 
scrambling the data themselves, a method which allows 
to avoid systematic uncertainties introduced by the use of 
Monte Carlo simulations. The scrambling method uses the 
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Fig. 3: Left: Cumulative autocorrelation distribution of the selected dataset (open markers) shown in comparison with 
the reference distribution (red solid line). Right: Performance of the algorithm to detect a single or multiple point-like 
sources. 



local coordinates (zenith and azimuth) and the detection 
time Ti of all selected data events. While keeping the 
pairs of zenith/azimuth for all events in order to avoid 
losing information about possible correlations between 
them, the detection time is drawn randomly from another 
event within the same detector configuration in order to 
keep track of the changing asymmetry of the detector. 
Using all selected events, a randomized sky map with the 
same number of events as in the data and naturally the 
same sky coverage is constructed. This randomized sky is 
then analysed in exactly the same way as the data to derive 
the autocorrelation function. The randomization process is 
performed about 10^ times and the derived autocorrelation 
distributions are averaged in order to suppress statistical 
fluctuations. The resulting reference distribution is shown 
as a red dotted line in the left plot of figure [3] 



C. Comparison between data and reference 

Structures in the sky distribution of our data will show 
up as differences between the autocorrelation distributions 
of the data and the reference distribution. The compari- 
son between them is performed by using the formalism 
introduced by Li&Ma f9l. This formalism results in the 
raw significances of the differences as a function of the 
cumulative angular scale which is shown in figure |4] As 
the comparison is performed bin-by-bin and as we scan 
over different angular scales, this result has to be corrected 
for the corresponding trial factor. We apply the method 
proposed by Finley and Westerhoff f^l and perform about 
10^ pseudo experiments in which the autocorrelation 
distributions of randomized sky maps are compared with 
the reference distribution. The probability to obtain the 
same or higher significance as the maximum deviation 
observed in the data is calculated and given as final p- 
value of the analysis. 



D. Performance and sensitivity 

The performance of the algorithm has been determined 
using mock data sets for which we scrambled the selected 
data events as described above. While keeping the total 
number of events in the toy model constant and taking 
into account the angular resolution, we added predefined 
source structures with various sizes and source luminosi- 
ties. These mock data sets where then analysed in exactly 
the same way as described above. As can be seen in the 
right plot of figure [3] the algorithm is sensitive enough 
to obtain a 3cr evidence in the exemplary case of 3 point 
like sources emitting each about 6 events. This source 
luminosity is at the detection threshold of the dedicated 
search for a point like excess in the same dataset |5|, |6|, 
which underlines the complementarity of the two methods. 

An important free parameter of the analysis is the binning 
of the autocorrelation distribution as it will determine the 
sensitivity to certain angular scales and which is connected 
to the angular resolution. For the used quality selection an 
average angular resolution of 0.5° has been determined 
from Monte Carlo simulations. Using toy simulations with 
various source scenarios an optimal binning of 0.1° has 
been determined. 

E. Results and discussion 

The described analysis has been applied to the 2090 
selected data events recorded by the ANTARES neutrino 
telescope between 2007 and 2008. The uncorrected sig- 
nificance as a function of the cumulative angular scale 
is shown in figure [4] A maximum deviation between the 
data and the reference distribution of 1.1 cr is found for an 
angular scale < 7° . Correcting for the scanning trial factor 
this corresponds to a p-value of 55 % and is therefore not 
significant. 

In the search for the sources of high energy cosmic rays, 
the detection of astrophysical sources of neutrinos may 
play a crucial role. Various experiments are currently 
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taking data or are in a preparatory phase to achieve this 
goal and the recorded data is scrutinized in numerous 
ways in order to extract a maximum of information. We 
presented here the first search for intrinsic clustering of 
data recorded with the ANTARES neutrino telescope. 
The data, taken during the deployment phase of the 
detector, do not show evidence for deviations from the 
isotropic arrival direction distribution expected for the 
background of atmospheric neutrinos and contamination 
by mis-reconstructed atmospheric muons. 
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Abstract. We present the search for the diffuse flux 
of astrophysical muon neutrinos using data coUected 
by the ANTARES neutrino telescope. We introduce a 
novel method to estimate the energy of high-energy 
muons traversing the ANTARES detector and discuss 
detaUed comparisons between data and Monte Carlo 
simulations. Using data recorded in 2008 and 2009 
a search for a high-energy excess over the expected 
atmospheric neutrino background is presented and 
stringent limits on the diffuse flux of astrophysical 
muon neutrinos in the energy range 20 TeV - 2.5 PeV 
are derived. 

I. Introduction 

Despite enormous efforts throughout the last century, 
the mystery of the origin of high-energy cosmic rays 
remains unsolved. Over the last years it became more and 
more obvious that multiple messengers will be needed to 
achieve this task. Fortunately fundamental particle physics 
processes like the production and subsequent decay of 
pions in interactions of high-energy particles predict clear 
links between high-energy cosmic rays and high-energy 
neutrinos as well as gamma rays. The detection of astro- 
physical neutrinos and the identification of their sources is 
one of the main aims of large neutrino telescopes operated 
in ice at the South Pole (IceCube) and in water at Lake 
Baikal and in the Meditemnean Sea (ANTARES). 

A. The ANTARES neutrino telescope 

The ANTARES detector HI became fully equipped and 
operational in 2008. The detector is composed of 12 
detection lines placed at a depth of 2475m off the French 
coast near Toulon. The detector lines are arranged on the 
seabed in an octagonal configuration, covering a base 
of 180 X 180 and are about 450m high. They hold 
a total of 885 optical modules (OM), 17" glass spheres 
housing each a 10" photomultiplier tube. The OMs look 
downward at 45° in order to optimize the detection of 
upgoing, i.e. neutrino induced, tracks. The geometry and 
size of the detector make it sensitive to neutrinos in the 




Fig. 1: Schematic view of the ANTARES detector 

TeV-PeV energy range. A schematic layout is shown in 
Fig.[T] 

B. Neutrino detection and dataset 

The neutrino detection relies on the emission of 
Cherenkov light by high-energy muons originating from 
charged current neutrino interactions near and inside the 
instrumented volume. All detected signals are transmitted 
via an optical cable to a shore station, where a farm 
of CPUs filters the data for coincident signals or hits 
in several adjacent OMs. The muon direction is then 
determined by maximising a likelihood which compares 
the times of the hits with the expectation from the 
Cherenkov signal of a muon track. 

We analyse here data taken with the ANTARES detector 
between 12/2007 and 12/2009. Related to the construction 
and maintenance efforts, this period includes data from a 
detector comprised of 9, 10 and 12 active detection lines. 
Data runs were selected according to a set of basic quality 
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Fig. 2: Left plot: Correlation between the energy estimator R and the generated muon energy after all quality cuts. 
Right plot: Distribution of the R energy estimator of data (markers) and simulated atmospheric muons normalized to 
the data (histograms) applying only loose quality cuts. Good agreement of the shape is found for both downgoing 
muons (a) and events mis-reconstructed as upgoing (b). 



criteria, which require for example low environmental 
background noise. The selection corresponds to a total 
live time of 334 days (136 days with 9 lines, 128 days 
with 10 lines and 70 days with 12 lines). 

II. Diffuse astrophysical neutrino flux 

The measured flux of high-energy cosmic rays has been 
used to derive upper bounds for the expected diffuse 
neutrino flux [2J, [3J. For the TeV to PeV energy range 
considered here, this flux is typically assumed to originate 
from particle interactions at or close to the cosmic ray 
acceleration sites. Although only weakly constrained, the 
neutrino energy spectrum is typically modelled by a 
simple E^^ power law. 

A. Background discrimination 

Two main backgrounds for the measurement of the flux of 
these astrophysical neutrinos can be identified: downgoing 
atmospheric muons which have been mis-reconstructed as 
upgoing and atmospheric neutrinos originating in cosmic 
ray induced air showers at the opposite side of the Earth. 
Both backgrounds can at least partially be discriminated 
using various parameters like the quality of the event 
reconstruction or an estimator for the energy of the 
muon. To optimize the selection criteria detailed Monte 
Carlo (MC) simulations have been used. The atmospheric 
muon flux has been simulated with the MURA^GE pack- 
age |4|. Generated atmospheric neutrinos are weighted 
corresponding to the 'Bartol' parametrisation Q. Due to 
the lack of information on the production of charm mesons 
in high-energy hadronic interactions, the presence of an 
additional component at high energies (above ^ 10 TeV) 
is possible. Among the models considered in |6| the 
Recombination Quark Parton Model (RQPM) was used. It 
gives the largest 'prompt' contribution to the atmospheric 
neutrino flux. Both event types, atmospheric muons and 
neutrinos, are processed with the full ANTARES detector 
simulation and reconstruction chain. Special care has been 
taken to reproduce the changing detector configuration 



during the analysed data taking period and the details of 
the data acquisition by including for example afterpulses 
in the PMT simulations. The simulated instrumental and 
environmental background noise has been extracted for 
each of the detector configurations from a representative 
real data-taking run. 

B. Atmospheric muon rejection 

Atmospheric muons are recorded with the ANTARES 
detector at a rate of several Hz and dominate the detector 
trigger rate. To remove a large majority of them from 
the dataset the selected events have to fulfil the following 
basic quahty criteria: 

• detection with at least two detector lines 

• more than 60 hits available for the reconstruction 

• reconstructed zenith angle < 80°, i.e. upgoing 
tracks 

To fully suppress mis -reconstructed atmospheric muons, 
an additional 2-dimensional cut has been derived. It com- 
bines the quality parameter A which is derived from the 
likelihood value of the track fitting algorithm and the 
number of hits used in the fitting procedure A^hit- The 
events have to pass the selection 
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Applied to Monte Carlo simulations, these cuts com- 
pletely remove the 2-10^ reconstructed tracks induced 
by atmospheric muons and reduce the contribution from 
atmospheric neutrinos (Bartol + RQPM) from 7 • lO'' to 
116 events Q. 

C. Energy estimator 

As the flux of astrophysical neutrinos is expected to follow 
a harder spectrum (oc E^^) than that of the atmospheric 
neutrino background it should become visible as an excess 
of high-energy events. This discrimination requires an 
estimation of the neutrino, or as best approximation. 
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Fig. 3: Left plot: Energy distribution for the Bartol+RQPM atmospheric neutrino flux and a E^^ astrophysical signal 
with arbitrary normalisation after all quality cuts. The horizontal arrow denotes the interval in which 90% of the 
signal events is expected. Right plot: Distribution of the energy estimator R for data (markers), the Bartol atmospheric 
neutrino flux (filled histogram) and the 'prompt' contribution (RQPM model, dashed histogram). The signal at the level 
of the derived upper limit is shown as the full, red line together with the cut at i? > 1.31 (vertical, dashed line). 



the muon energy. Various energy estimators are under 
study within the ANTARES collaboration [81. Here we 
exploit the structure of the arrival times of photons created 
along the muon track at the OMs. This time structure 
is sensitive to the energy as higher energy muons have 
a higher probability to create electromagnetic showers 
along the track. The light emitted by these showers is 
responsible for delayed hits in the OMs with respect 
to the detected Cherenkov photons. A robust parameter 
sensitive to delayed photons is the mean number of hit 
repetitions R within an event. It is calculated by averaging 
the number of hits Ri recorded by an OM within a 
500 ns time window over all OMs contributing to the 
reconstructed muon track. The clear correlation with muon 
energy is shown for simulated events in Fig. |2j left plot. 
The estimator has been extensively studied both on MC 
and on (atmospheric muon dominated) data (e.g. Fig. |2j 
right plot) and an average HWHM resolution after the 
discussed quality criteria of log(£'i.cc/-Etruc) = 0.4 could 
be determined. 

D. Optimisation of the event selection 

The final discrimination between the atmospheric neutrino 
background and the astrophysical signal is achieved by 
a cut in the R variable. This cut has been optimized 
before un-blinding the data by minimizing the Model 
Rejection Factor (MRF) |9|. Taking into account all 
possible fluctuations of the number of background events 
riB, the average upper hmit /igo%("^B) of background-only 
pseudo-experiments is derived. The optimal value of the 
cut in R is then determined by minimizing the MRF given 
as /29o%(?ib)/"-s> where ns is the number of expected 
signal events from an E^"^ test spectrum. The optimum 
has been found when selecting events above i?cut = 1.31. 

The determination of the energy estimator cut allows to 
define the energy range to which the analysis will be sen- 
sitive. We define this range as the interval containing 90% 



of the signal events and obtain 20 TeV < E^ <2.h PeV 
(see Fig. [3] left plot). 

E. Low energy region and systematic uncertainties 

Before un-blinding the data in the defined region R > 
1.31 the distribution of the selected data events below 
this cut has been compared to the corresponding Monte 
Carlo simulations. Whereas 125 events are selected from 
data, the Bartol flux parametrisation combined with the 
RQPM model for the prompt contribution predicts 105 
events. The observed difference is well within the sys- 
tematic uncertainty of the atmospheric neutrino flux given 
as 25 — 30% |10| and has been corrected before un- 
blinding the high-energy region by applying a scale factor 
k = 105/125 = 1.19 to the MC predictions. 

Further systematic uncertainties on the number of ex- 
pected events in the high-energy signal region include 

• the contribution of prompt events derived as maximal 
deviation between the models discussed in |6J of 
events. The maximal value (1.7 events) is used in the 
following. 

• uncertainties in the spectral shape of the atmospheric 
neutrino flux which is related to uncertainties in the 
primary cosmic ray energy spectrum. By varying the 
spectral index of the neutrino spectrum by ±0.1 in- 
dependently below and above « 10 TeV we obtained 
a sys. uncertainty of ±1.1 events. 

• uncertainties in the details of the description of the 
detector and environmental parameters like the angu- 
lar acceptance of the OMs, details of PMT afterpulses 
and water properties like absorption and scattering 
lengths which lead to a total uncertainty of 5%. 

III. Results 

Applying the discussed quality and energy selection crite- 
ria to the simulated dataset gives a background estimation 
of riB — 10.7 ± 2 events in the high-energy region. 
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Fig. 4: The upper limit for a E ^ diffuse high-energy i/^ + flux derived from data of the ANTARES neutrino 
telescope compared to previous measurements and phenomenological upper bound predictions. 



In the analysed dataset of the ANTARES detector 9 
events were selected. Including the systematic uncertainty 
on the background expectations following the method 
discussed in ifTTIl we derive the 90% c.l. upper limit on 
the number of signal events as ^go%{nB) — 5.7. The 
corresponding upper limit on the neutrino flux is given by 
^90% = 4' ■ f^90%/^S^ where no, is the number of events 
expected from the flux <j). We obtain 1(121 : 

E'^(t)gQ% = 5.3 X 10"^ GeV^^s^^sr-i (5) 

As can be seen in Fig. |4] the derived limit is competitive 
with previous results and further constrains models of the 
diffuse flux high-energy z^^ + flux. 
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Abstract. Measurement of the atmospheric muon neu- 
trino spectrum is important since it is a test on 
the expected conventional atmospheric neutrino fluxes 
as well as on the prompt contribution from charm 
decay. In addition to that, it forms the irreducible 
background for high energy cosmic neutrinos searches. 
A number of methods have been developed within 
the ANTARES collaboration for the estimation of the 
energy of muon tracks traversing the detector. These 
include a maximum likeUhood method, an approach 
based on neural networks as weU as an estimator based 
on the energy loss of the muon along its track. The 
reconstructed muon energy can be used to calculate 
the energy spectrum of the atmospheric neutrinos by 
using a singular value decomposition approach on 
regularized unfolding. The most recent results on the 
energy reconstruction and spectrum deconvolution will 
be presented. 

I. Introduction 

The ANTARES detector is a high energy neutrino tele- 
scope deployed in the Mediterranean sea, 40 km off the 
coast of Toulon in France. It consists of a three dimen- 
sional array of photomultiplier tubes designed to detect 
light produced by charged particles crossing the instru- 
mented volume. Its main scientific goal is the detection 
of astrophysical high energy neutrinos. 

When a neutrino interacts with a nucleon in the area 
close to or inside the instrumented volume, it produces 
a detectable signature. ANTARES is optimized to detect 
muons produced by the charged current neutrino interac- 
tion with nuclei in the surrounding medium. The detection 
principle is based on the Cerenkov light produced by 
charged particles when they traverse a medium with a 
velocity higher than the speed of light. This light is 
detected by the photomultiplier tubes, housed into glass 
spheres called optical modules (OMs). The OMs record 
information about the charge and time of each hit. The 
position of each OM is provided by an acoustic position- 
ing system. This information allows for both directional 
H], ||2l and energy reconstruction. 



The majority of neutrino candidate events detected with 
ANTARES is expected to be due to atmospheric neutrinos, 
coming from the decay of pions and kaons produced 
by cosmic ray interactions in the atmosphere. Neutri- 
nos produced by this mechanism constitute the so-called 
"conventional" component of the atmospheric neutrino 
spectrum. At higher energies, the decay lengths of these 
mesons are longer than their path lengths in the atmo- 
sphere thus leading to a suppression of produced neutri- 
nos. The "prompt" decay of charmed mesons, produced 
at high energies, provides an additional contribution to 
the atmospheric neutrino flux at the higher energies end 
of the spectrum fi\. The prompt neutrino spectrum is 
less steep than the conventional atmospheric neutrino 
spectrum. Measurement of the atmospheric neutrino spec- 
trum can therefore provide valuable information on the 
prompt component which is expected to dominate the 
atmospheric neutrino flux at energies above ^ lOTeV. 
Downgoing atmospheric muons, produced by cosmic rays 
interactions in the atmosphere constitute an enormous 
background to neutrino candidate events. This background 
is reduced by selecting well-reconstructed upgoing events 
since the muons cannot cross the Earth. The atmospheric 
neutrinos cannot be rejected this way, resulting in an 
irreducible background for high energy neutrino searches 
of galactic or extragalactic origin. Fortunately, a cosmic 
diffuse neutrino flux is expected to fall as £'^^ while the 
atmospheric neutrino spectrum is steeper, proportional to 
i?^'^^. This allows the extraction of the cosmic neutrino 
flux statistically by measuring eventual changes in the 
slope of the measured neutrino energy spectrum. The 
presence of a prompt component can be established in 
the same way. 

In the present paper, three methods developed within the 
ANTARES collaboration aiming at the reconstruction of 
the muon energy are presented in section |ll] At the end of 
this section their performance is presented and discussed. 
Section [III] contains a description of the unfolding pro- 
cedure that will be used to reconstruct the atmospheric 
neutrino energy spectrum. A test of the method on a 
neutrino Monte Carlo sample is performed. 
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II. Energy reconstruction 

Relativistic muons passing through matter lose energy 
by means of various processes. The most common one 
is ionization of water molecules. Above 1 TeV radiative 
processes start dominating, leading to an almost linear 
dependence of the energy loss per unit length on the 
energy of the muon. Pair production, bremsstrahlung 
radiation and photonuclear interactions are the processes 
that are responsible for the increase of energy losses at 
higher energies. The stochastic nature of these phenomena 
requires the use of a mean energy loss per unit length 
approximation described by 

- ^ = a{E) + h{E)E. (6) 
ax 

The first term in this equation is almost constant and 
accounts for the ionization energy losses while the second 
term describes the stochastic energy loss phenomena that 
dominate the energy losses above 1 TcV. 

The three strategies presented in the following are at- 
tempting to reconstruct the muon energy by basically 
inverting numerically equation |6] and using the amount 
of light detected by the telescope as an estimate for the 
muon energy loss. Technically, the first strategy uses a 
maximum likelihood method, the second a neural network, 
whereas the third uses only analytical approximations. 
In addition to these three energy reconstruction strate- 
gies, an energy estimator based on the average number 
of multiple hits recorded on the same OM during the 
event has been used for the diffuse flux search by the 
ANTARES collaboration 0, ||5]. The results presented in 
this section are expectations from a realistic simulation of 
the detector, including detailed OM electronics response, 
photon tracking and optical background generation. 

A. Maximum likelihood method 

The first method presented here is based on a maximum 
likelihood estimation of the energy of the muon. We define 
a function which gives the likelihood that the observed 
time t and amplitude A of each hit on the OMs are the 
result of a given muon track. The likelihood function is 
defined as 

NoM 

^E) ^ Y[ P,{E), (7) 

i 

with the product running over all OMs within 300 m from 
the track and Pi{E) the probability for individual OMs 
to record a hit of a certain amplitude or not. The form of 
Pi{E) is given by 

nmax 

P{A-{n)) ^ £Pp(n;(n)).Pg(A;n), (8) 

n=l 

P(0;(n)) = e-<"> +Pthrcshoid((n)). (9) 

Equation [8] describes the probability of observing a certain 
amplitude A, given that the expected number of photo- 
electrons is (n). Pp{n\ (n)) is the Poissonian probability 



of having n photoelectrons given that the expectation is 
(n), while the probability of n photoelectrons in the pho- 
tocathode producing an amplitude A is given by Pg{A\ n) 
and is assumed to be a Gaussian. Equation |9] represents 
the Poissonian probability of having no photoelectrons, 
when the expectation is (n) photoelectrons, to which 
the probability P^/j that the produced photoelectrons give 
amplitudes below the threshold of the electronics is added. 
The best energy estimate E is found by minimizing the 
negative logarithm of the likelihood function C{E). The 
energy dependence of the likelihood function is included 
in the expected number of photoelectrons (n) [E), which 
is calculated using the probability density functions of 
the photon arrival times on the OMs. In addition to the 
muon energy, these density functions depend on the track 
geometry i.e. distances of each OM to the track and 
the OM orientation compared to the track. A linear fit 
is applied to determine the relation between the energy 
estimate E and the true muon energy E^. 

B. Artificial neural networks 

A second method used to determine the muon energy is 
based on neural networks. An artificial neural network 
(ANN) is a simulated collection of interconnected nodes 
where each node produces a certain response to a set 
of input signals. A number of input parameters Xi are 
given to the input layer of the network. The values of the 
input parameters are transmitted to the nodes of the next 
layer, where a connection weight is assigned for each 
input Xi connected to a node j. The method described 
here implements a feed-forward neural network i.e. the 
nodes are grouped in layers and information flows only 
in one direction, from the input layer to the output node. 
Three hidden layers are being used in the ANN energy 
reconstruction method. The input to the j*'' node of a layer 
is given by 

Zj = ^ WijXi + Wo, (10) 

i 

where wq is the bias of the node. The summation is 
performed on the nodes of the previous layer, and its 
output is given by an activation function gj{zj). The back 
propagation algorithm is used to determine the weights 
Wij and the bias wo that optimize the performance of the 
neural networks method by minimizing an error function. 
The training is performed on Monte Carlo simulation 
samples. The network is a mapping from the space of 
input variables Xi onto a space of output variables y^. 
The only output in the present method is the energy 
of the muon. The space, time and charge distributions 
of the recorded hits, the number of hits per OM at 
different distances from the muon track as well as various 
geometrical parameters such as the track direction and 
the distance from the track to the center of gravity of 
the detector are used as input parameters. The optical 
background rate has also been proven to be a useful input 
parameter. The ANN energy reconstruction method uses 
a total of 56 input parameters. 
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C. dE/dx estimator 

This method |5 1 is based on the construction of an energy 
loss dE/dx estimator p. The estimator is constructed as 



(11) 



where the sum of all hit amplitudes in the event is divided 
by the muon path length in the detector sensitive vol- 
ume and the detector acceptance e. The sensitive volume 
extends 2.5 effective attenuation lengths A away from the 
instrumented volume. The acceptance e is given by 



1 
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NoM 

E 



a{0) 



(12) 



where Nqm is the number of OMs in the detector and 
a{6) is the photomultiplier tube's angular efficiency. The 
term — describes the number of photons reaching 
the OM after traveling distance rj in water The detector 
acceptance e measures the fraction of light which can 
be seen by the OMs. An interpolation between tabulated 
values of the estimator p and the true muon energy i?^ is 
used to determine the relation between them. 

D. Energy reconstruction performance 

In order to study the performance of an energy recon- 
struction algorithm, one has to look at the mean and 
standard deviation of the log'io e'""" distribution. Since 
the performance is dependent on the muon energy these 
distributions are examined as function of the true energy. 
For all three methods presented here, these distributions 
are well described by a Gaussian fit. The mean of the 
Gaussian fit shows how far away is the energy estimate 
from the true energy while the standard deviation repre- 
sents the resolution. A mean of the fit very close to zero 
and a small standard deviation indicate a high performance 
of the method. The bias and the resolution of the three 
methods as a function of the muon energy are shown in 
figures [T] and [2] For this study the true Monte Carlo track 
has been used, in an attempt to decouple the performance 
of the energy reconstruction from the quality of the track 
fit. The overestimation of the true energy by the maximum 
likelihood and dE/dx methods for lower energies, as 
shown in figure [T] is due to the fact that for muon 
energies lower than a few TeV the energy loss is almost 
constant, making it difficult to distinguish between e.g. a 
100 Gev and a 500 GeV muon. The optical background 
from potassium decay in sea water and bioluminescence 
can affect the performance below the critical energy since 
these low energy events do not produce sufficient light 
and the effect of the optical background noise is more 
pronounced. The resolution of the energy reconstruction 
methods presented here is below 0.5 in logigE [GeV] at 
energies above 5-10 TeV. The maximum likelihood and 
artificial neural networks methods have a more stable 
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Fig. 1: Mean of the gaussian fit on log 
as a function of the true muon energy for the different 
methods. A perfectly unbiased energy reconstruction is 
characterized by zero mean. 
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Fig. 2: Standard deviation of the gaussian fit on logig %^/^ 
distribution as a function of the true muon energy for the 
different methods. A small standard deviation corresponds 
to a good resolution. 



behavior with the neural network approach reaching as 
low as 0.25 at higher energies. 

The impact of using the reconstructed track was esti- 
mated for the maximum likelihood method, using tracks 
reconstructed with the standard ANTARES reconstruction 
strategy IJ) and selected to increase the population of well 
reconstructed tracks and minimize the atmospheric muon 
background events. The selection efficiency, estimated 
from simulation, varies between 10% at 1 TeV and 20% 
at 100 TeV with almost 99% purity. The median angular 
resolution is 0.25° above 1 TeV. The performance of the 
energy reconstruction using the reconstructed track is 
shown in figures [T] and |2j indicated as "ML after cuts" 
in the legend. 

III. Atmospheric neutrino spectrum 

DECONVOLUTION 

The atmospheric neutrino spectrum cannot be recon- 
structed by just assigning each event to the corresponding 
estimated energy bin. This is due to the limited resolution 
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of the energy reconstruction. The atmospheric neutrino 
spectrum is steeply falHng so events with overestimated 
energy will introduce a large distortion on the spectrum 
at high energies. This can be overcome by the use of 
unfolding techniques [7|. The unfolding is performed 
following the method presented in 0, using the package 
RooUnfold [7|. 

The unfolding problem can be formulated as follows. The 
true neutrino spectrum x{E) is distorted by a number of 
different factors. The energy that is measured is only part 
of the muon energy, depending on the part of its track that 
is visible by the detector. An additional distortion comes 
from the fact that the muon carries only a fraction of the 
parent neutrino energy. Finally, the limited resolution of 
the energy reconstruction itself contributes further to the 
distortion of the spectrum. Let y{E) be the vector of the 
measured spectrum, then in matrix notation the problem 
is expressed as y{E) — Rx{y). R is called the response 
matrix and describes the effect of the detector and the 
reconstruction on the true spectrum. The response matrix 
is constructed by means of Monte Carlo simulation. A 
simple matrix inversion turns out to be a naive attempt 
to solve the problem since statistical fluctuations in the 
data vector x{E) will give a highly oscillating solution 
y{E) with large errors. The problem can be located and 
addressed by performing a singular value decomposition 
on the matrix R — USV'^. This is equivalent to ex- 
pressing the solution vector as a sum of terms weighted 
by the inverse singular values of the matrix S. The 
small singular values are responsible for enhancing the 
statistically insignificant coefficients in the solution. These 
coefficients can be damped out by imposing an additional 
constraint on the smoothness of the solution y{E). Since 
there is no reason to expect abrupt irregularities in the 
spectrum, the solution is not allowed to exhibit significant 
bin to bin variations. Various methods exist to decide 
the amount of regularization one should impose on the 
solution in order to have the optimal trade off between 
bias towards the expected solution and the size of the 
covariance matrix. 

The unfolding method was tested using the results from 
the maximum likelihood energy reconstruction on a realis- 
tic 12 line configuration Monte Carlo dataset, equivalent 
to 3 years of data taking livetime. The response matrix 
was constructed assuming the atmospheric neutrino flux 
parametrizations from 0, ifTOl while the fluxes used to 
generate the test data are taken from [TlJ, [12J. After 
unfolding to the neutrino spectrum at the detector level, 
the atmospheric neutrino flux (figure [3]l can be determined 
by including the effects of neutrino propagation through 
the Earth and detector efficiency. The method succeeds in 
reconstructing the "true" Monte Carlo neutrino spectrum. 
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Fig. 3: Application of the unfolding method on a realistic 
Monte Carlo dataset for the determination of the atmo- 
spheric neutrino flux weighted with E^. Systematic effects 
are not included. See text for explanation of the fluxes. 



described. Their precision varies from 0.25 to 0.5 in 
logj^o-^ [Gs^] above 5-lOTeV when the true track geom- 
etry is used as an input. This resolution may be improved 
for particular classes of events, especially in the case of 
the maximum likelihood and dE/dx methods. The results 
of the maximum likelihood reconstruction on a simulated 
dataset have been used to test the unfolding method. 
Systematic uncertainties are under study and have not been 
included in these proceedings. 
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IV. Conclusions 



Three different methods for muon energy reconstruc- 
tion currently under evaluation in ANTARES have been 
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Abstract. Multiple searches for neutrinos in correlation 
with gamma-ray bursts (GRBs) using data from the 
ANTARES neutrino telescope have been performed. 
One search uses data taken during 2007, at which time 
the telescope consisted of 5 detector lines, to measure 
the number of neutrino-induced muons in correlation 
with 37 GRBs that occurred during the data-taking 
period. While no correlations are observed, upper 
limits have been placed on the flux of neutrinos for 
different GRB models. A second search uses an alter- 
native method to identify the shower at the neutrino- 
interaction vertex. This search is particularly sensitive 
to electron-neutrinos, but is able to observe neutrinos 
of any flavor. The sensitivity of this method to measure 
correlations between neutrinos and prompt photons 
emitted by GRBs is presented for a typical neutrino 
flux model. 

I. Introduction 

First discovered in the 1960's, a gammy-ray burst (GRB) 
is observed as an extremely luminous flash of photons 
from a seemingly random point in the sky. The en- 
ergy spectrum of these photons can be parametrized by 
a broken power law of the form dN/dEocE"", where 
a w — 1 below the break energy (typically around an 
MeV) and a ~ — 2 above it [IJ. The large energy output 
of a typical GRB, around lO'^^ erg, combined with the 
GRB occurrence rate, about one per galaxy per million 
years, makes gamma-ray bursts a candidate source of high 
energy cosmic rays |2|. If GRBs are indeed accelerating 
hadrons, then these hadrons can be expected to suffer 
inelastic collisions with particles found in the acceleration 
regions, typically shock waves in relativistic jets |3|, |4|, 
0. Such hadronic interactions would ultimately yield 
mesons that decay to produce high-energy neutrinos. 

The ANTARES telescope seeks to observe such cosmic 
neutrinos. For this purpose, the sea water of the Mediter- 
ranean Sea is used both as a target and Cerenkov medium. 
Neutrinos interacting with a nucleus in or around the 
instrumented volume will produce a hadronic shower at 
the collision vertex and will also produce an energetic 
lepton in the case of a charged-current interaction. These 



relativistic charged particles will produce Cerenkov ra- 
diation which can then be measured by photo-sensitive 
detection units in the ANTARES telescope. Neutrinos 
with an energy of around 10 TeV or higher will yield 
leptons traveling along the same trajectory as the neutrino. 
The source of such a neutrino can then be identified 
through its correlation with the neutrino's arrival time 
and/or direction. 

To measure the light produced by charged particles gen- 
erated in neutrino interactions, the ANTARES telescope 
consists of 12 detector lines that hold photo-sensitive 
units. Each line is anchored to the sea floor at a depth 
of 2475 m and supports 25 triplets of photo-multiplier 
tubes (PMTs) spaced evenly along the line in 14.5 m 
intervals. Further details of the detector configuration can 
be found in Ref. IQ. The detector lines are not rigid 
and are therefore affected by sea currents. The position 
and orientation of each PMT is determined by multiple 
positioning systems, described in reference |7|. 

While the sea water around the telescope is extremely 
dark, ANTARES can and does detect photons that are not 
due to neutrino interactions. The main sources of such 
background signals are (a) down-going muons produced 
by the interactions of cosmic rays with the Earth's at- 
mosphere, (b) the decay of the radioactive isotope ''''K 
present in the sea water and (c) light emitted by Uving 
organisms in the water The latter two sources produce 
random signals in each photo-detection unit at a rate 
which varies with the sea current and is around 60- 
150 kHz per PMT. This background is reduced through the 
application of triggering algorithms that select potential 
neutrino events The so-called atmospheric muons 
form a background which is reduced after the trigger 
as part of a neutrino analysis procedure, as described 
in the following sections. With the search for transient 
sources such as gamma-ray bursts, the vast majority of 
the atmospheric-muon background can be reduced simply 
by demanding that a potential neutrino event correlates in 
time to the a priori known time of the source flare. 
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Fig. 1: The angular resolution for muon tracks passing 
the selection described in the text, as a function of the 
neutrino energy. 



II. Muon Track Search 

Data taken by the ANTARES telescope during 2007 
has been analyzed to search for neutrino-induced muons 
arriving in correlation with prompt photons from gamma- 
ray bursts. The prompt photon signals were detected 
by external satellite experiments, and were mainly long- 
duration GRBs detected by Swift S, flOl. During this 
data-taking period, the ANTARES telescope was still 
under construction and consisted of only 5 detector lines. 
Muon candidates are extracted from raw data using a 
trigger algorithm that searches for a collection of causally- 
connected hits consistent with photons emitted by a muon 
traversing the detector 1 8 1 . The average trigger rate during 
the data-taking period was 1 Hz. 

The muon trajectory is determined by a reconstruction 
package that uses probability density functions (PDFs) 
describing the expected photon arrival times, of both 
signal and background photons, to find the local muon po- 
sition and direction that maximizes the likelihood that the 
observed hit times correspond to the expected times ifTTl . 
lfT2l . Potential neutrino- induced muons are required to 
have an up-going trajectory (to reject muons from atmo- 
spheric cosmic-ray showers) and to have a good fit quality. 
The fit quality is estimated using two parameters, one 
based on the (log of the) value of the likelihood function, 
A, and one based on the angular precision calculated by 
the fit procedure, ct/. 

The distribution of these parameters observed in the data is 
well reproduced by simulations of the background due to 
atmospheric muons and neutrinos. While signal neutrinos 
have been simulated for several different spectra, a single 
set of cuts on the quality parameters is found to provide 
a good separation of signal from background events. The 
angular resolution for simulated muons passing these cuts 
is shown in Fig. [T] 

The good angular resolution for muon tracks allows the 
direction, in addition to the time, of neutrino candidate 



events to be used in the search for correlations with 
gamma-ray bursts. The candidate tracks are required to 
point back to the GRB position to within 2°, and to occur 
during the arrival of prompt photons, i.e. during the T90 
observed by the satellite. With these criteria, the expected 
background rate is found to be about 1.0 x 10^^ Hz. Of 
the remaining background events, 90% are due to (misre- 
constructed) atmospheric muons from showers above the 
detector and 10% to atmospheric neutrinos producing up- 
going muons. 

This search for muons produced by neutrinos in corre- 
lation with gamma-ray bursts has been applied to the 
data taken during the flares of 37 GRBs that occurred 
during 2007. No neutrino candidate events are observed 
in correlation with the GRBs. The limits placed on the 
total Vf^i + flux of these bursts at the 90% confidence 
level, using the Feldman-Cousins recipe [13], for three 
different GRB models are presented in Fig. |2] In addition 
to a general E^^ spectrum, two GRB neutrino spectra 
have been explored, both based on the fireball model. The 
Waxman-Bahcall model is that described in Ref. |[T4l . and 
all GRBs are assumed to follow the same spectrum. The 
model of Guetta et al. is described in Ref. |15|. For 
this model, the neutrino spectrum of each GRB has been 
calculated separately using the (satellite) measured GRB 
parameters. 

The reliability of the simulations used in this analysis has 
been studied in depth. The uncertainty with the largest 
impact arises from the efficiency and angular acceptance 
used for PMTs in the simulations. Reducing the efficiency 
of the tubes by 15% results in a 12% reduction of 
the acceptance of the telescope to E^^ muon-neutrinos, 
without affecting its angular resolution. A degradation 
of the PMT timing of only a few nanoseconds in the 
simulations is found to be incompatible with muon data. 
Such studies allow a conservative estimate of 15% to be 
placed on the uncertainty of the angular resolution of the 
telescope. Uncertainties on PMT timing can also reduce 
the acceptance of the telescope. The aggregate effect of 
the PMT timing and efficiency results in an uncertainty of 
15% on the relative acceptance of ANTARES to cosmic 
muon-neutrinos. 

HI. Collision Vertex Shower Search 

An additional search method has been developed to search 
for neutrinos of any flavor in correlation with gamma-ray 
bursts. This is done by locating the shower(s) produced 
at the vertex of the collision between the neutrino and a 
nucleus in or around the detector. Such a deep inelastic 
scattering will produce a hadronic shower irrespective of 
the neutrino flavor. In the case of electron neutrinos un- 
dergoing a charged current interaction, an electromagnetic 
shower is also produced at the vertex by the resulting 
electron or positron. Such showers occurring during the 
prompt photon emission of a GRB are then sought. This 
analysis is the first to seek to measure neutrinos of 
multiple flavors in the ANTARES data. 
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Fig. 2: The upper limits on the total i^^ 
37 GRBs, obtained by the muon track search, are shown 
by the solid lines for different flux models. Two of the 
flux models are shown by the dashed lines. The expected 
Waxman-Bahcall flux llT4l is shown in black and the 
expected Guetta et al. flux fTS'l is shown in gray. The 
limits have been placed using data taken during 2007, 
when the telescope consisted of 5 detector lines. 
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Fig. 3: The spacial resolution of the shower reconstruction 
algorithm, determined using special calibration data in 
which LED beacons were flashed. 



As the extent of the vertex shower is small compared 
to the spacing between detection units in the telescope, 
photons radiated by charged particles in the shower appear 
to originate from the same point in space and time. A 
simple procedure to reconstruct the position and time of 
the vertex shower has been developed. This reconstruction 
is applied to candidate neutrino events that have passed 
at least one trigger criteria. In addition to the trigger 
algorithm described in Sect. |ll] a second trigger condition 
in which clusters of time-correlated hits in neighboring 
PMT triplets are sought has been applied. This second 
trigger algorithm was implemented during the 2008 data 
taking period. 



The shower vertex position and time is then determined 
by minimizing the difference between the expected photon 
arrival times and the measured hit times. The arrival time 
of photons is assumed to depend only on the distance 
traveled from the shower vertex. The resolution of this 
reconstruction procedure has been studied using special 
calibration data taken by ANTARES in which LED bea- 
cons were used to obtain a burst of photons originating 
from a single point. The spacial resolution of the shower 
reconstruction is found to be around 1 m, as shown in 
Fig. |3] and the timing resolution around 2.4 ns. 

A sample of reconstructed showers from triggered events 
will contain many electromagnetic showers produced by 
the stochastic energy loss of atmospheric muons traversing 
the telescope. Neutrino-induced showers are extracted 
from this background by placing cuts on quality pa- 
rameters. The most discriminating variable is found to 
be the number of direct hits, Noir Hits' the number of 
selected hits having a time residual better than 15 ns. In 
addition, requiring the number of detector lines used in 
the fit to be greater than two is found to further reduce 
the atmospheric muon background. The cut on NjjirHits 
has been determined by finding the value that minimizes 
the flux of a GRB required to discover a source with 
Tgo =100 s at the Sa level in 50% of test cases. The GRB 
flux used during cut optimization is assumed to follow the 
Waxman-Bahcall flux lfT4l . 

While a rare occurrence, a sparking PMT could produce a 
signal with a similar profile to that of a neutrino-induced 
shower. To eliminate such events, showers are rejected 
if they are reconstructed too close to any PMT. This cut 
effectively reduces the fiducial volume of the telescope by 
about 2%, but since the detector is capable of measuring 
high-energy neutrino-induced showers from far outside the 
instrumented volume, the ultimate effect on the sensitivity 
is much less. 

The expected rate of background events passing the full 
set of quality cuts has been estimated from data taken 
during 2007-2008 in which no GRB was observed by 
a satellite experiment. During this time, the ANTARES 
telescope was under construction and data was taken while 
the telescope consisted of 5, 9, 10 and 12 detector lines. 
The average rate of surviving background events over the 

2007 (5 line) and 2008 (9-12 lines) periods are found to 
be 4.4 X 10^^ Hz and 6.8 x 10~^ Hz, respectively. Note 
that the background is much larger than that obtained by 
the muon track search method described in Sect. |ll] due 
to the lack of directional information. With the simple 
shower reconstruction, the field of view of the telescope 
cannot be restricted to the region of sky around the GRB. 

The (Neyman) sensitivity of the vertex shower search 
method to observe the total flux of neutrinos and anti- 
neutrinos of all three flavors from GRBs |14| has been 
calculated. The sensitivity obtained for the 2007 and 

2008 data taking periods, averaged over the viewable sky, 
is shown in Fig. [4] For the 2007 period, the detector 
consisted of 5 lines and the sensitivity is shown for the 
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Fig. 4: The sensitivity of the shower vertex search method 
appUed to the average ANTARES detector as it was during 
the 2007 and 2008 data taking periods to the total neutrino 
and anti-neutrino flux of GRBs. The neutrino spectra are 
assumed to foflow the Waxman-Bahcall spectrum lfT4l . 



same 37 GRBs studied in Sect. For the 2008 period, 
the sensitivity is shown for 65 GRBs and averaged over 
the 9, 10 and 12 line detector configurations. 

The impact of systematic uncertainties in the detector 
simulations are under study. A significant effect comes 
from the possibility that the efficiency of each PMT 
is not accurately reproduced in the simulations. This is 
being studied by reducing the efficiency with which signal 
photons are converted to hits. Other simulation properties 
have also been studied by varying the Cerenkov light 
propagation, altering the high energy shower generation 
and using different models of the neutrino-nucleus cross 
sections. These effects lead to a reduction of the sensitivity 
of the shower analysis that is currently estimated to be 10- 
20%. 



37 selected GRBs of 2007. Limits on the maximum flux 
that would have only a 10% chance to produce no events 
in the telescope have been calculated. 

The second method seeks to measure the shower generated 
at the neutrino collision vertex and is the first analysis 
of ANTARES to be sensitive to neutrinos of any flavor. 
Showers produced by the stochastic energy loss of atmo- 
spheric muons are rejected by studying the topology of the 
shower event. Due to the lack of directional information 
available from the shower reconstruction, the background 
rate of this method is higher than that of the muon 
analysis. The sensitivity of the ANTARES telescope as it 
was during 2007-2008 to showers from a flux of neutrinos 
from GRBs has been determined. 
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IV. Summary 

Two methods have been developed that can be used 
to search for a detectable flux of neutrinos emitted by 
gamma-ray bursts. Both methods attempt to measure 
neutrinos aiTiving at the Earth in conjunction with the 
prompt photons emitted by GRBs (as observed by exter- 
nal satellite experiments). For both methods, muons and 
neutrinos generated by the interaction of cosmic -rays with 
the Earth's atmosphere form a substantial background that 
must be reduced. 

The first method reduces this background by measuring 
the trajectory of the muon produced by muon-neutrino 
interactions. The muon is demanded to be well recon- 
structed and to come from within 2° of the GRB. It has 
been applied to data taken by the ANTARES telescope 
during 2007, when the experiment consisted of 5 detector 
lines. No neutrino events are observed in correlation with 
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Abstract. The ANTARES telescope is well suited to 
detect neutrinos produced in astrophysical transient 
sources as it can observe a full hemisphere of the sky at 
all the times with a duty cycle close to unity. The back- 
ground and point source sensitivity can be drastically 
reduced by selecting a narrow time window around 
the assumed neutrino production period. Radio-loud 
active galactic nuclei with their jets pointing almost 
directly towards the observer, the so-called blazars, 
are particularly attractive potential neutrino point 
sources, since they are among the most likely sources 
of the observed ultra high energy cosmic rays and 
therefore, neutrinos and gamma-rays may be produced 
in hadronic interactions with the surrounding medium. 
The gamma-ray light curves of blazars measured 
by the LAT instrument on-board the Fermi satellite 
reveal important time variability information. A strong 
correlation between the gamma-ray and the neutrino 
fluxes is expected in this scenario. 

An unbinned method based on the minimization of 
a likelihood ratio was appUed to a subsample data 
coUected in 2008 (61 days live time). By looking for 
neutrinos detected in the high state periods of the 
AGN light curve, the sensitivity to these sources has 
been improved by about a factor 2 with respect to 
a standard time-integrated point source search. First 
results on the search for ten bright and variable Fermi 
sources are presented. 

I. Introduction 

The production of high-energy neutrinos has been pro- 
posed for several kinds of astrophysical sources, such 
as active galactic nuclei, gamma-ray bursters, supernova 
remnants and microquasars, in which the acceleration of 
hadrons may occur. Neutrinos are unique messengers to 
study the high-energy universe as there are neutral and 
stable, interact weakly and travel directly from their point 
of creation in the source without absorption. Neutrinos 
could play an important role in understanding the mecha- 
nisms of cosmic ray acceleration and their detection from 



a source would be a direct evidence of the presence of 
hadronic acceleration in that source. 

Radio-loud active galactic nuclei with their jets pointing 
almost directly towards the observer, the so-called blazars, 
are particularly attractive potential neutrino point sources, 
since they are among the most likely sources of the 
observed ultra high energy cosmic rays and therefore, neu- 
trinos and gamma-rays may be produced in hadronic in- 
teractions with the surrounding medium [l]. The gamma- 
ray light curves of blazars measured by the LAT instru- 
ment on-board the Fermi satellite reveal important time 
variability information on timescale of hours to several 
weeks, with intensities always several times larger than the 
typical flux of the source in its quiescent state |2 1. A strong 
correlation between the gamma-ray and the neutrino fluxes 
is expected in this scenario. 

In this paper, the results of the first time-dependent search 
for cosmic neutrino sources in the sky visible to the 
ANTARES telescope are presented. The data sample used 
in this analysis is described in Section 2, together with 
a discussion on the systematic uncertainties. The point 
source search algorithm used in this time-dependent anal- 
ysis is explained in Section 3. The results are presented in 
Section 4 for a search on a list of ten selected candidate 
sources. 

II. ANTARES 

The ANTARES collaboration has completed the construc- 
tion of a neutrino telescope in the Mediterranean Sea with 
the connection of its twelfth detector line in 2008 IS). 
The telescope is located 40 km on the southern coast 
of France (42°48'N, 6°10'E) at a depth of 2475 m. It 
comprises a three-dimensional array of photomultipliers 
housed in glass spheres (optical modules), distributed 
along twelve slender lines anchored at the sea bottom 
and kept taut by a buoy at the top. Each line comprises 
up to 25 storeys of triplets of optical modules (OMs), 
each housing a single 10" PMT Since lines are subject to 
the sea current and can change shape and orientation, a 
positioning system comprising hydrophones and compass- 
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tiltmeters is used to monitor the detector geometry. The 
main goal of the experiment is to search for neutrinos 
of astrophysical origin by detecting high energy muons 
(>100 GeV) induced by their neutrino charged current 
interaction in the vicinity of the detector. 

The arrival time and intensity of the Cherenkov light 
on the OMs are digitized into hits and transmitted to 
shore, where events containing muons are separated from 
the optical backgrounds due to natural radioactive decays 
and bioluminescence, and stored on disk. A detailed 
description of the detector and the data acquisition is given 
in lUl |4|. The arrival times of the hits are calibrated 
as described in reference |5|. The online event selection 
identifies triplets of OMs that detect multiple photons. At 
least 5 of these are required throughout the detector, with 
the relative photon arrival times being compatible with the 
light coming from a relativistic particle. Independently, 
events were also selected which exhibit multiple photons 
on two sets of adjacent, or next to adjacent floors. 

The data used in this analysis corresponds to the period 
from September 6th to December 31st, 2008 (54720- 
54831 modified Julian day), taken with the full detector. 
Some filtering has been applied in order to exclude periods 
in which the bioluminescence-induced optical background 
was high. The resulting effective life time is 60.8 days. At- 
mospheric neutrinos are the main source of background in 
the search for astrophysical neutrinos. These neutrinos are 
produced from the interaction of cosmic rays in the Earth's 
atmosphere. Only charged current interactions of neutrinos 
and antineutrinos were considered. An additional source 
of background is due to the mis-reconstructed atmospheric 
muons. The track reconstruction algorithm derives the 
muon track parameters that maximize a likelihood func- 
tion built from the difference between the expected and 
the measured arrival time of the hits from the Cherenkov 
photons emitted along the muon track. This maximization 
takes into account the Cherenkov photons that scatter in 
the water and the additional photons that are generated by 
secondary particles (e.g. electromagnetic showers created 
along the muon trajectory). The algorithm used is outlined 
in |6|. The value of the log-likelihood per degree of 
freedom (A) from the track reconstruction fit is a mea- 
sure of the track fit quality and is used to reject badly 
reconstructed events, such as atmospheric muons that 
are mis-reconstructed as upgoing tracks. Neutrino events 
are selected by requiring that tracks are reconstructed 
as upgoing and have a good reconstruction quality. In 
addition, the error estimate on the reconstructed muon 
track direction obtained from the fit is required to be less 
than 1°. 

The angular resolution can not be determined directly 
in the data and has to be estimated from simulation. 
However, the comparison of the data and MonteCarlo 
simulations from which the time accuracy of the hits has 
been degraded has yielded a constrain on the uncertainty 
of the angular resolution of the order of 0.1° |8|. Figure [T] 
shows the cumulative distribution of the angular difference 



between the reconstructed muon direction and the neutrino 
direction with an assumed spectrum proportional to E:j'^, 
where is the neutrino energy. For this period, the 
median resolution is estimated to be 0.4 +1- 0. 1 degree. 
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Fig. 1: Cumulative plot of the distribution of the angle 
between the true Monte Carlo neutrino direction and the 
reconstructed muon direction for 
events selected for this analysis. 
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III. Time-dependent search algorithm 

This time dependant point source analysis is done using 
an unbinned method based on a likelihood ratio max- 
imization. The data is parameterized as a two compo- 
nents mixture of signal and background. The goal is to 
determine, at a given point in the sky and at a given 
time, the relative contribution of each component and to 
calculate the probability to have a signal above a given 
background model. The likelihood ratio A is the ratio of 
the probability density for the hypothesis of background 
and signal (Hsig+bkg) over the probability density of only 
background (Hbkg)- 



N 
i=l 



(13) 



TV 
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log 



Pbkgiaiyti) 
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where Usig and N are respectively the unknown number 
of signal events and the total number of events in the 
considered data sample. Psig{ai,t.i) and Pbkg{cti,ii) are 
the probability density function (PDF) for signal and 
background respectively. For a given event i, ti and ai 
represent the time of the event and the angular difference 
between the coordinate of this event and the studied 



The probability densities Pgig and Pbkg are described by 
the product of two components: one for the direction 
and one for the timing. The shape of the time PDF for 
the signal event is extracted directly from the gamma- 
ray light curve assuming the proportionality between the 
gamma-ray and the neutrino fluxes. For the signal event, 
this directional PDF is described by the one dimension 
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point spread function, which is the probability density of 
reconstructing an event at an angular distance from the 
true source position. The directional and time PDF for the 
background are derived from the data using respectively 
the observed declination distribution of selected events in 
the sample and the observed time distribution of all the 
reconstructed muons. Figure |2] shows the time distribution 
of all the reconstructed events and the selected upgoing 
events for this analysis. Once normalized to an integral 
equal to 1, the distribution for all reconstructed events 
is used directly as the time PDF for the background. 
When data is at 0, it means that there are no data taken 
during these periods (ie detector in maintenance) or data 
with a very poor quality (high bioluminescence or bad 
calibration). 
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Fig. 3: Average number of events required for a Scr 
discovery (50 % C.L.) produced in one source located at a 
declination of -40° as a function of the width of the flare 
period. These numbers are compared to the one obtained 
without using the timing information. 
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Fig. 2: Time distribution of the reconstructed events. 
Black: distribution for all reconstructed events. Red filled: 
distribution of selected upgoing events (A > —5.4 and 
/? < 1°). 

The null hypothesis is given with Usig = 0. The obtained 
value of Xdata on the data is then compared to the 
distribution of A given the null hypothesis. Large values of 
Xdata compared to the distribution of A for the background 
only reject the null hypothesis with a confident level equal 
to the fraction of the scrambled trials above Xdata- This 
fraction of trials above Xdata is referred as the p-value. 
The discovery potential is then defined as the average 
number of signal events required to achieve a p-value 
lower than 5ct in 50 % of trials. Figure[3]shows the average 
number of events required for a 5a discovery (50 % 
C.L.) produced in one source located at a declination 
of -40° as a function of the total width of the flare 
periods. These numbers are compared to the one obtained 
without using the timing information. Using the timing 
information yields to an improvement of the discovery 
potential by about a factor 2-3 with respect to a standard 
time-integrated point source search pSl. 

IV. Search for neutrino emission from 

GAMMA-RAY FLARE 

This time-dependent analysis has been applied to bright 
and variable Fermi blazar sources reported in the first 
year Fermi LAT catalogue |9| and in the LBAS catalogue 
(LAT Bright AGN sample |10J). The sources located in the 



visible part of the sky by Antares from which the averaged 
1 day-binned flux in the high state is greater than 20 10~* 
photons. cm^^.s^^ above 300 MeV in the studied time 
period and with a significant time variability are selected. 
This list includes six flat spectrum radio quasars (FSRQ) 
and four BLlacs. Table U lists the characteristics of the ten 
selected sources. 



Name 


OFGL name 


Class 


redshift 


-F300 


PKS0208-512 


J0210.8-5100 


FSRQ 


1.003 


4.43 


AO0235+164 


J0238.6-F1636 


BLLac 


0.940 


13.19 


PKS0454-234 


J457. 1-2325 


FSRQ 


1.003 


13.56 


OJ287 


10855.4-^2009 


BLLac 


0.306 


2.48 


WComae 


J1221.7+28.14 


BLLAc 


0.102 


2.58 


3C273 


11229. U0202 


FSRQ 


0.158 


8.68 


3C279 


J1256. 1-0548 


FSRQ 


0.536 


15.69 


PKS 15 10-089 


J15 12.7-0905 


FSRQ 


0.36 


28.67 


3c454.3 


12254.0-^1609 


FSRQ 


0.859 


24.58 


PKS2155-304 


J2158.8-3014 


BLLac 


0.116 


7.89 



TABLE I: List of bright variable Fermi blazars selected for 
this analysis. F300 is the gamma-ray flux above 300 MeV 
(10~^ photons.cm~^.s~^). 



The light curves published in Fermi web page for the 
monitored sources ifTTI are used for this analysis. These 
light curves correspond to the one-day binned time evo- 
lution of the average gamma-ray flux above a threshold 
of 100 MeV from August 2008 to August 2010. The high 
state periods are defined using a simple and robust method 
based on three main steps. First, the baseline is determined 
with an iterative linear fit. After each fit, the points where 
the flux is above a given threshold are suppressed. When 
the baseline is computed, all the points (green dots) where 
the flux minus its error are above the baseline plus two 
times its fluctuation and the flux is above the baseline 
plus three times its fluctuation are used as priors from 
which the flares are defined. The last step consists on, 
for each selected point, adding the adjacent points for 
which the emission is compatible with the flare. Finally, 
an additional delay of 0.5 day is added before and after 
the flare in order to take into account that the precise 
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time of the flare is not known (1-day binned LC). With 
this definition, a flare has a width of at least two days. 
Figure |4] shows the time distribution of the Fermi LAT 
gamma-ray light curve of 3C454 for almost 2 years of data 
and the determined high state periods (blue histogram). 
With the hypothesis that the neutrino emission follows 
the gamma-ray emission, the signal time PDF is simply 
the normahzed de-noise light curve. 




54700 54800 5491 



55000 55100 55200 55300 

MJD (days) 



Fig. 4: Gamma-ray light curve (red dots) of the blazar 
3C454.3 measured by the LAT instrument onboard the 
Fermi satellite above 100 MeV for almost 2 years of data. 
Blue histogram: high state periods. Green line and dots: 
baseline and significant dots above this baseline used for 
the determination of the flare periods. 



The most significant source is 3C279, which has a pre-trial 
p- value of 1.03 %. The unbinned method finds one high- 
energy neutrino event located at 0.56° from the source 
location during a large flare in November 2008. Figure |5] 
shows the time distribution of the Fermi gamma-ray light 
curve of 3C279 and the time of the coincident neutrino 
event. This event has been reconstructed with 89 hits 
spread on 10 lines with a track fit quality A = —4.4 
and an error estimate (3 = 0.3°. The post-trial probability 
is computed taking into account the ten searches. The 
final probability, 10 % is compatible with background 
fluctuations. 




Fig. 5: Gamma-ray light curve (black dots) of the blazar 
3C279 measured by the LAT instrument onboard the 
Fermi satellite above 100 MeV. Blue histogram: high 
state periods. Green dashed line: fit of a baseline. Red 
histogram: time of the ANTARES neutrino event in coin- 
cidence with 3C279. 



V. Summary 

This paper discusses the first time-dependent search for 
cosmic neutrinos using the data taken with the full 12 
lines ANTARES detector during the last four months 
of 2008. Time-dependent searches are significantly more 
sensitive than standard point-source search to variable 
sources thanks to the large reduction of the background of 
atmospheric muons and neutrinos over short time scales. 
This search has been applied to ten very bright and vari- 
able Fermi LAT blazars. The most significant observation 
of a flare is 3C279 with a p-value of about 10 % after 
trials for which one neutrino event has been detected 
in time/space coincidence with the gamma-ray emission. 
Limits have been obtained on the neutrino fluence for the 
ten selected sources. The most recent measurements of 
Fermi in 2009-11 show very large flares yielding a more 
promising search of neutrinos |12|. 
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Abstract. The ANTARES telescope is well suited to 
detect neutrinos produced in astrophysical transient 
sources as it can observe a full hemisphere of the 
sky at all the times with a duty cycle close to unity 
and an angular resolution better than 0.5 degrees. 
Potential sources include gamma-ray bursts (GRBs), 
core collapse supernovae (SNe), and flaring active 
galactid nuclei (AGNs). To enhance the sensitivity of 
ANTARES to such sources, a new detection method 
based on coincident observations of neutrinos and 
optical signals has been developed. A fast online muon 
track reconstruction is used to trigger a network of 
small automatic optical telescopes. Such alerts are 
generated one or two times per month for special 
events such as two or more neutrinos coincident in 
time and direction or single neutrinos of very high 
energy. Since February 2009, ANTARES has sent 37 
alert triggers to the TAROT and ROTSE telescope 
networks, 27 of them have been followed. First results 
on the optical images analysis to search for GRB and 
core-collapse SNe will be presented. 

I. Introduction 

The detection of high energy cosmic neutrino from a 
source would be a direct evidence of the presence of 
hadronic acceleration within the source and provide im- 
portant information on the origin of the high energy 
cosmic rays. Powerful sources of transient nature, such 
as gamma ray bursts or core collapse supernovae, offer 
one of the most promising perspectives for the detection 
of cosmic neutrinos as, due to their short duration, the 
background from atmospheric neutrinos and muons is 
strongly reduced. For example, several authors predict the 
emission of neutrinos in correlation with multi-wavelength 



signals, e.g. the Fireball model of GRBs As neutrino 
telescopes observe a full hemisphere of the sky (even 
the whole sky if downgoing events are considered) at all 
times, they are particularly well suited for the detection 
of transient phenomena. 

The ANTARES neutrino telescope JZl is located in the 
Mediterranean sea, 40 km South of the french coast of 
Toulon, at a depth of about 2500 m below sea level. 
The detector is an array of photomultipliers tubes (PMTs) 
arranged on 12 slender vertical detection lines. Each string 
comprises up to 25 floors, i.e. triplets of optical modules 
(OMs) housing one PMT each. Data taking started in 
2006 with the operation of the first line of the detector. 
The construction of the 12 line detector was completed 
in 2008. The main goal of the experiment is to search 
for neutrinos of astrophysical origin by detecting high 
energy muons (>100 GeV) induced by their neutrino 
charged current interaction in the vicinity of the detector. 
Due to the large background from downgoing cosmic ray 
induced muons, the detector is optimised for the detection 
of upgoing neutrino induced muon tracks. 

In this paper, the implementation and the first results of a 
strategy for the detection of transient sources is presented. 
This method, earlier proposed in ||3|, is based on the 
optical follow-up of selected neutrino events very shortly 
after their detection by the ANTARES neutrino telescope. 
The alert system is known as "TAToO" (Telescopes and 
ANTARES Target of Opportunity) ||41, 

II. ANTARES NEUTRINO ALERTS 

The criteria for the TAToO trigger are based on the 
features of the neutrino signal produced by the expected 
sources. Several models predict the production of high 
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energy neutrinos greater than 1 TeV from GRBs f5\ 
and from Core Collapse Supernovae [6J. Under certain 
conditions, multiplet of neutrinos can be expected Q. 

Two online neutrino trigger criteria are currently imple- 
mented in the TAToO alert system: 

• the detection of at least two neutrino induced muons 
coming from similar directions within a predefined 
time window; 

• the detection of a single high energy neutrino induced 
muon. 

A basic requirement for the coincident observation of a 
neutrino and an optical counterpart is that the pointing 
accuracy of the neutrino telescope should be at least 
comparable to the field of view of the TAROT (S) and 
ROTSE L9J telescopes (« 2° x 2°). 

To select the events which might trigger an alert, a fast and 
robust algorithm is used to reconstruct the calibrated data. 
This algorithm uses an idealized detector geometry and is 
independent of the dynamical positioning calibration. As 
a result, the hits of the three OMs of a storey are grouped 
and their location assigned to the barycenter of the storey. 
The storey orientations as well as the line-shape deviations 
from straight lines are not considered in the online recon- 
struction. A detailed description of this algorithm and its 
performances can be found in ifTOl . The principle is to 
minimize a which compares the times of selected hits 
with the expectation from a Cherenkov signal of a muon 
track. The resulting direction of the reconstructed muon 
track is available within about 10 ms and the obtained 
minimal is used as fit quality parameter to remove 
miss-reconstructed tracks. 

Atmospheric muons, whose abundance at the ANTARES 
detector fTTl is roughly six orders of magnitude larger 
than the one of muons induced by atmospheric neutrinos, 
are the main background for the alerts and have to 
be efficiently suppressed. Among the surviving events, 
neutrino candidates with an increased probability to be 
of cosmic origin are selected 1121 . 

In order to establish the criteria used for our neutrino 
selection, we have analysed a subsample of data taken by 
ANTARES after the completion of the 12-line detector, 
corresponding to a livetime of 70.3 days. During this 
period, around 350 upgoing neutrino candidates were 
reconstructed and have been compared to a Monte Carlo 
(MC) simulation of atmospheric muons and neutrinos 
using the same lifetime. 

Figure [T] shows the angular resolution of the online algo- 
rithm as a function of the neutrino energy. This resolution 
is defined as the median of the space angular difference 
between the direction of the incoming neutrino and the 
reconstructed neutrino-induced muon. For neutrinos with 
an energy higher than a few tens of TeV, an angular 
resolution of 0.4 degree is achieved, despite of the approx- 
imations related to the detector geometry. For example, the 
inclination of the ANTARES line for a typical sea current 



of 5 cm/s induces a systematic angular deviation of less 
than 0.2 degree. 




5 5.5 6 
log, JE^, (GeV)] 

Fig. 1: Angular resolution obtained for both online and 
offline reconstructions as a function of the neutrino energy. 



A. Multi-neutrino trigger 

The typical signature of the transient emission of high 
energy neutrinos is a neutrino burst, i.e. a multiplet of 
neutrino events originating from the source in a short time 
window. A trigger for this event type is implemented as 
the detection of two upgoing events reconstructed with 
at least two lines in a 15 minutes time window with 
a maximum angular difference of 3°. The time window 
was optimized to include most predictions of the neutrino 
emission by various models for transient sources. The 3° 
angular window was selected to match the convolution of 
the track reconstruction angular resolution and the field 
of view of the robotic optical telescopes (« 2° x 2°). 
The accidental coincidence rate due to background events, 
from two uncorrelated upgoing atmospheric neutrinos, is 
estimated to be 7 x lO^'^ coincidences per year with the 
full ANTARES detector. With such a small background, 
the detection of a doublet (triplet) in ANTARES would 
have a significance of about 3 (5) sigma. 

B. High energy event trigger 

Since the neutrino energy spectrum for signal events is 
expected to be harder than for atmospheric neutrinos, a 
cut on the reconstructed energy efficiently reduces the 
atmospheric neutrino background while most of the signal 
events are kept. The selection of the alert candidates is 
based on two simple energy estimators: the number of 
storeys used in the track fit and the total amplitude (in 
photoelectrons) of the hits in the storeys. 

The event selection for the high energy trigger has been 
tuned on atmospheric neutrinos in order to obtain a false 
alarm rate of about 25 alerts per year This rate was agreed 
between ANTARES and the optical telescope collabora- 
tions. A requirement of at least 20 storeys on at least three 
lines and an amplitude greater than 180 photoelectrons 
will select around 25 high energy events per year with 
the full 12 line configuration of the ANTARES detector. 
The TAToO alert criteria select neutrinos of energies above 
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10 TeV for the single high energy trigger (calculated with 
a neutrino Monte Carlo generated with an E^^ energy 
spectrum). 

Figure [2] shows the estimation of the point spread function 
for a typical high energy neutrino alert. Around 70% of the 
events are contained in the field of view of a typical robot- 
ical telescope (« 2° x 2°). With a larger delay (few tens 
of minutes after the time of the burst), we are able to run 
the standard reconstruction tool which provides a much 
better angular resolution using the dynamical positioning 
of the detector lines |13|. Simulations indicate that, with 
this algorithm, ANTARES reaches an angular resolution 
smaller than about 0.3-0.4° for neutrino energies above 
10 TeV (curve labeled 'offline algorithm' in Figure [TJ. 




Fig. 2: Bi-dimensional angular resolution. The black 
square corresponds to the TAROT telescope field of view 
(w 2° X 2°). 



III. Observation strategy of the robotical 

TELESCOPES 

ANTARES is organizing a follow-up program in col- 
laboration with the TAROT and ROTSE telescopes The 
TAROT |8| network is composed of two 25 cm optical 
robotic telescopes located at Calern (France) and La Silla 
(Chile). The ROTSE |91 network is composed of four 
45 cm optical robotic telescopes located at Coonabarabran 
(Australia), Fort Davis (USA), Windhoek (Namibia) and 
Antalya (Turkey). The main advantages of these instru- 
ments are the large field of view of about 2x2 square 
degrees and their very fast positioning time (less than 10s). 
These telescopes are perfectly tailored for such a program. 
Thanks to the location of the ANTARES telescope in the 
Northern hemisphere (42.79 degrees latitude), all the six 
telescopes are used for the optical follow-up program. 
Depending on the neutrino trigger settings, the alert are 
sent at a rate of abour one or two times per month. With 
the current settings, the connected telescopes can start 



taking images with a latency of the order of one minute 
with respect to the neutrino event (TO). 

As it was said before, the rolling search method is 
sensitive to all transient sources producing high energy 
neutrinos. For example, a GRB afterglow requires a very 
fast observation strategy in contrary to a core collapse 
supernovae for which the optical signal will appear several 
days after the neutrino signal. To be sensitive to all 
these astrophysical sources, the observational strategy is 
composed of a real time observation followed by few 
observations during the following month. For the prompt 
observation, 6 images with an exposure of 3 minutes 
and 30 images with an exposure of 1 min are taken 
respectively by the first available TAROT and ROTSE 
telescopes. The integrated time has been defined in order 
to reach an average magnitude of about 19. For each 
delayed observation, six images are taken at T0h-1,h-2, +3, 
+4, +5, +6, +7, +9, +15, +21 days after the trigger for 
TAROT (8 images for ROTSE the same days plus TOh-16 
and TOh-28 days). 



IV. Optical image analysis 

Once the images are taken, they are automatically dark 
subtracted and flat-fielded at the telescope site. Once the 
data are copied from the telescopes, an offline analysis is 
performed combining the images from all sites. This off- 
line program is composed by three main steps: astrometric 
and photometric calibration, subtraction between each 
image and a reference one and light curve determination 
for each variable candidates. 

Curently, two offline analysis pipelines are used: the 
ROTSE automated pipeline and one specially adapted 
to the TAROT and ROTSE image quality based on a 
program originally developed for the supernova search 
in the SuperNova Legacy Survey (SNLS) project. Cases 
like variable PSF due to the atmospheric conditions or 
the lower quality images on the CCD edges have to be 
optimized in order not to loose any optical information. 
The choice of the reference is based on quality criteria 
such as the limiting magnitude and the seeing. For the 
GRB search, the reference is picked among the follow- 
up observations (few days after the alert) where no GRB 
signal is expected anymore while for SN search, either 
we consider the first night observation or we order it few 
months later to have a better quality reference in absence 
of a SN signal. It is also planned that the image analysis 
step will be included at the end of the automatic detection 
chain. 

The ROTSE pipeline has been applied to five alerts from 
which optical images have been recorded during the first 
24 hours after the neutrino alert sending. The minimum 
delay between the neutrino detection and the first image is 
around 70 s. No object has been found for which the light 
curve is compatible with a fast time decreasing signal. 
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V. Conclusion 

The method used by the ANTARES collaboration to 
implement the search for coincidence between high energy 
neutrinos and transient sources followed by small robotic 
telescopes has been presented. Of particular importance 
for this alert system are the ability to reconstruct online the 
neutrino direction and to reject efficiently the background. 
With the described ANTARES alert sending capability, the 
connected optical telescopes can start taking images with a 
latency of the order of one minute, which will be reduced 
to about 15 s in the near future. The precision of the 
direction of the alert is much better than one degree. The 
quasi-online availability of a refined direction obtained 
using the measured geometry of the ANTARES detector 
further improves the quahty and efficiency of the alert 
system. 

The alert system is operational since February 2009, and 
as of May 2011, 37 alerts have been sent, all of them 
triggered by the high energy selection criterium. No dou- 
blet trigger has been recorded yet. After a commissioning 
phase in 2009, almost all alerts had an optical follow-up 
in 2010, and the live time of the system over this year 
is strictly equal to the one of the ANTARES telescope, 
namely 87%. These numbers are consistent with the 
expected trigger rate, after accounting for the duty cycle of 
the neutrino telescope. The image analysis of the 'prompt' 
images has not permitted to discover a GRB afterglow 
associated to the high energy neutrino. The analysis of 
the rest of the images to look for the Ught curve of a core 
collapse SN is still on-going. 

The optical follow-up of neutrino events significantly 
improves the perspective for the detection of transient 
sources. A confirmation by an optical telescope of a 
neutrino alert will not only provide information on the 
nature of the source but also improve the precision of 
the source direction determination in order to trigger 
other observatories (for example very large telescopes 
for redshift measurement). The program for the follow- 
up of ANTARES neutrino events is already operational 
with the TAROT and ROTSE telescopes and results based 
on analysis of the optical images will be presented in a 
forthcoming paper. This technique could be extended to 
observations in other wavelength regimes such as X-ray 
or radio. 
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Abstract. Neutrinos of all flavours are produced during 
a Supernova (SN) explosion. The antineutrinos can 
interact in the medium surrounding the underwater 
ANTARES telescope and produce positrons emitting 
Cerenkov radiation which is recorded by the optical 
modules (OMs) of the detector. The signature of this 
event is a simultaneous increase of the counting rate 
in the detector. The aim of this work is to develop 
several methods to detect the SN signal in ANTARES 
detector. In particular this analysis, using the specific 
geometrical configuration of the OMs, minimizes the 
noise effect introduced by the bioluminescence during 
the measurement. The significance as a function of the 
number of active OMs and of the distance to SN is 
also evaluated. 

I. Introduction. 

The ANTARES neutrino telescope HI is a large water 
Cerenkov detector that could provide the opportunity to 
detect the Cerenkov light produced by positrons from SN 
antineutrinos via reaction v,, + p ^ n + . 

To investigate this possibility we evaluate the response 
of the ANTARES optical modules (OM) to SN neutrinos 
with a help of Geant4 simulations and compare the 
possible SN signal rates to the ANTARES OMs rates, 
which are recorded during the normal run operations. 

The ANTARES detector is made of detection lines that 
consist of 25 storeys, a group of a 3 OMs that are ar- 
ranged close to each other (distance is less than Im). 
This configuration provides an opportunity to reduce the 
background using only hits that are in coincidence be- 
tween OMs in the same storey. 

The methods to minimize the effect of the biolumines- 
cence burst are discussed in the Section 3 while the 
statistical analysis based on singles rates, double coinci- 
dences and triple coincidences in one storey is presented in 
Section 4, 5 and 6. The performances of the three different 
meth-ods are also compared, and significance is evaluated. 



II. Simulations. 

We performed the simulations using Geant4 with the 
detailed OM description. In these simulations we included 
the full 3 OMs ANTARES storey configuration while 
for the positrons we used the energy spectrum and the 
flux obtained averaging in time the distributions from the 
model 57 of A.Burrows f2l for a SN1987A like event at 
a 10 kpc distance. 

As a result we have obtained an excess of 14 hits in 
a 105ms time interval on each OM. The time interval 
of 105ms was chosen because it both corresponds to 
the standard time slice of the data flow in ANTARES 
detector and to the period of SN explosion where the 
highest neutrino flux is produced according to f2\. To 
fully exploit the geometry of ANTARES detector, we have 
also obtained the number of coincidences produced by 
the positrons light between two or three OMs of the same 
storey in a 25ns time window. The possibility to use these 
multi coincidences will be described in the next sections. 
To evaluate the number of hits from an event at a distance 
R different from the 10 kpc we used a simple proportion, 
knowing that $ * AuR — const . 

III. The bioluminescence filter. 

In addition to the Cerenkov light produced by muons, 
neutrino telescopes also detect Cerenkov light from the 
decay of radioactive elements and light from luminescent 
organisms. Fig[T] represents the time variation of the 
optical background rate detected in one OM. One of the 
main contributions comes from the radioactive decay of 
""^K. The frequency of this noise is evaluated by MC 
calculations to be around 40kHz and does not change 
since the concentration depends on the salinity that is 
almost constant on time and sea depth. However in fig{T] 
the mean rate is about 80kHz, the additional contribution 
being due to a constant bioluminescent activity. 

Above the slow variable noise one can see high sharp 
peaks due to the emission of light bursts from macro 
organisms in transit close to the OM. Later in this note 



V.KULIKOVSKIY et al. SN NEUTRINO DETECTION IN THE ANTARES 



iiiO 



1^ 




Fig. 1: Example of singles rate (0.3pe threshold) events 
in one OM of the detector. 



we will term these peaks as "bioluminescent burst" to 
distinguish them from the continuous or slowly variable 
bioluminescence noise. 

To exclude the bioluminescence bursts on one OM we 
have fitted the distribution of the number of hits in 
a 105ms period. The fit is performed with a Poisson 
function. There are three parameters in the fit: A is the 
mean value of the Poisson distribution, j™"^ is the right 
margin of the fit and A is the scaling factor (height of the 
distribution). 

The filter for the bursts is very simple: if in one time 
frame the number of hits for a particular OM is higher 
than j™"^, then the OM is considered as affected by a 
bioluminescence burst and for this time frame OMi is not 
used in the calculation of the total detector hit counts. 

IV. 4 Supernova signal with the single OMs 

RATES. 

Detection significance can be defined as the ratio between 
the measured signal and the uncertainty of the measure- 
ment 

S = ^^g""^ (15) 

^Tneasurement 

In case of 14 signal hits (hgignai) and of 80kHz of back- 
ground for 900 active OMs (Nqm) the significance can 
be estimated by the formula: 

hsignalNoM 
background NoM 

The main idea of all methods described here being to mea- 
sure the increase of the total hits in the detector in one time 
slice and to compare it with the previous measurements. 
For single hits in the specific case of ANTARES detector, 
this cannot be directly applied due to the bioluminescence 
bursts that have been therefore excluded using the filter 
as described above. The procedure we have used is the 
following: 



5* = 



(16) 



a) The data collected in a 45min time interval, sufficiently 
large to have good statistics, were used to determine the 
fit h^^^ of the hit counts distribution for each z*'* OM in 
the detector, h^^^ can be seen as the probability density 
function (PDF) of the hits in OM^. The results from the 
fit were also used to define quality cuts. In particular 
only OMs for which the hit distribution follows the next 
criteria where accepted for the analysis: the fit converges, 
A >2000, the area of the Poisson component is larger than 



30% of the total area (in other words: there is less than 
70% of the bio-burst hits). 

b) The data are then analysed time slice by time slice. At 
each time slice the hits of all the OMs that have passed 
the bioluminescence cut are summed together to determine 
the total hit count of the detector in that particular time 
slice. 

c) The total hit count of the detector should be described 
by the "detector PDF" distribution, which is different for 
every time slice because OMs which have passed the 
bioluminescence cut are different. The detector PDF dis- 
tribution can be evaluated using the PDFs for every OM, 
but this is a very time consuming procedure. To simplify 
we have made a simulation using all the parameters from 
fits of the single PDFs which proved that detector PDF 
is still a Gaussian with an average M = J^rrii where 

is mean of PDF for i*'' working OM and with a sigma 
S = \/^crf where ai is a variance for i*'* PDF. 

To estimate the significance two runs were selected. One 
of them was #39656, recorded on 16 May 2009, where bi- 
oluminescence conditions were good (rate: 80kHz for the 
bottom, 60kHz - middle, 60kHz - top), another one was 
#40154 recorded on 13 April 2009 when bioluminescence 
was higher (rate 200-160-140kHz). 

The parameters of the fit (A, A, j™"^) were then extracted 
from the experimental distributions as deduced from the 
analysis of the 2 runs. In order to compare with the other 
methods described in the note, we have calculated the 
significance assuming 900 working PMTs and keeping the 
same mean characteristics of the background of the two 
previous runs. 

The results of these simulations are shown in fig|2] where 
the significance is plotted as a function of the distance 
of the SN. 900 OMs in the same conditions of the runs 
39656 and 40154 were simulated. The decrease of the 
significance with the increase of the SN signal (from 
distance 3.5 kpc and closer) is due to the bioluminescence 
filter: when the SN signal is very high a lot of OMs are 
excluded from the analysis. The dotted lines with a l.Vcr 
label in the same figure are obtained when the effect due 
to the long-term fluctuations is included in the calculation. 
This effect is probably ascribed to the baseline rate which 
is slightly changing during the analysis time of 45min and 
which increases the width of the distribution by a factor 
of about 1.7. 

A maximum significance of 5 for SN at the 4 kpc is 
obtained at low bioluminescence level (run #39656) for 
the full working detector when the long-term fluctuation 
is included. 



V. Supernova signal for double coincidences 

BETWEEN OMs 

The ANTARES configuration allows the light from ''''K 
decay to be simultaneously seen by a couple of OMs in a 
storey. This fact is rather well known and it is used in the 
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#39656 B0-G0-6OkHi 

a396S6SO-60-60kHi 1.7fl 

(l401S4 2[l0.1fi0.140kHi 

■ — #40154200-160-140kHi 




a 10 
distance(kpc] 

Fig. 2: Simulated significance for 900 OMs detector in 
dependence of the distance from SN explosion. OMs 
signal is simulated according to run #39656(solid) and 
#40154(dotted) conditions. Dashed and dot dashed lines 
represent calculations with the 1.7 increase of sigma to 
account for long term fluctuations. 
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RMS 

■»} J ndf 
Prob 



22ai79 
Q.01244 
11.26 
1155/137 


4.2246+04 ±1.4-056^02 
2955 ± 109.8 
-0.5992 < 0.1213 
4.224 ±0.149 




Fig. 3: Example of the time coincidence distribution for 
storey 1635, OM 1 and 2. Dashed line represents random 
coincidence level. The area of the Gaussian ac-counts for 
the number of hits from ^^K coincidences. 
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Fig. 4: Uncertainty on the coincidence rate extracted from 
the fit of the simulated distribution for different time sUce 
numbers (solid curve). Dotted curve represents error in 
absence of the random background. 



couple of OMs assuming an 80kHz random background 
on each OM and a Gaussian distribution with 16Hz for 
the coincidences. These simulations gave the uncertainty 
on the extracted true coincidence rate as reported in fig|4] 
for various number of time slices Nts- 
In this method we calculate the global rate of the coin- 
cidences in the detector at any time slice and compare 
it with the sum of the rates of each couples. The global 
detector rate was determined from the time differences 
between all coincidence hits of every OM couple that were 
collected in a common distribution. This distribution was 
then fitted with the same procedure described above and 
the error evaluated as in fig |4] where the number of couples 
instead of the number of time slices was assumed. The true 
coincidence detector rate extracted from the fit is stable in 
time and equal to the sum of the every OM couple rates 
in case of normal work or higher in case of SN explosion. 



detector calibrations [3]. Later on by the couple of OMs 
we assume only OMs from the same storey. 

For every couple of OMs the plot of the time difference 
between every two hits from both OMs can be deter- 
mined. One example of this time distribution is shown 
in figj3] 

We can distinguish a plateau due to the random coinci- 
dences in the time window r of 20ns and a Gaussian 
shaped distribution due to "'"K decays seen by the both 
modules (later we call them true coincidences to distin- 
guish from the random coincidences). To extract these two 
components the distribution was fitted with a Gaussian 
plus a constant. The rate of the true coincidences extracted 
from the fit is about 16Hz for a couple of OMs and it is 
independent from the bioluminescence activity. The rate 
is thus stable in time, but it depends from the couple of 
OMs since efficiency of the OMs is different. 

To understand how the uncertainty on the true coincidence 
rate from the fit depends on the statistics i.e. on the 
number of the summed time slices, we performed a simple 
simulation of the time difference distribution for one 



We verified this procedure on real data. To this purpose 
one period with a 450 OM couples was analysed and 
the rate for every couple was calculated from 20000 
time slices of the data to provide as seen from fig|4] 
an uncertainty less than 5%. The global detector true 
coincidence rate was then obtained from the fit for every 
time slice. The difference between this rate and the sum of 
the true coincidence rates of every couple was evaluated 
and normalized to the detector rate. As expected, the 
distribution of this normalized difference has a mean 
and sigma 17%, in good correspondence with the expected 
error from fig|4]for 450 Nts- 

Simulations with Geant4 estimates a value of 22Hz for the 
true coincidence rate, the difference with respect to the 
measured 16Hz being ascribed mostly to the description 
of the lateral region of the OM. Geant4 simulations also 
show that the light produced by the positrons induced 
by SN neutrinos may simultaneously hit two and even 
three OMs in the same storey. These additional events 
increase the coincidence rate during the SN explosion. 
The Geant4 simulation of this effect gives a rate of about 
2.8Hz in the first 100ms time window i.e. an increase of 



V.KULIKOVSKIY et al. SN NEUTRINO DETECTION IN THE ANTARES 



5.5 

5 
4.5 
4 




\. . . 














\ 
























































3.5 
















3 
2.5 
















2 












-. 






- 




1 1 1 1 


1 . 1 1 










5 


1DQ 150 200 250 300 




200 250 300 

Background rate(kHz) 



Fig. 5: Simulated 5a significance as a function of SN 
distance and background rate for the double (up) and a 
triple coincidence (bottom) methods. 



0.28 events for each pair of OMs in the considered 105 
ms time slice for a 10 kpc SN1987A. The measurement 
of this small increase, which was reduced down to 2Hz to 
take into account the difference between MC and data, can 
extend the significance to the SN explosion in ANTARES 
detector 

The 5(7 significance obtained from the simulations of 900 
OMs detector is reported in fig|5]as a function of the SN 
distance for various bioluminescence rates. For an event 
occurring at 4 kpc the significance is comparable to the 
previous method for a background rate of the lOOkHz. 

Given the ANTARES geometry, a triple coincidence in 
each storey can occur whenever a single event (light from 
a track, a decay) is simultaneously seen by the 3 OMs. 
In our approach we require that 3 hits, one for each OM, 
are detected during a 20ns time window. This means also 
that in every pair of the hits in the triple coincidence, the 
time difference is less than 20ns. In order to obtain the 
time difference distribution like for doubles it was decided 
to use three dependent distributions: ti — t-z, ^2 — ^3 and 
ts — 12, where ti is time of hits from OMi. In this case the 
random coincidences do not form a plateau in the time 
difference distribution. This is because, for example, the 
difference <i — 12 is constrained by the condition that both 
ti and t2 must be close to <3 by less than 20ns. Using 
Bayes' theorem one can evaluate the shape of the random 
coincidences part: 



The fit is thus done adding the Gaussian of the true coin- 
cidences to the random part (4). True coincidences rate 
for the storey is calculated as a mean value of the true 
coincidence rates extracted from the fit of the three dis- 
tributions ti — t2, t2 — ts and ^3 — t2- 
To find the true triple coincidence rate, one month of 
data was analysed. It was found that the mean storey 
rate is about 0.065Hz and it is stable during 3 months. 
Simulations with Geant4 give a higher value: 0.2Hz. This 
large discrepancy can be explained considering the akeady 
observed behaviour in the double coincidence case where 
the simulations give a value around 22Hz against the 
observed 16Hz. In case of a triple coincidence this inef- 
ficiency is additionally amplified. The Geant4 simulation 
also estimates the total number of triple coincidences from 
SN to about 0.035 per storey in a time slice. Assuming 
the same damping in the detection efficiency as for '"'K 
rates (0.2Hz MC against 0.065Hz experimental), the SN 
triple coincidence rate is expected to be about 0.1 IHz. 

Similarly to the double, the coincidences rate of the triple 
coincidences in the detector during one time slice is cal- 
culated. Unfortunately the number of triples is much less 
than the doubles, the fitting procedure is not applicable 
so that the significance is affected by the random coinci- 
dences ([TtJi. The simulated significance as a function of 
the SN distance and background is shown in the fig|5] 

VI. Conclusion. 

We have investigated different possibilities for detecting 
SN neutrinos with the ANTARES telescope. The signifi- 
cances estimated for the three methods we have described 
are comparable (fig|2j fig|5]l. Double and triple coin- 
cidences seem preferable because significance increase 
with the SN signal and they do not strongly suffer from 
the bioluminescence activity. Of course the simultaneous 
acquisition of an external SN trigger like the one provided 
by the SNEWS collaboration will strength our analysis. 
The study to implement this possibility is now in progress. 
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Abstract. According to the theory of hadronic acceler- 
ation, ultra high energy cosmic rays are expected to 
be accompanied by gamma-rays and neutrinos from 
pion decays formed in the interactions of protons 
with photons. While gamma-rays have been linked to 
astrophysical sources by many experiments (H.E.S.S, 
MAGIC, Fermi), no point source of UHECRs or neu- 
trinos have been found so far. In this contribution, we 
present the results of multimessenger stacking sources 
analysis developed to investigate the correlation of 
arrival directions of neutrino candidate events and 
UHECRs. This analysis has been applied on neu- 
trino candidate events detected during 2007-08 by the 
ANTARES telescope (ANTARES collaboration 2011), 
and 69 ultra-high energy cosmic rays observed by the 
PAO (The Pierre Auger Collaboration 2010). 

I. Introduction 

Finding point sources of both UHECRs and neutrinos is 
a very challenging task, as both messengers have many 
down-sides. UHECRs are rare and do not point back 
to their sources, since they are scrambled by galactic 
and intergalactic magnetic fields lUl, EJ, S, [4}, |5|, 

El, Q, El, 0, nni, im, hh, m, mi, na, m- 

Also, due to the interaction with the cosmic microwave 
background photons, their range may be limited to the 
distance of about lOOMpc or less HT), (H, ||T9l, EOl . 
li2Tl . However, they are detectable with large shower 
arrays, like the Pierre Auger Observatory (PAO), which so 
far reported the observation of few tens of events above 
55EeV 1211, |l22|. On the other side, cosmic neutrinos, as 
neutral and weak interacting particles, should point back 
to their sources and their traveling distances should not 
be limited, but at the same time they are very difficult to 
detect. Currently operating neutrino telescopes ANTARES 
and IceCube have not yet observed excess above the 
atmospheric neutrino flux coming from air showers |23|, 
I24J. 



Previously, the Pierre Auger Observatory (PAO) reported 
an anisotropy in the arrival directions of UHECRs ETIl and 
indicated a correlation with Active Galactic Nuclei (AGN) 



from the VCV catalog (25). The correlation was the most 
significant for 27 cosmic rays with energies higher than 
57EeV and AGNs at distances less than 75Mpc. The 
suggested correlation with the nearby AGN sources mostly 
following the location of the supergalactic plane decreased 
in the subsequent analysis ll22]| with 69 events at energies 
above 55EeV, observed until 31st December 2009. 

In this paper, we investigate the correlation of arrival 
directions of 2190 neutrino candidate events detected by 
5-12 Une ANTARES neutrino telescope, and 69 UHECRs 
observed by the PAO. This stacking sources analysis is 
developed for a blinded set of neutrino events. Blinding 
was performed by scrambling in right ascension. At the 
end, the analysis was done with the unblinded neutrino 
dataset. 

II. The ANTARES telescope and data samples 

The ANTARES neutrino telescope is located in the 
Mediterranean Sea, about 40km off the southern coast of 
France (42 48N, 6 lOE), at a depth of 2475m. It was 
completed in 2008, and its final configuration is a three- 
dimensional array of 885 photomultipliers in glass spheres 
(optical modules), distributed along twelve lines. These 
photomultipliers detect Cherenkov photons, from relativis- 
tic muons produced in neutrino interactions nearby the 
detector. The total instrumented volume of the ANTARES 
telescope is about lO^m'^. The data acquisition system of 
the detector is based on "all-data-to-shore" concept, in 
which signals from the photomultipliers above a given 
threshold are digitized and sent to shore for processing. 
The data used in this analysis were collected during 2007 
and 2008, while the detector was operating with 5 to 12 
lines. For a part of that period, data acquisition was inter- 
rupted for the employment of new lines, and in addition, 
some periods were excluded due to high bioluminescence- 
induced optical background. The resulting effective live 
time is 304 days. The final data sample consists of 2190 
up-going (zenith angle up to 90°) neutrino candidate 
events. No selection was done based on the energy re- 
construction. The angular resolution was estimated to be 
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0.5±0.1 degrees. Details of data reconstruction, effective 
area, angular resolution and the dataset are given in [24J. 

This sample of neutrino candidate events was correlated 
with UHECR events recorded with the PAO surface detec- 
tors between January 1st 2004, and 31st December 2009. 
Those events have zenith angles below 60 degrees, and 
reconstructed energy above 55EeV. 69 events satisfy these 
cuts, and they are all in the ANTARES telescope field of 
view, as the field of view for the ANTARES telescope and 
the PAO greatly overlap. 

III. Background simulations 

A Monte Carlo set with 10^ skies is generated, each with 
2190 neutrinos and 69 cosmic rays. In each of the million 
skies the position of the UHECRs is fixed according to 
the PAO dataset, and the positions of 2190 neutrinos are 
obtained from the blinded 5-12 line dataset scrambled in 
right ascension. The numbers of neutrinos within bins of 
any size, centered on 69 UHECR events can be counted, 
probability density functions can be calculated and fitted 
with Gaussian distributions, giving the mean neutrino 
count expected from the randomized background samples, 
and the standard deviation of the neutrino count. 

Figure [T] shows an example count of neutrinos within 
bins of 1-10 degree size. The count of events is done by 
adding neutrinos for which the minimum angular distance 
to UHECRs is smaller than the bin size. In this way, 
a double counting of neutrino events is avoided. After 
optimizing the bin, the significance of observed number of 
neutrino events within 69 bins can be calculated by com- 
parison with the distributions for the pure background MC 
samples. After optimizing the bin, and unblinding the data, 
the significance of observed number of neutrino events 
within 69 bins can be calculated by comparison with the 
distributions for the pure background MC samples. 

IV. Signal simulations 

Simulation of sources with E^'^ spectrum and equal 
flux strength from each 69 UHECR direction was per- 
formed. Flux values from O.SxlO^^GeV cm^^ s^^ to 
10xlO~*GeV cm~^ s^^ are considered. Further, for each 
flux, the event rate per source is determined, using the 
effective area for 5-12 lines and the live time given in 
[24]. For every source, the amount of signal neutrinos is 
determined according to the Poisson distribution with the 
event rate per source R as mean value. 

^.lO^GeV 

R = t <^{E,)A,ff{E,)dE, 

JlO^GeV 

where i is a live time of the ANTARES detector, $ is the 
flux from each individual UHECR direction, and A^fj is 
the ANTARES telescope effective area. 

Signal neutrinos are randomly generated following a 
Gaussian distribution centered on the measured UHECR 
coordinates, and the ANTARES telescope resolution [24 1. 



exp {-{xi{a, S) - fi)"^ / {2<tI)) X2{a, 6) 
exp {-{x2{a, S) - ^ifl {2al)) Xf^^ia, 6) 

where Xi = Xi (a, S) are neutrino coordinates, Xfin{a,6) 
are the final signal neutrino coordinates, ai is the tolerance 
value for the magnetic deflection, and (T2 is the ANTARES 
telescope angular resolution. The width of the tolerance 
value accounting for a magnetic deflection Gaussian is 
chosen to be 3 degrees, based on the PAO results [22,25]. 

The same amount of background neutrinos is removed 
from a declination band of 10 degrees centered on each 
UHECR to ensure that every random sky has 2190 events, 
and to keep the neutrino declination distribution profile 
close to the observed profile. The angular resolution 
of PAO air shower reconstruction is about 0.9 degrees, 
less than the expected magnetic deflection, so it is not 
considered additionally. 

V. Sensitivity and bin optimization 

Although it is not possible to know the actual upper limit 
that will result from an experiment until looking into 
unblinded data, Monte Carlo predictions can be used to 
calculate the average upper limit, or so called Feldman- 
Cousins sensitivity lf26l . Il27l . that would be observed after 
hypothetical repetition of the experiment with expected 
background n^, and no true signal Ug = 0. Over an en- 
semble of experiments with no true signal, the background 
nfc will fluctuate to different Uobs values with different 
Poisson probabilities, and upper limits /igo. 

The "mean upper limit" is the sum of these expected 
upper limits, weighted by their Poisson probabihty of 
occurrence: 



"o6s=0 

Over an ensemble of identical experiments, the strongest 
constraint on the expected signal flux corresponds to a 
set of cuts that minimizes the so called "model rejection 
factor" /igo/ris and at the same time minimizes the 
mean flux upper limit that would be obtained over the 
hypothetical experimental ensemble. 
The described Feldman-Cousin's approach with the Rolke 
extension [28 1 was used to calculate the mean upper limit 
on E^^ flux per source, for a 90% confidence level, from 
10^ background samples, as a function of a search bin, as 
shown on a Figure |2] 

The search bin that minimizes the mean upper limit 
for 3 degrees tolerance value accounting for a magnetic 
deflection is found to be 4.9 degrees. 

VI. Discovery potential 

With the bin size optimized and fixed, it is possible to 
estimate the probability of making a 3cr or a 5(7 discovery 
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given a certain signal flux. This is done by taking the 
chosen significance from the background MC samples, 
and comparing this to MC samples with signal of a certain 
flux. The number of skies with signal, that have more 
neutrinos in the given bin than the chosen significance 
from background only, is counted and gives a direct 
measure of the discovery potential for that particular flux. 

Figure |3] shows the discovery potential for 3ct (dashed 
line) and 5a (solid line) discovery, for an optimized bin 
of 4.9 degree bins, and a tolerance value accounting for 
a magnetic deflection of 3 degrees. Around 125 signal 
events on the whole sky are needed for a 5 a discovery 
in 50% of trials, and around 75 events are needed for a 
3(7 hint. Those values correspond to source flux values of 
about 2.16x10-^ GeV cm'^ s^^ and 1.29x10"^ GeV 
cm^^ s^^ respectively. 

VII. Results from the unblinded data 

To analyze the level of correlation between the distribution 
69 UHECRs reported by the Pierre Auger Observatory, 
and ANTARES neutrino candidates, we unblinded 2190 
neutrino events. The significance of observed correlation 
is determined with the help of randomized background 
samples, using the bin size of 4.9 degrees. 

The most probable count for the optimized bin size of 4.9 
degrees, or the mean background expectation from the 
randomized samples is 310.49 events (in all 69 bins, i.e. 
on the whole sky), with the standard deviation of 15.22 
events. 

After unblinding 2190 ANTARES neutrino candidate 
events, a count of 290 events within 69 bins is ob- 
tained. On Figure [4] neutrino candidates are represented 
with crosses and neutrino candidates correlating with 
observed UHECRs are highlighted as large triangles. The 
observed number of correlated neutrino events is beneath 
expected (negative correlation), with the significance of 
about 1.35(7. This result is compatible with a background 
fluctuation. The corresponding upper flux limit, assuming 
the equal flux from all UHECR sources, is 4.96x10^^ 
GeV cm-2 ^-i^ 
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Fig. 1: The probability density functions (PDFs) of the 
number of neutrino events in 1-10° bins centered on 69 
UHECRs observed by the Pierre Auger Observatory. The 
counts were obtained from the background only Monte 
Carlo simulations, each with 2190 neutrino events blinded 
by scrambling the observed events in right ascension. 
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Fig. 2: The mean upper flux limits as a function of the 
search bin size. The minimum flux value for magnetic 
field deflection tolerance value of 3 degrees is obtained 
for 4.9° search bin. 
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Fig. 3: The discovery potential of 3(7 (dashed line) and 
5(7 (solid line) as a function of number of signal neutrino 
events on the whole sky derived from the signal source 
flux. 
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Fig. 4: Crosses represent ANTARES neutrino events outside of 4.9 degree bins centered on UHECRs observed by the 
PAO, and triangles represent ANTARES neutrino events correlating with observed UHECRs. 
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Abstract. Cataclysmic cosmic events can be plausible 
sources of both gravitational waves (GW) and high- 
energy neutrinos (HEN). Both GW and HEN are 
alternative cosmic messengers that may traverse very 
dense media and travel unaffected over cosmological 
distances, carrying information from the innermost 
regions of the astrophysical engines. Such messengers 
could also reveal new, hidden sources that have not 
been observed by conventional photon astronomy. 

A neutrino telescope such as ANTARES can determine 
accurately the time and direction of high energy 
neutrino events, and a network of gravitational wave 
detectors such as LIGO and VIRGO can also provide 
timing/directional information for gravitational wave 
bursts. Combining these informations obtained from 
totally independent detectors can provide original ways 
of constraining the processes at play in the sources, 
and also help confirming the astrophysical origin of a 
HEN/GW signal in case of concomitant observation. 

This contribution describes the motivations and reach 
of a joint GW+HEN search using concomitant data 
taken with the ANTARES, VIRGO and LIGO detec- 
tors. The specific strategies developed for the selection 
of HEN candidates in ANTARES and the combination 
of HEN and GW data are presented with a focus on the 
currently ongoing analysis of the concomitant datasets 
taken in 2007 during the VIRGO VSRl and LIGO 
S5 science runs, while ANTARES was operating in a 
5-line configuration. 

I. Introduction and motivations 

The recent years have witnessed the development and 
operation of a new generation of detectors offering un- 
precedented opportunities to observe the universe through 
all kind of cosmic radiations. In particular, both high- 
energy (3>GeV) neutrinos (HENs) and gravitational waves 
(GWs), which have not yet been directly observed from 
astrophysical sources, are considered as promising tools 



for the development of a multi-messenger astronomy (see 
e.g. fT] and |2| for recent reviews). Contrary to high- 
energy photons (which are absorbed through interactions 
in the source and by the photon backgrounds) and charged 
cosmic rays (which are deflected by ambient magnetic 
fields), both HENs and GWs can escape from the core 
of the sources and travel over large distances through 
magnetic fields and matter without being altered. They are 
therefore expected to provide important information about 
the processes taking place in the core of the production 
sites and they could even reveal the existence of sources 
opaque to hadrons and photons, that would have remained 
undetected so far 

It is expected that many astrophysical phenomena pro- 
duce both GWs, originating from the cataclysmic event 
responsible for the onset of the source, and HENs, as a 
byproduct of the interactions of accelerated protons (and 
heavier nuclei) with ambient matter and radiation in the 
source. The detection of coincident signals in both these 
channels would then be a landmark event and sign the 
first observational evidence that GW and HEN originate 
from a common source. The most plausible astrophysical 
emitters of GWh-HEN are presented in Section |ll] 

The feasibility of common GW and HEN astronomies is 
also related to the concomitant operation of both types 
of detectors, which is summarized in the time chart of 
Fig. [T] Section III briefly describes the detection princi- 
ples and the performances achieved by the ANTARES 121 
neutrino telescope as well as by the GW interferometers 
VIRGO [4] and LIGO [6], that are currently part of 
this joint search prograrrQ As both types of detectors 
have completely independant sources of backgrounds, 
the correlation between HEN and GW significances can 
also be exploited to enhance the sensitivity of the joint 



L Similar studies involving the IceCube neutrino telescope l5l, 
currently under deployment at the South Pole and looking for 
neutrinos from astrophysical sources in the Northern Hemi- 
sphere, are also under way|j24|. 
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Fig. 1: Time chart of tiie data-taking periods for the ANTARES, 
VIRGO and LIGO experiments, indicating the respective up- 
grades of the detectors (as described in the text). The deployment 
of the KM3NeT neutrino telescope is expected to last three 
to four years, during which the detector will be taking data 
with an increasing number of PMTs before reaching its final 
configuration. 1231 . 



channel, even in absence of detection. The combined false 
alarm rate is indeed severely reduced by the require- 
ment of space-time consistency between both channels. 



In Section IV we describe in more detail the strategies 
being developed for a joint GWh-HEN searches between 
ANTARES and the network of GW interferometers using 
the currently available dataset^ 

II. ASTROPHYSICAL EMITTERS OF GRAVITATIONAL 
WAVES AND HIGH-ENERGY NEUTRINOS 

Potential sources of GWs and HENs are likely to be 
very energetic and to exhibit bursting activity. Plausible 
GWh-HEN emission mechanisms include two classes of 
galactic sources which could be accessible to the present 
generation of GW interferometers and HEN telescopes. 

Microquasars are believed to be X-ray binaries involving 
a compact object that accretes matter from a companion 
star and re-emits it in relativistic jets associated with 
intense radio (and IR) flares. Such objects could emit 
GWs during both accretion and ejection phases |7|; and 
the latter phase could be correlated with a HEN emission 
provided the jet has a hadronic component |8|. Soft 
Gamma Repeaters (SGRs) are highly magnetized X- 
ray pulsars with a soft 7-ray bursting activity. In the 
popular magnetar model [91, the outbursts are caused by 
star-quakes associated to large-scale rearrangements of the 
magnetic field. The deformation of the star during the 
outburst could lead to the emission of GWs within range 
of current interferometric GW detectors |10|. Sudden 
changes in the large magnetic fields would also accelerate 
protons or nuclei interacting with thermal radiation and 
generating a flux of HENs [1 1 1 . 

The most promising class of known extragalactic bursting 
sources are Gamma-Ray Bursts (GRBs). In the prompt 
and afterglow phases, HENs (10^ — 10^° GeV) are ex- 
pected to be produced by accelerated protons in relativistic 
shocks and several models predict detectable fluxes in 
km'^-scale detectors |12|, although no evidence for GRB 
neutrinos has been observed yet by IceCube 40 |13|. 
While gamma-ray and HEN emission from GRBs are 
related to the mechanisms driving the relativistic outflow. 



GW emissions is closely connected to the central engine 
and hence to the GRB progenitor. 

Short-hard GRBs are thought to originate from coalesc- 
ing binaries involving black holes and/or neutron stars; 
such mergers could emit GWs detectable from relatively 
large distances, with significant associated HEN fluxes 
[14|. Long-soft GRBs are most probably induced by 
"collapsars", i.e. collapses of a massive star into a black 
hole, with the formation of an accretion disk and a jet that 
emerges from the stellar envelope 1151 . This model also 
involves the emission of a strong burst of GWs during the 
gravitational collapse of the (rapidly rotating) progenitor 
star and in the pre-GRB phase; however this population is 
distributed over cosmological distances so that the associ- 
ated HEN signal is expected to be faint [161. The subclass 
of low-luminosity GRBs, with 7-ray luminosities a few 
orders of magnitude smaller, are believed to originate from 
a particularly energetic type Ibc core-collapse supernovae, 
although the mechanism responsible for their onset is still 
debated I.17J . They could produce stronger GW signals 
together with significant high- and low-energy neutrino 
emission; moreover they are more frequent than typical 
long GRBs and often discovered at shorter distances li 1 81 . 
Finally, the choked GRBs are thought to be associated 
with supernovae driven by mildly relativistic, baryon-rich 
and optically thick jets, so that no 7-rays escape [19J. Such 
"hidden sources" could be among the most promising 
emitters of GWs and HENs, as current estimations predict 
a relatively high occuiTence rate in the volume probed by 
current GW and HEN detectors 



III. The DETECTORS 

The ANTARES detector ||2l is the first undersea neu- 
trino telescope; its deployment at a depth of 2475m in 
the Mediterranean Sea near Toulon was completed in 
May 2008. It consists in a three-dimensional array of 
884 photo multiplier tubes (PMTs) distributed on 12 lines 
anchored to the sea bed and connected to the shore 
through an electro-optical cable. Before reaching this final 
(12L) setup, ANTARES has been operating in various 
configurations with increasing number of lines, from one 
to five (5L) and ten (lOL). 

ANTARES detects the Cherenkov radiation emitted by 
charged leptons (mainly muons, but also electrons and 
taus) induced by cosmic neutrino interactions with matter 
inside or near the instrumented volume. The knowledge 
of the timing and amplitude of the light pulses recorded 
by the PMTs allows to reconstruct the trajectory of the 
muon and to infer the arrival direction of the incident 
neutiino. The cuiTent reconstruction algorithms achieve an 
angular resolution (defined as the median angle between 
the neutrino and the reconstructed muon) of about 0.4° for 
neutrinos above 10 TeV|'9'l. The design of ANTARES is 
optimized for the detection of up-going muons produced 

2. More information can be obtained on the website of the first 
Workshop on Gravitational Waves and High-Energy Neutrinos: 
http://www.gwhen-2009.org . 
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Fig. 2: Instantaneous common sky coverage for VIRGO + LIGO 
+ ANTARES in geocentric coordinates. This map shows the 
combined antenna pattern for the gravitational wave detector 
network (above half-maximum), assuming that ANTARES has 
100% visibility in its antipodal hemisphere and 0% elsewhere. 
The colour scale is from 0% (left, blue) to 100% (right, red). 



by neutrinos which have traversed the Earth and interacted 
near the detector; its field of view is ^ 2ti sr for neutrino 
energies 100 GeV < < 100 TeV. 
The data acquisition system of ANTARES is based on 
the "all-data-to-shore" concept, which allows to operate 
different physics triggers to the same data in parallel. 
Satellites looking for GRBs can also trigger the detector 
in real time via the GCN (Gamma-Ray Burst Coordinate 
Network) alert system. ANTARES has also implemented 
the possibility to trigger an optical telescope network 
on the basis of "golden" neutrino events selected by 
a fast, online reconstruction procedure I22J. All these 
characteristics make this detector especially suited for the 
search of astrophysical point sources, and transients in 
particular ANTARES is intended as the first step towards 
a km"^ -sized neutrino telescope in the Mediterranean Sea, 
currently under R&D in the framework of the KM3NeT 
Consortium [23 J and whose construction could start by 
2014. 

The GW detectors VIRGO PI, with one site in Italy, 
and LIGO (see e.g. |6|), with two sites in the United 
States, are Michelson-type laser interferometers. They 
consist of two light storage arms enclosed in vacuum 
tubes oriented at 90° from each other. Suspended, highly 
reflective mirrors play the role of test masses. Current 
detectors are sensitive to relative displacements (hence 
GW amplitude) of the order of lO^^o to lO^^s Hz^^/^. 
Their detection horizon is about 15 Mpc for standard 
binary sources. 

The first concomitant data-taking phase with the whole 
VIRGO/LIGO network (VSR1/S5) was carried on in 
2007, while ANTARES was operating in 5L configuration 
(see Fig. [T]i. A second data-taking phase was conducted 
between mid-2009 and end 2010 with upgraded detectors 
(S6/VSR2 and VSR3), in coincidence with the operation 
of ANTARES 12L. Another major upgrade for both 
classes of detectors is scheduled for the upcoming decade: 
flie Advanced VIRGO/Advanced LIGO and KM3NeT 
projects should gain a factor of 10 in sensitivities respect 



to the presently operating instruments. As can be seen 
from Figure |2j the VIRGO/LIGO network monitors a 
good fraction of the sky in common with ANTARES: 
the instantaneous overlap of visibility maps is about 4 sr 
(- 30% of the sky). 

IV. Joint GWh-HEN searches using ANTARES 

DATA 

GW interferometers and HEN telescopes share the chal- 
lenge to look for faint and rare signals on top of abun- 
dant noise or background events. The GWh-HEN search 
methodology involves the combination of GW/HEN can- 
didate event lists, with the advantage of significantly 
lowering the rate of accidental coincidences. Preliminary 
studies on the feasibility of such searches 11241 . ||25]| 
indicate that, even if the constituent observatories provide 
several triggers a day, the false alarm rate for the combined 
detector network can be kept at a very low level (^1/(600 

yr)). 

An important ingredient of these searches is the definition 
of an appropriate coincidence time window between HEN 
and GW signals hypothetically arriving from the same 
astrophysical source. A case study that considered the 
duration of different emission processes in long GRBs, 
based on BATSE, Swift and Fermi observations, allowed 
to derive a conservative upper bound tew — tHEN € 
[—500s, +500s] on this time window ll26ll . For short 
GRBs, this time-delay could be as small as a few seconds. 

Different strategies for a joint search are currently being 
developed. The first one consists in an event-per-event 
search for a GW signal correlating in space and time with 
a given neutrino event considered as an external trigger. 
Such a search is rather straightforward to implement 
as it allows to make use of existing analysis pipelines 
developed e.g. for GRB searches. It is currently being 
applied to the concomitant set of data taken between 
January 27 and September 30, 2007 with ANTARES 
5L/VSR1/S5. 

Fig. |3] shows the schematic flowchart of such a HEN- 
triggered search. The ANTARES 5L data are filtered 
according to quality requirements similar to those select- 
ing the well-reconstructed events that are used for the 
standalone searches for HEN point sources. The list of 
candidate HENs includes their arrival time, direction on 
the sky, and an event-by-event estimation of the angular 
accuracy, which serves to define the angular search win- 
dow for the GW. For the purpose of this joint search, the 
angular accuracy is defined as the 90% quantile (and not 
the median) of the error distribution on the reconstructed 
neutrino direction, obtained from Monte Carlo studies. 
The on-source time window is taken to be [—500s, +500s] 
around the neutrino arrival time. 

The list of HEN triggers is then passed to the X- 
pipeline fZTl, an algorithm which performs coherent 
searches for unmodelled bursts of GWs on the combined 
stream of data coming from all active interferometers 
(ITFs). The background estimation and the optimization 
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Fig. 3: Schematic flow diagram of a HEN-triggered search for 
GWs, as described in the text. 

of the selection strategy are performed using time-shifted 
data from the off-source region in order to avoid con- 
tamination by a potential GW signal. Once the search 
parameters are tuned, the box is open and the analysis is 
applied to the on-source data set. If a coincident event 
is found, its significance is obtained by comparing to 
the distribution of accidental events obtained with Monte- 
Carlo simulations using time-shifted data streams from the 
off-source region. 

Alternatively, comprehensive searches for time- 
coincidences between independent lists of neutrino 
and GW events can also be performed, followed by a 
test of spatial correlation using the combined GW/HEN 
likelihood skymap. This second, more symmetrical, option 
requires the existence of two independent analysis chains 
scanning the whole phase space in search for interesting 
events. A major issue for both types of analyses lies in 
the combined optimisation of the selection criteria for 
the different detection techniques. 

V. Conclusions 

The joint GWh-HEN analysis program described here 
could significantly expand the scientific reach of both GW 
interferometers and HEN telescopes. Beyond the bene- 
fit of a potential high-confidence discovery, coincident 
GW/HEN (non-)observation shall play a crucial role in 
our understanding of the most energetic sources of cosmic 
radiation and in constraining existing models. A first 
search for GW signals in coincidence with HEN neutrinos 
detected by ANTARES 5L in 2007 is being performed, 
and optimisation studies are ongoing to exploit at best the 
concomitant datasets taken in 2009-2010 with upgraded 
detectors. 
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Abstract. Using the ANTARES neutrino telescope, the 
largest neutrino telescope in the Northern hemisphere, 
with its first configuration with 5 lines of photode- 
tectors to the actual nominal one corresponding to 
a total of 12 Unes, we have studied our ability to 
search indirectly for an evidence of Dark Matter 
annihilations in heavy astrophysical objects as the 
Sun and the Galactic centre. First results have been 
obtained using the data recorded by ANTARES in 
2007 and 2008, and compared with neutrino fluxes 
predicted within a minimal supersymmetric extension 
of the Standard Model with supersymmetry-breaking 
scalar and gaugino masses constrained to be universal 
at the GUT scale, the CMSSM, as well as a minimal 
Universal Extra-Dimensions scenario with one extra 
compact dimension where all the Standard Model 
fields propagate into the bulk, the mUED. The current 
limits over the neutrino/muon fluxes coming from 
Dark Matter annihilations, and the spin-dependent 
cross-section with protons, as well as the expected 
sensitivities predicted after several years of data taking 
with ANTARES will be presented for each source. 

I. Indirect dark matter search and the 

ANTARES NEUTRINO TELESCOPE 

As the very old problem about the unseen planets, measur- 
ing a deviation from the known laws of gravitation and the 
theory of general relativity in large astrophysical systems, 
we assume the existence of a huge amount of unseen 
matter, the dark matter 1 1 ], which represents 83 % of the 
matter in the universe. Massives, and weakly interacting 
with the matter, the WIMPs (Weakly Interacting Massive 
Particles), defining the dark matter nature whatever a 
particular framework in Particle Physics, can be elastically 
scattered in a medium to be decelerated, and gravitation- 
ally attracted in heavy objects as the galactic center, or the 
stars like our Sun which is the point source of the present 
analysis. It's through the self-annihilation of WIMPs into 
the Sun's core, that muon neutrinos (independently of 



electron, and tau neutrinos) can be produced in an energy 
range reachable to a neutrino telescope as ANTARES. 

Since the 29th of May 2008, the neutrino telescope 
ANTARES has been completed at 40 km offshore from 
Toulon at 2475 m depth, making it the largest neutrino 
telescope in the northern hemisphere, and the first to 
operate in the deep sea |2|. It is made up of 12 mooring 
detection lines holding ten-inch photomultiplier tubes, 
distributed into 25 storeys on each line, designed for 
the measurement of neutrino induced charged particles 
based on the detection of Cerenkov light emitted in water. 
With a low energy threshold of about 20 GeV for well 
reconstructed muons, and an angular resolution of about 
1° at low energy {E^ < 1 TeV), ANTARES has a 
great potential to detect neutrinos induced by the self- 
annihilation of WIMPs into the Sun's core. 



II. Dark Matter signal, and efficiency of 
ANTARES 

Model-independently, the differential amount of neutrinos 
dN^/dz as a function of the ratio z — E^/Mwimp, 
where Mwimp is the mass of the WIMP, which can be 
created from the Sun's core, coming to the surface of the 
Earth, is computed using the software package WIMPSIM 
Q. With a high statistics of production (12 x 10'^ of 
generated events inside the Sun) for all possible self- 
annihilation channels which can be reached as qq, II, (f)(f>*, 
and vD. Fig. [T] gives a glimpse of the main channels in 
the two common frameworks studied here: CMSSM, and 
mUED. No branching ratios have been taken into account 
for the three first channels (WW, bb, and rf), and all 
the contributions have been merged in the mUED case 
taking into account the following branching ratios ID: 
0.23 for rf, 0.077 for tt, 0.005 for bb, and 0.014 for 
i^D. All of these are shown for Mwimp = 350 GeV as 
an example. All of the physical contributions have been 
taken into account as the three-flavors oscillations, or the 
T lepton regeneration through the Sun medium. 
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Fig. 1: Self-annihilation channels WW, bb, tt, and the 
mUED particular case (from the top left to the bottom 
right). 



Then, an averaged effective area A^f /{Mwiaip), with 
the Sun as a point source, can be defined as a function 
of the WIMP masses, and a couple of cuts (Q,^'), where 
Q is a track fit quality cut based on a fast reconstruction 
algorithm (see |5| for more details), and 5* is the half- 
cone angle around the source. The following expression 
defines A^f f{MwjMp) as: 



Aeff{MwiMp) — 



Mw 



^ np.r).. 



X A, 



where Nrec. {Q, is the number of reconstructed neu- 
trino events considering the cuts couple (Q,^), iVge„. 
is the number of isotropically generated events, Ag^n. 
is the surface used in the event generation, Psun is 
the visibility of the Sun in the data taking period, and 
P{dN^ I dz) is the dark matter muon neutrinos signal 
energy spectrum normalized to one. Fig. |2] shows (top) 
the simple multiplication, over the neutrino energy range 
e [10 GtN;MwiMp] (Mwimp = 350 GeV here), of 
the effective area by the normalized differential amount of 
neutrino events at the surface of the Earth. As an example, 
a Q < 1.4 has been chosen with a half-cone angle 5" < 3°. 
Five different cases are presented here. First the effective 
area "alone" ^e^y in the direction of the Sun without 
any differential flux taken into account, then the effective 
area multiplied by the normalized differential amount of 
neutrinos for the W^W~ channel, for the bb channel, 
for the TT, and for the mUED merged channels. Even 
if the contribution in number of neutrinos is very large 
at very low energy regime for all of these channels, the 
final number of events is going to be smoothed by the 
weakness of the effective area at very low energy. Then, 
the major channels, in sense of neutrinos production, are 
going to be the WW, and the tt, especially for the 
CMSSM framework. For the mUED, as this model is 
dominated by the production of leptons/anti-leptons due 





Fig. 2: Simple multiplication of the effective area ^e// 
(m^) by the normalized differential amount of neutrino 
events at the surface of the Earth, over the neutrino 
energy range E,y £ [10 GeiV;M\YiMp], on the top. And, 
normalized effective area A^f f (Mwimp) {m?) for a 
Q < 1.4, and a half-cone angle 5* < 3° around the Sun, 
as a function of different WIMP masses, on the bottom. 



to a branching ratio closed to 0.23 H, the number of 
neutrinos produced is going to be very closed to the one 
from the rf channel, isolated in this plot for the CMSSM. 
Finally, a bump localized on the WIMP mass, caused 
by the direct production of neutrinos, is visible, direct 
production which is particular to the mUED framework. 

The integration of the previous distributions A^ff x 
P{dN^/dz) concludes on the effective areas dedicated 
to the WIMP masses A^ff (Mwimp), and this for 
each channel considered here in fig. |2] (bottom). These 
Aef f (Mwimp) have been computed for each (Q,^), and 
for a set of thirteen different WIMP chosen masses (10, 
25, 50, 80.3, 100, 150, 200, 250, 350, 500, 750, and 
1000 GeV). As an example, the same set of cuts and 
style of lines has been chosen than the previous plot, 
for the WW, bb, tt, and mUED merged channels. The 
bin centered around a mass Mwimp = 80.3 GeV has 
been isolated of the others because of its value of mass, 
closed to the W boson mass (Mw = 80.399± 0.023 GeV, 
PDG [6 1), corresponding to the beginning of the WW 
channel contribution (explaining the lacks in the mass 
bins centered on the 10, 25, and 50 GeV), particular to 
a MSSM framework, and helicity vanished in the general 
UED framework. Main and equivalent contributions from 
the WW and the rf channels (usually called "Hard" 
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Fig. 3: Sentivity (^^^^" in neutrino flux as a function of the 
WIMP masses, in the range imp €[10 GeV;l TeV] 
and from the data 2007-2008, for the channels W+W', 
bb, TT, and mUED. 



channels) still appeared, and a lower contribution from the 
bb channel (usually called "Soft" channel). Concerning the 
mUED case, its integrated effective areas are very closed 
to the TT ones, main channel in this framework. Finally, a 
lack in statistics from the effective area and the dark matter 
simulations didn't allow to compute all the channels for 
the WIMP masses at 10 (no channels), and 25 GeV (bb is 
missing). 

III. Results and discussions 

Following the Model Rejection Factor (MRF) tech- 
ruque Q to optimize the (Q,VE') cuts couple in the Sun 
direction, the sensitivity in neutrinos cf)^^'^" of ANTARES 
must be maximized for each value of WIMP mass, 
and each considered channel in the dark matter self- 
annihilation. For this, tf)^^'^° can be expressed as: 



790% 



90% 



^eff{MwiMp) X Teff ' 



(19) 



where is the average upper limit considering a Pois- 
son statistics in the Feldman-Cousins approach [SJ, which 
has been computed from a scrambled data set extracted 
from the 2007-2008 data taking period of ANTARES, 
where the configuration of the detector has evolved from a 
5 detection lines structure to the final one at 12 lines. Then, 
Tef / ~ 292.9 days is the effective time of 2007-2008 data 
taking period, and A^f f{MwiMp) is the effective area 
defined inllll 

As a direct consequence of the (Q,'^) couple optimization, 
a representation of the sensitivity 0^"^" can be produced 
in fig. [3] The most stringent sensitivities come from the 
"Hard" channels as W^W~ and tt, for the Supersym- 
metry, explained by the amount of neutrinos produced 
through these channels (see injlljl, and by extension due to 
the effective area Aeif{M-wiMp) (fig- 12] on the bottom) 
developped from them. For the particular case of the 
mUED framework, its sensitivity is very similar to the 
"Hard" channels previously mentioned because of the tt 
channel branching ratio supremacy (see in 
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Fig. 4: Sentivity (\?^'^'' in muon flux as a function of 
the WIMP masses, in the range MwiMp/inllO GeV;l 
TeV] and from the data 2007-2008, for the CMSSM 
framework, on the left side, with Baksan 1978 — 1995 (in 
grey). Macro 1989 — 1998 (in black), SuperKamiokande 
1996 - 2001 (in yellow), IceCube-22 for the bb chan- 
nel in 104 days (in green dashed line), IceCube-22 for 
the W^W~ channel in 104 days (in green solid line), 
ANTARES 2007-2008 for the bb channel (in red dashed 
line), ANTARES 2007 - 2008 for the VF+H^" channel 
(in red solid line), and ANTARES 2007-2008 for the tt 
channel (in red dot-dashed line). And the right side, the 
equivalent for the mUED framework. 



After this, to compare properly the sensitivity in neutrino 
flux from ANTARES 2007-2008 data to the other experi- 
ments, the sensitivity in muon flux must be produced, us- 
ing the cross section of neutrinos with the Earth medium, 
the muon range, the nucleon density in the vicinity of the 
detector, and the neutrino transmission probability through 
the Earth, all derived from the atmospheric neutrinos (anti- 
neutrinos) Monte-Carlo. Fig. |4] (top) shows it for the 
CMSSM framework. 

In the same way, the equivalent can be reached for 
the mUED framework in fig. |4] (bottom). Unfortunately, 
any sensitivities from other experiments are not officially 
available in time for this proceeding. Concerning the 
present one for ANTARES, the shape is conformed to 
what is expected. This sensitivity in muon flux is very 
similar to the W^W^ one, or tt one, from the CMSSM 
framework due to the fact that the mUED neutrino flux 
construction is dominated by the tt channel through 
its branching ratio (see in [ll]). More the effective area 
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Aef f{MwiMp) is large for high mass Mwimp values 
(see fig. 2|i, more the resulting is weak. 

All of these sensitivities, for both frameworks studied in 
this paper, can be compared to the theoretical parameter 
spaces allowed by the experimental constraints, and de- 
rived from the SuperBayes simulation for the CMSSM 
framework |9|, and for mUED IITOl . In this way, the spin- 
dependent cross-section with protons d-^'gj^ is developped 
from the sensitivity using the method described 

in ifm . as it appears in fig. |5] respectively for the 
CMSSM (top), and the mUED (bottom). For CMSSM 
and mUED, more the theoretical parameter space has 
a clear color more this one is allowed by the actual 
experimental contraints, with a comparison to the present 
direct detection experiment as KIMS |131, CDMS lfT2l . 
COUPP 1 14], or Picasso 1 15 1, and the usual indirect exper- 
iment like IceCube-22 [16J, and SuperKamiokande ifTTl . 
These curves show a bb channel sensitivity very closed to 
the SuperKamiokande one, and better than the IceCube- 
22 one for the same channel with more than one order of 
magnitude for the low energy regime (m^ < 180 GeV), 
and a W~^W~ channel sensitivity better than the IceCube- 
22 one until ^ 440 GeV, to be very closed to this 
last one for higher mass values. Concerning the rf 
channel, no equivalent has been found for IceCube or 
others, which places ANTARES as the actual best one for 
this channel in the CMSSM framework. This fact can be 
explained by the energy threshold of ANTARES around 
20 GeV, lower than for IceCube ~ 50 GeV, which has an 
effect over all the WIMP mass range due to the integration 
over the neutrino energy range, E^, €[10 GeV;A/vi/7Afp]- 
For the mUED case which appears in fig. |5] (bottom), we 
took the same line style code than the previous one con- 
cerning the CMSSM (top) for the direct experiment. The 
contributions of the indirect search experiments is always 
a great complement to the direct detection experiments, 
with a best sensitivity with ANTARES on all the Mlkp 
mass range. 
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Abstract. Although the first evidence for neutrino 
oscillations came from measurements on atmospheric 
neutrinos in underground experiments, neutrino oscil- 
lations have yet to be demonstrated in high energy 
neutrino telescopes, whose energy threshold is signifi- 
cantly higher. Recent studies have shown that a clean 
sample of atmospheric neutrinos with energies as low 
as 20GeV can be isolated in the ANTARES neutrino 
telescope. Such a threshold is low enough to allow the 
observation of neutrino oscillation features. A robust 
analysis method is presented which allows the extrac- 
tion of atmospheric neutrino oscillation parameters. 



possible to extract a clean sample of atmospheric neutrinos 
with energies around the first minimum in the survival 
probability of muon neutrinos propagating through the 
Earth. 

Even though dedicated experiments are more sensitive for 
the measurements of neutrino oscillations, such an anal- 
ysis is an important check for ANTARES. The analysis 
being very sensitive not only to the telescope efficiency 
but also to the quality of the angular and energy recon- 
struction, it becomes a benchmark for the understanding 
of the detector. 



I. Introduction 

It is now well established that neutrinos can switch 
from a flavour to another, neutrino flavour eigenstates 
for the weak interaction being different from neutrino 
mass eigenstates. This phenomenon, hypothesized in 1957 
by Pontecorvo and described more precisely in the early 
sixties yj, is known as neutrino oscillations, and has been 
measured by a number of experiments 

Although part of the scientific program of the ANTARES 
neutrino telescope, oscillations studies have been delayed 
because of the difficulty to reliably reconstruct muons 
at a sufficiently low energy. Indeed, the granularity of 
ANTARES is rather coarse compared with the range 
of GeV muons in water : its 885 photomultiplier tubes 
are distributed by triplets placed 14.5 m apart along the 
detector lines, themselves spaced by about 60 m |2|. As 
a consequence, the energy range where the effect of 
the atmospheric neutrino oscillations should be visible 
is situated at the very edge of the detector sensitivity, 
which makes the neutrino oscillations measurement very 
challenging. 

The topic experienced a regain of interest since the 
completion of the detector in 2008. As progress was 
steadily made in the detailed understanding of the detector 
and of its environment, it became clear that it should be 



II. Purpose 

Neutrinos oscillations are commonly described in terms 
of E/L dependence, where E is the neutrino energy and 
L its oscillation path length For a neutrino telescope 
such as ANTARES, detecting neutrinos crossing the Earth, 
L can be translated as 2R sin 9, R being the Earth radius 
and 9 the anti-elevation (that is, the angle between the 
neutrino direction and the horizontal axis : 9 = 7r/2 
for a vertical upgoing neutrino). Within the two-flavour 
approximation, the i/^ survival probability can then be 
written 

P{i^t, ^ i^p) 1 - sin^ 26*23 sin^ ^2.54i?ATO^2^^^ > 

6*23 and Amgj being respectively the mixing angle and the 
squared mass difference of the involved mass eigenstates 
(with R in km, E in GeV and Am|2 in eV^). 
According to recent results from the MiNOS experi- 
ment fA\, the first minimum in the muon neutrino survival 
probability E/s\\\9 spectrum occurs at about 24 GeV 
(figure [T]). This interesting region is in principle accessible 
for ANTARES. Indeed, Monte-Carlo simulations have 
shown that it is possible to extract clean samples of 
atmospheric neutrinos with energies in this range, which 
should enable the observation of the first minimum. 
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Fig. 1: Survival probability of muon neutrinos as a 
function of E/sm6 in the two-flavour approximation, 
with sin^ 26*23 = 1 and l^rnl^ = 2.32 x lO"^ (the 
line width corresponds to MiNOS 1 a uncertainty on this 
parameter Q), neglecting matter effects. 



III. Method 

A simple and robust method to measure neutrino os- 
cillations in ANTARES consists in taking advantage of 
the different granularity of the detector along and across 
the detection lines. Indeed the range of muons below 
about 20 GeV makes them unlikely to give signal on 
different lines. They can be well reconstructed however 
using only one line, especially if they are close to the 
vertical. These low energy vertical tracks should be sup- 
pressed in the current oscillation scenarios compared to 
the more energetic and less vertical ones. Consequently, 
the observed ratio between the number of tracks giving 
signal on a single detector line (ID tracks) and the number 
of tracks giving signal on several detector lines (3D tracks) 
should differ from the expected ratio in the case of no 
oscillations. Many systematic effects such as the global 
flux normalisation or absolute detection efficiencies should 
cancel in this ratio, leading to a small residual systematic 
effect. The observed ratio can be expressed in terms of 
Monte-Carlo events as 

/id 



R = 
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Ajno osc 
^^3D 



23 • 



Sin^ 26I23 • /3D '' 



(20) 



where -^"disd'^ is the number of expected ID (3D) tracks 
without neutrino oscillations and 
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is summed over each event i. 

For each value of Am^s the mixing angle can be calcu- 
lated analytically from the measurement of R according 



to equation 20 The statistical and systematic errors from 
R propagate directly into an error on the mixing angle, 
which leads to a band in the (^23, Am^j) plane. 
Tracks are reconstructed using a fast reconstruction algo- 
rithm (51, which has the advantage of being more effi- 



cient at lower energies than other existing reconstruction 
strategies. A cut on the estimated quality of the fitted 
track and on the reconstructed angle is performed, in 
order to select well-defined events and to increase the 
upgoing neutrinos purity. The ID tracks are required to be 
reconstructed using hits from at least 8 detector storeys. 
Lowering this cut would considerably enhance the number 
of ID events and thus the sensitivity to oscillations, but 
would also increase dramatically the contamination from 
misreconstructed atmospheric muons. Such a hard cut 
ensures a very low contamination and leads to a minimal 
track length of about 100 m. The energy threshold is 
thus roughly 20 GeV, which is below the first oscillation 
minimum in the E / sin 6 spectrum. 

This method is robust as it does not depend on many 
assumptions, nor does it rely on an energy estimator. 
Furthermore the reconstruction algorithm is used only to 
assert the quality of the selected events. Finally, assuming 
they affect similarly ID and 3D events, the ratio should 
cancel most systematics, and in particular the large un- 
certainties on the atmospheric neutrino flux normalization 
and the uncertainties on the detector simulation. The strict 
hit selection criteria and the strong cut on the number of 
hit storeys minimize the sensitivity of the ratio to potential 
biases in the optical background simulation. Remaining 
systematic uncertainties are expected to be within a few 
percents, smaller than the current statistical errors. 

IV. Expected sensitivity 

Table |l] shows, according to MiNOS results H, the number 
of expected events for a 170 days Monte-Carlo sample 
(this corresponds to 2008 active time). The atmospheric 
neutrino flux is weighted to match Bartol parametriza- 
tion ill. The contamination of misreconstructed downgo- 
ing atmospheric muons is negligible. The number of ID 
events is suppressed by 16% in the case of oscillations 
while the number of 3D events is suppressed only by 
3.6 %. The effect to be observed is however small : the 
suppression of ID tracks concerns only 26 events for this 
sample. It is clear that the E/ smO spectrum cannot be 
reliably extracted under such conditions : a larger statistics 
is needed. 





no osc. 


osc. 


contamination 


ID tracks 


186 


160 


4.7 


3D tracks 


522 


504 


1.2 



TABLE I: Number of expected Vi_i/v^ charged current 
events, after 170 days of lifetime, with and without 
oscillations, for ID and 3D tracks, and number of mis- 
reconstructed atmospheric muons surviving the cuts. 

The sensitivity after 1000 days, which may be achieved 
after 4 or 5 years of real data taking depending on various 
external conditions, is presented in figure |2j extrapo- 
lating from the numbers given in table [l] The number 
of expected ID events is roughly one per day, which 
leads for 1000 days to a statistical error of about 3 % 
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ANTARES 1000 days sensitivity [PRELIMINARY] 




Fig. 2: ANTARES expected sensitivity to atmospheric 
neutrino oscillations after 1000 days. 

(~^30 events). Studies where PMT and water parameters 
have been varied within their tolerance have shown that 
a 3 % systematic error on R is realistic. Consequently, a 
total standard deviation of about 5 % on i? has been used 
to draw the measurement contours of figure |2] 

According to these results, although not competitive with 
dedicated experiments, ANTARES should be sensitive to 
neutrino oscillation parameters through disappearance of 
atmospheric muon neutrinos with the simple and robust 
analysis presented here. Current preliminary results using 
a restricted data sample are compatible with world oscil- 
lation data. Processing of the whole ANTARES data set 
is foreseen. 

V. Prospects 

Several directions are currently being investigated to im- 
prove ANTARES sensitivity to neutrino oscillations. In a 
near future, it should be possible to increase significantly 
the number of low energy reconstructed muons thanks 
to a dedicated reconstruction algorithm coupled to an 
optimized hit selection, which would reduce statistical 
uncertainties using the same amount of data. Furthermore 
the use of a different reconstruction algorithm would be an 
important cross-check of such an analysis. Additionally, 
the proportion of low energy events can be improved by 
selecting events contained in the detector 

Such improvements should enhance ANTARES sensitivity 
to neutrino oscillation parameters. They might also allow 
the extraction of the E / sin 9 spectrum, which would be a 
great opportunity to cross-check the understanding of the 
detector 

VI. Conclusions 

A robust method to extract atmospheric muon neutrino 
oscillation parameters from ANTARES data has been pre- 
sented. Although not competitive with dedicated experi- 
ments, ANTARES should be able to reach some sensitivity 
to these parameters, which would be a demonstration of 



the understanding of the ANTARES detector Preliminary 
analysis of a restricted data sample is compatible with 
existing constraints on neutrino oscillation parameters. 
The complete analysis of ANTARES data is ongoing. 
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Abstract. The ANTARES underwater neutrino tele- 
scope, located in the Mediterranean Sea at a depth 
of 2475 m, 40km off the Provencal coast, is composed 
by an array of 885 photomultipliers distributed on 12 
vertical lines. The detector is fully operational since 
May 2008. Besides the detection of neutrino-induced 
muons, the telescope is more generally sensitive to 
particles which emit Cherenkov light, and, thanks to 
its large volume, offers new opportunities to improve 
the sensitivity to exotic particles, such as magnetic 
monopoles. Magnetic monopoles are stable magneti- 
cally charged particles first introduced in quantum 
mechanics by Dirac in 1931, which showed that their 
existence would give an explanation to the quanti- 
zation of the electric charge. They would have been 
produced in the Early Universe, and would bring a 
first proof to the existence of grand unified models. 
The paper introduces magnetic monopoles and their 
signal caracteristics in a neutrino telescope, and then 
describes the analysis which uses a dedicated algo- 
rithm able to reconstruct the velocity of the incoming 
particles. A new upper limit on the upward going 
magnetic monopole flux, extracted from ANTARES 
data taken in 2008, is presented. It is, at present, 
the best worldwide constraint in the velocity range 
(3 = [0.625, 1], and provides, for the first time with a 
neutrino telescope, a constraint for monopole velocities 
below the Cherenkov threshold. 

I. Introduction 

Magnetic monopoles are hypothetical magnetically 
charged particles first introduced by Dirac in 1931. He 
demonstrated that the existence of magnetic monopoles 
naturally leads to the quantization of the electric 
charge [IJ. In 1974, 't Hooft and Polyakov discovered 
independently that, in certain spontaneously broken gauge 
theories, magnetic monopoles are not only a possibility, 
but a requirement Q, fJI. However, despite intensive 
search efforts, no magnetically charged particles were 
detected up to now, and stringent experimental flux lim- 
its beyond the original Parker bound of flux Fp ^ 
10^^^ cm^^s^^sr^^ [4| were provided. The development 



of huge detectors for performing neutrino astronomy gives 
rise to new expectations for the search for magnetic 
monopoles, and ANTARES offers new interesting oppor- 
tunities in this field. 

In this paper a search for upward going relativistic mag- 
netic monopoles for one year of data taking in 2008 with 
the ANTARES detector |5| is presented. Section |Il] will 
first introduce magnetic monopoles and their signature 
in the ANTARES neutrino telescope will be briefly pre- 
sented. Then the simulation and reconstruction algorithm 



will be detained in section III and finally the analysis 



method will be explained in section IV and the final 
results will be given in section [V] 

II. Magnetic monopoles and their signal in the 
ANTARES neutrino telescope 

A. Magnetic monopoles 

Most of the Grand Unified Theories (GUTs) predicts the 
creation of magnetic monopoles in the early Universe. 
Indeed, in 1974, 't Hooft and Polyakov showed indepen- 
dently that elements caracterising well a magnetic charge, 
as introduced by Dirac in 1931 |1|, occur as solitons in 
Gran Unified Theories (GUT) in which a larger gauge 
group breaks down into a semi-simple subgroup contain- 
ing the explicit U(l) group of electromagnetism O. 
These particles are topologically stable and carry a mag- 
netic charge defined as a multiple integer of the Dirac 
charge gp, = where e is the elementary electric 
charge, c the speed of light in vacuum and h Planck's 
constant. Depending on the specific model, the masses 
inferred for magnetic monopoles range over many orders 
of magnitude, from 10"* to lO^^ GeV [6|. 

Since fast monopoles have a large interaction with matter, 
they can loose large amounts of energy in the terres- 
trial environment. The total energy loss of a relativistic 
monopole with one Dirac charge is of the order of 
10^^ GeV after having crossed the full diameter of the 
Earth Q. Because magnetic monopoles are expected to 
be accelerated in galactic coherent magnetic field domains 
to energies of about 10^^ GeV |8|, calculations suggest 
that monopoles with mass below ^ 10^^ GeV would be 
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able to cross the Earth and reach the ANTARES detector 
as upward going signals. 

B. Signal in the ANTARES neutrino telescope 

The ANTARES detector is an underwater neutrino tele- 
scope immersed in the Western Mediterranean Sea at a 
depth of 2475 m 1 9 1 . In its full configuration the detector 
consists in 12 mooring Unes, each comprising 25 storeys 
separated by 14.5 m. A storey consists of 3 glass spheres 
housing a photomultiplier. Finally each line is connected 
to a Junction Box, which is connected to the shore station 
at La Seyne-sur-Mer by an electro-optical cable of 40 km 
length. 

As for electric charged particles, magnetic monopoles 
crossing ANTARES with a velocity higher than the speed 
of light in the sea water of /? ^ 0.74, are expected to 
emit Cherenkov light. The expected light is however much 
more intense, with about 8550 times more Cherenkov 
photons than muons of the same velocity [lOJ . 

In addition to the direct photon emission above the 
Cherenkov threshold, magnetic monopoles can produce 
indirect light below. Indeed, by ionizing the sea water 
along its path, a monopole with a velocity (5 > 0.51 
will knock out electrons ((5— rays), which will get enough 
energy for emitting Cherenkov light ifTTIl . Fig. [T] indicates 
the amount of light produced during the crossing of a 
monopole of charge g]j as a function of its velocity, 
compared to the emission of a through going muon. 
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Fig. 1: Number of Cherenkov photons in the 300 — 600 nm 
wavelength range per centimeter emitted by a monopole with 
g ~ go (solid line), by the 5— rays produced along its path 
(dashed line) and by a muon (dotted line), as a function of the 
particle velocity. 



III. Simulation and reconstruction 
A. Simulation 

Upward going magnetic monopoles with one unit of Dirac 
charge (g — go) have been simulated using a Monte 
Carlo program based on GEANT3 1X2^ for ten ranges 
of velocities in the region /3 = [0.550, 0.995] and with 
the incoming direction distributed isotropically over the 
lower hemisphere. The simulation of emitted photons is 
processed inside a cylindrical volume surrounding the 
instrumented volume. A radius of 480 m, four times 



larger than that used for the standard ANTARES muon 
simulation, is chosen in order to take into account the 
large amount of light emitted by a magnetic monopole. 

When searching for upward going magnetic monopoles, 
the main physical backgrounds are upward going at- 
mospheric neutrino-induced muons and downward going 
atmospheric muons. For simulating atmospheric muons, 
the CORSIKA program |fT3l in combination with the 
QGSJET llT4l code was used, employing the Horandel 
model for the cosmic ray spectrum, whereas atmospheric 
neutrinos were simulated according to the Bartol flux lfT6l . 
[17| and the RQPM model|18|. 

In order to match the real detector conditions, the sim- 
ulations have been performed using the dead channel 
information and the optical background rates extracted 
from the data. 

B. Reconstruction algorithms 

Before the reconstruction step, the different trigger algo- 
rithms are applied to data before events are stored [191 . 
These conditions lead to less than 15 % inefficiency for 
monopoles with a velocity greater than /3 — 0.58 that 
produced at least six hits in the detector. 

The standard track reconstruction assumes that particles 
travel at the speed of light. In order to improve the 
sensitivity for magnetic monopoles travelling with lower 
velocities, the velocity was introduced as a free parameter 
in the track fit. The modified algorithm performs two 
independent fits, a track fit and a bright point fit. The 
former reconstructs the track of a particle crossing the 
detector at a velocity /?, introduced as a free parameter, 
while the latter reconstructs the event as a point-like light 
source. Both fits minimize the same quantity defined 
as : 

(^7 ^i) 



X 



N 

E 



A,, 



(21) 



The first term on the right hand side is the sum for N 
hits of the square of the time residual, where is the 
expected arrival time of a hit, ti the measured time and 
CTi the time uncertainty for each hit. The second term Ai 
corresponds to the contribution of hit charge weighted 
by its minimal approach distance. It brings a penalty to 
hits with large amplitude combined with a large minimal 
approach distance. 

This algorithm yields a resolution 5/3 on the reconstructed 
velocity of magnetic monopoles of about 5/3 ~ 0.025 for 
velocities lower than the Cherenkov threshold, Pc = 0.74, 
improving to 5/3 ~ 0.003 for higher velocities. 

IV. Analysis 

The analysis has been performed for data taken from De- 
cember 2007 to December 2008. Quality requirements are 
applied to the data to ensure low levels of bioluminescent 
activity and a well calibrated detector After this selection 
the data is equivalent to a total of 136 days of live time : 
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44 days with 12 lines, 46 days with 10 lines and 47 days 
with 9 lines. 

For checking the agreement between simulated atmopheric 
events and real data, 15 % of the whole dataset, corre- 
sponding to 20 days of data taking, was taken and used at 
each step of the analysis. To avoid any bias, the sample 
of 15 % of data has at the end been removed, and so 
has not been considered in the final result. Despite the 
good agreement of shape between data and MonteCarlo 
atmospheric distributions observed, a normalisation factor 
of 1.8 is needed to be applied on simulated atmopheric 
muons, which is consistent with the expected uncertainties 
on the optical module angular acceptance for downgoing 
particles and on the input parameters of the model- 
dependent atmospheric muon flux. 




nhit 

Fig. 2: Distributions of nhit for events reconstructed with 
Prec ~ [0.775, 0.825]. Atmospheric down-going muons and up- 
going neutrinos (solid histograms) are compared to the 15 % data 
sample (points). The dashed line stands for magnetic monopoles 
generated in the same velocity bin. 



A. Preliminary cuts 

Only magnetic monopoles crossing the Earth are con- 
sidered in the analysis, and a natural precut is then to 
only keep particules reconstructed as upward going, i.e. 
with a reconstructed zenith angle less than 90°. This 
cut is motivated by the large amount of background of 
downward going events coming from atmospheric muons. 

Moreover a second preliminary cut consists in choosing 
only particles for which the track was reconstructed by 
using at least 2 lines, to improve the reconstruction quality. 

Finally, the last preliminary cut consists in removing 
events for which the quality factor of the fit is better with 
the bright point hypothesis than with the track one. Elec- 
tromagnetic and hadronic showers dominate such events, 
and a large fraction of poorly reconstructed atmospheric 
muons are then removed. 



B. Discriminating variables 

The event selection for optimising the Model Discovery 
Factor was performed by using two discriminating vari- 
ables. As shown fig. [T] the first discriminant between 
atmospheric events and magnetic monopoles is the number 
of hits nhit produced in the detector Fig.|2]represents dis- 
tributions of the number of hits used by the reconstruction 
algorithm for events reconstructed in the velocity range 
(3rec — [0.775,0.825]. A good agreement is observed 
between simulated atmospheric events and the 15 % data 
sample, whereas the distribution for magnetic monopoles 
simulated in the same velocity bin confirms the potential 
of discrimination for the nhit variable. 

However, for the lowest monopole velocities below the 
Cherenkov threshold, using only the number of hits is not 
sufficient for discrimating monopoles from atmospheric 
background events, and a new variable has thus been 
defined by applying two reconstruction algorithms. The 
first one is the initial muon fit algorithm in which the 
velocity /? is fixed at 1, whereas the second modified 
algorithm allows /3 as a free parameter in the procedure. 



The discriminating parameter A is then defined as 



A = log 



X?(/3 



I) 



XtiPfre 



(22) 



where XtiP — 1) Xt (/3/ree) are the track quality pa- 
rameters for a fixed and free velocity /3 respectively. This 
new variable is then expected to be positive for monopoles 
and negative for atmospheric muon and neutrino events. In 
fig.|3]are shown the A distributions of events reconstructed 
in the velocity range /3 — [0.775,0.825], for atmospheric 
simulated muons and neutrinos, and for the 15% data 
sample, as well as for magnetic monopoles simulated in 
the same range of velocity. 




Fig. 3: Distributions of A for events reconstructed with /?,.ec = 
[0.775,0.825]. Atmospheric down-going muons and up-going 
neutrinos (solid histograms) are compared to the 15 % data 
sample (points). The dashed line stands for magnetic monopoles 
generated in the same velocity bin. 



Finally, the optimisation of the model discovery factor for 
a 5(7 discovery at 90% probability was performed for 8 
bins of monopole velocity between /3 — [0.625, 0.995] by 
scanning the n/iit- A-plane. 

V. Results 

The [nhit, A] cuts for every velocity bins, defined pre- 
viously were then applied to the 85 % data sample. The 
number of atmospheric background events expected after 
116 days of data taking is indicated in Table 1 for each 
bin of reconstructed velocity, as well as the number of 
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observed events when the selection cuts are applied to 
the 85 % data sample. No significant deviations from the 
background-only hypothesis are observed by comparing 
the number of expected events and the number of observed 
events. All the observed events are therefore considered as 
background, with no claim for a monopole discovery. The 
Feldman-Cousins 90 % C.L. upper limit on the upward 
going magnetic monopole flux is reported in Table 1. 





Number of exp. 


Number of 


90% C.L. flux u. 1. 


range 


background events 


observed events 





[0.625, 
[0.675, 
[0.725, 
[0.775, 
[0.825, 
[0.875, 
[0.925, 
[0.975, 



0.6751 
0.7251 
0.7751 
0.8251 
0.8751 
0.9251 
0.9751 
1.0251 



2.2 X 

1.3 X 
4.5 X 

1.1 X 

8.2 X 
6.9 X 

2.3 X 
1.3 X 



10"^ 
10-1 
10"^ 
10-** 
10-^ 
10-^ 
10"^ 
10~^ 



7.5 X 10" 
8.9 X 10" 
4.0 X 10" 
2.4 X 10" 
1.8 X 10" 
1.7 X 10" 

1.6 X 10" 
1.3 X 10" 



TABLE I: Number of background events expected after 116 
days, compared to the number of observed events from the 85 
% of data. The Feldman-Cousins 90 % C.L. flux limit is also 
reported. 



In the flux upper limit calculation, systematic uncertainties 
were considered. They mainly depend on uncertainties 
on the detector efficiency. The major contributions to the 
uncertainties are the OM angular acceptance which con- 
tributes to 15 % and the uncertainty in the light absorption 
length in water which gives an uncertainty of 10 % \2\i . 
To estimate the effect of the loss of detection efficiency, 
18 % of hits per event were removed randomly in Monte 
Carlo monopole simulation. This leads to a deterioration 
of the upper limit of 3 % above the Cherenkov threshold 
and between 30 % and 7 % ffor velocities below. 

The flux upper limit for upward going magnetic 
monopoles is shown in Fig. |4] as a function of the 
monopole velocity, and is compared to previous limits. 



'E 10 " - 




VI. Conclusion 

The analysis performed on 1 16 days of data taking in 2008 
by the ANTARES detector, for the search for magnetic 
monopoles yields to new constraints on the upward going 
magnetic monopole flux. This new upper limit improves 
by a factor of about 3 the AMANDA'S one over the range 
0.75 < /? < 0.995 (7 = 10), whereas for the first time for 
a neutrino telescope, an analysis is sensitive and improves 
the flux upper limit below the Cherenkov threshold. 
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Fig. 4: The ANTARES flux upper limit (90 % C.L.) on up- 
going magnetic monopole with g = go for relativistic velocities 
0.625 < P < 0.995 is compared to the published upper limits 
set by MACRO |22| for an isotropic flux, and Baikal 123) 
and AMANDA [24 1 for up-going monopoles, as well as the 
theoretical Parker bound (4|. 
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Abstract. We discuss the search for down-going nucle- 
arites with the ANTARES neutrino telescope. After 
a very brief description of ANTARES, we explain 
the detection mechanism of slowly moving nuclearites. 
The search strategy, based on a blind analysis is 
then discussed. The search for nuclearites in data 
collected in 2007 and 2008 with different ANTARES 
configurations is ongoing. 

I. Introduction 

The ANTARES neutrino telescope was completely de- 
ployed in 2008 at a depth of 2475 meters, about 40 
km South of Toulon, France. It is composed of 12 
vertical lines, bearing a total of 885 optical modules 
(OMs) distributed in 25 storeys, at depths between 2050 
and 2400 meters. A smaller instrumentation line is also 
present. The OMs, consisting of a glass sphere housing 
a 10" Hamamatsu photomultiplier (PMT) are arranged in 
triplets per storey, with the axes looking 45° downwards. 
Each line is connected by an electro-optical cable to the 
main junction box, and then to the shore station 
A secondary junction box is also present, ensuring the 
connection of various sea science experiments. All data 
above a certain threshold are sent to shore, where they 
are filtered by various triggers before being stored f2l|. 
This is the so called "all data to shore" strategy. 

ANTARES is optimized for the detection of upward going 
relativistic muons, but is also sensitive to a variety of 
exotic particles, as fast magnetic monopoles and slow 
nuclearites lO, dU. In this paper we report on the status 
of nuclearite search with the ANTARES detector. 

II. Nuclearite properties 

If Strange Quark Matter (SQM) is the ground state of 
Quantum Chromodynamics (QCD), nuggets of SQM (so 
called "nuclearites") have to be present in the penetrating 
cosmic rays Q. They could have been produced in the 
early Universe, or in violent astrophysical processes as 
binary strange star collapses [6 J. The only phenomeno- 
logical bound to the nuclearite flux in the galaxy is 



derived from the dark matter density The nuclearite 
detection in very large volume neutrino telescopes is 
possible through the black body emission from their over- 
heated path in water |5|. Nuclearites are known also as 
"strangelets"; in most cases this denomination refers to 
nuclear SQM nuggets with mass similar to that of heavy 
nuclei, "nuclearites", as used in this paper, refer to much 
heavier objects (M > 10^'' GeV). The density of SQM 
is estimated to be a little larger than that of ordinary 
nuclear matter, p^v — 3.6 x 10^^ g cm^'^. According 
to [Tl, the chemical potential difference between s and 
u or d quarks should induce a small residual positive 
charge for the finite size SQM. Nuclearites moving with 
typical velocities of gravitationally trapped objects in the 
Galaxy (/? = w/c ~ 10^'^) should have the residual charge 
compensated by an electron cloud, or/and by electrons in 
weak equilibrium inside the SQM ijS). 

A. Nuclearite energy loss in matter 

For slowly moving nuclearites (/3 ~ 10^"^) the dominant 
energy loss is through elastic collisions with the atoms in 
the traversed medium: 



where a is the nuclearite cross section, and p the density 
of the medium. The cross section has the typical atomic 
value a ~ IT y. 10"^^ cm^ for nuclearite masses M < 
8.4 X 1014 GeV, and 7r(3Af/47rpjv)2/3 for larger masses. 
Assuming that outside the atmosphere the velocity of a 
nuclearite is (3 = 10^^, in order to cross the Earth its 
mass should be larger than about 10^^ GeV As we 
expect that the nuclearite flux in the cosmic rays decreases 
with increasing mass (as for heavy nuclei), we concentrate 
for this analysis only on downgoing nuclearites. 

B. Principle of nuclearite detection in ANTARES 

The energy release through elastic collisions (Eq. 1) 
results in the over-heating of the nuclearite track in 
matter. In the case of water, it is estimated that a fraction 
77 ~ 3 X 10^'' of the energy loss is dissipated as visible 
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Fig. 1: Light yield from a vertically down-going nuclearite 
at the level of the upper, middle and lower ANTARES 
storeys. 



black body radiation emitted by the expanding cylindrical 
shock wave |5|. The number of visible photons emitted 
per unit path length may be estimated as: 



dx 



dE/dx 



(24) 



where {E-^) ~ tt eV is the average energy of visible 
photons. 

As a simple exemplification Fig. [T] shows the light yield 
at the ANTARES depth assuming a vertically down- 
going nuclearite, with the initial velocity (above the Earth 
atmosphere) /3 ~ 10^"^, as a function of its mass. The 
curves correspond (from up to down) to the light yield 
at the level of the upper, middle and respectively lower 
ANTARES storeys. Nuclearites with masses larger than 
few 10^'^ GeV produce enough light to be detectable 
with the ANTARES telescope. A complete Monte Carlo 
simulation describing the light production along nuclearite 
paths, considering initial isotropic arrival direction from 
the upper hemisphere and distances from the detector axis 
expanding to 5 absorbtion lengths in water was developed 
and used in the analysis described bellow. 

III. Nuclearite search strategy in ANTARES 

The basic information in ANTARES is the "hit", the time 
and charge (amplitude) information of a photon detected 
by a PMT. If the charge is over a pre-defined threshold, 
the hit is called "LO" hit and buffered. A local coincidence 
("LI" hit) is defined as two LO hits in the same storey 
within 20 ns, or a single hit with large amplitude (3 
photoelectrons (pe) or more, collected by a single PMT). 
The "directional trigger" (DT) requires at least 5 LI 
casually connected hits anywhere in the detector, within a 
2.2 /IS window. A second trigger, implemented in a later 
phase, is the "cluster trigger" (CT). It is based on the 
"T3" cluster, a combination of two LI hits in adjacent or 
next-to-adjacent storeys, and requires two of such clusters 
in a 2.2 /is window. All PMT pulses in a 2.2 /is time 
window are conserved in a buffer. If a trigger occurs, 
all hits with charges above a given threshold are saved 



5L,3pe MC nuclearites 
Entries 17G905 
Overtiow 1 .07e-03 
5L,3pe MC 
Entries 4. 14545e+06 




1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 

dt (ns) 

Fig. 2: Snapshot duration distributions for simulated nu- 
clearites and muons, broader and narrower curves respec- 
tively. The distributions are normalized to unity, and were 
obtained for a 5 line ANTARES configuration, DT only 
and a 3 pe threshold. The vertical line marks the optimal 
CI cut for this configuration. 



for the actual time window, together with the previous 
buffered information. If, in the next time window another 
trigger fires, the data are merged together, so the minimum 
duration of an event "snapshot" is 4.4 /is. 

The data recorded in 2007 and 2008 were obtained 
during ANTARES completion, so they refer to various 
detector configurations. The variations in bioluminescence 
imposed also different LO threshold values. In order to 
maximize the analysis efficiency, the ANTARES approach 
is to start with a "blind analysis". This consists in defining 
the search strategy using Monte Carlo simulations, and 
validating them on a small fraction (15%) of the available 
data. 

The search for nuclearites in ANTARES is based on the 
fact that the duration of a nuclearite event should be 
much longer (about 1 ms) than in the case of a muon. 
Furthermore, the luminosity produced by the passage of 
a nuclearite is expected to be orders of magnitude larger 
than from muons. The isotropic light emitted along the 
heated path should induce a long succession of fake muon 
signals, so we can use the same triggers as those im- 
plemented for relativistic ionizing particles. Monte Carlo 
simulations have been performed for different nuclearite 
masses, in all ANTARES configurations during 2007 and 
2008. The background level was extracted from real runs. 
The simulated events were processed using the DT and CT 
(for the time period in which they were active) triggers. 
Triggers select only the hits compatible with the passage 
of a muon, so nuclearite events may expand on multiple 
adjacent snapshots. Downwards going muons are the main 
source of background for nuclearite searches, so we inves- 
tigated the snapshot duration dt for both simulated muons 
and nuclearites, defined as the time difference between the 
last and the first LI hits that produce a trigger Muons 
were simulated using the MUPAGE code [7|. Fig. |2] 
shows a comparison between the distributions of snapshot 
durations for simulated nuclearites and muons. This result 
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corresponds to a 5 line ANTARES configuration, DT 
trigger and a 3 pe LO threshold. The vertical line (at a 
snapshot duration of 3250 ns) represents the optimized cut 
for this configuration. All configurations were analyzed in 
a similar way, and optimal cuts (ranging between 2500 and 
4750 ns) were defined as to minimize the upper flux limits 
obtainable if no signal is found. In the following we refer 
to those first level selection, based only on the snapshot 
duration, as "CI" cuts. We tested the CI cuts on 15% of 
real data, selected to reproduce the proportion of different 
experimental configurations. 

For the few events passing the CI cut, the evolution 
in time of the event charge barycenter was studied. As 
the light emitted by nuclearites is isotropic, the charge 
barycenter offers an estimation of the position of the 
source at a certain moment. All events were consistent 
with static light sources, so they were interpreted as 
bioluminescence bursts. Characteristic to those events is 
that they were all single-snapshot events. In order to 
reduce the bioluminescence background, we introduced a 
second level cut C2, to be applied only to single-snapshot 
events, requiring an event duration twice the CI cut value. 
With this second cut, no event in the 15% data sample 
survived. 

The efficiencies of the cuts for different nuclearite masses 
were computed and are presented in Table 1. They are 
defined as the fraction of the nuclearite events detectable 
after filtering the data. 

In terms of nuclearite mass, the threshold of our search 
(defined as the minimum mass of a nuclearite entering the 
Earth atmosphere with (i — 10~^) is about 3 x 10^^ GeV. 
This is in agreement with the results presented in Fig. [T| 
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Fig. 3: The ANTARES 90% CL upper limit for a flux of 
downgoing nuclearites, obtained from the 2007 and 2008 
data. The final limits reported by MACRO and SLIM are 
also shown. 



such slow moving objects should have masses larger than 
some 10^2 GeV |5|. 

We recall that MACRO was an underground detector op- 
erated in the National Gran Sasso Laboratories of INFN, 
and SLIM was a passive detector operated at high altitude, 
at the Chacaltaya Cosmic Ray Laboratory. 

The shape of our limit reflects the narrowing around the 
vertical direction of our solid angle acceptance as the 
nuclearite mass becomes smaller The ANTARES preUmi- 
nary limit improves significantly the MACRO result in all 
the common mass range; the SLIM limit remains better 
for nuclearite masses lower than about 5 x 10^"^ GeV. 

V. Conclusions 



IV. Results 

After unblinding, we performed the analysis of the ex- 
perimental data collected during the 2007 and 2008 runs. 
The analysis followed the pre-defined strategy: data are 
filtered applying the corresponding CI cuts, or the C2 cuts 
for the single-snapshot events. No residual contamination 
was found from downward going atmospheric muons after 
this selection. 

Very few events survived our cuts, but they are not com- 
patible with the hypothesis of slowly downgoing particles; 
they might be safely interpreted as due to biolumines- 
cence. Their topology was investigated using the event 
charge barycenter versus time correlation. This approach 
is equivalent to a fast but less accurate geometrical recon- 
struction of the event. Consequently, we could compute a 
90% confidence level upper limit for a flux of downward 
going nuclearites, from the ANTARES 2007 and 2008 
data. Our result is presented in Fig. [3] compared with the 
MACRO limit |8|, as updated in [9J, and with the SLIM 
limit |10|. All limits in Fig. [3] are obtained in the same 
hypothesis: nuclearites are supposed to enter the Earth 
atmosphere with a velocity of 300 km/s, and refer only 
to downgoing nuclearites. In order to traverse the Earth, 



We have presented the strategy to search for downward go- 
ing nuclearites as well as the preliminary results obtained 
from the 2007 - 2008 data analysis with the ANTARES 
neutrino telescope. The results discussed were obtained 
after unblinding the data corresponding to the detector 
configurations during the 2007 and 2008 data tacking. 
No candidate event survived all criteria; the 90% C.L. 
ANTARES upper flux limit improves significantly the 
MACRO result ID, In the lowest mass region, the 
SLIM limit |10l still remains the most stringent published 
result. 

The same analysis procedure will be applied for the 
ANTARES data recorded after the completion of the 
detector (2009 and later on). 
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Abstract. The composition of the cosmic rays in the 
knee region (ss 10^ TeV/nucleus) of the all particle 
spectrum is considered to be the result of the particle 
acceleration and propagation from the astrophysical 
sources. The steeply falling cosmic ray spectrum makes 
a direct measurement of the composition difficult, 
but it can be inferred from the measurements of the 
showers generated by the interaction of the primary 
cosmic ray with the Earth atmosphere. In particular 
the characteristics of the muon bundles produced in 
the showers depend on the primary cosmic ray nature. 
The ANTARES telescope is situated 2.5 km under the 
Mediterranean Sea off the coast of Toulon, France. It 
is taking data in its complete configuration since 2008 
with nearly 900 photomultipliers installed on 12 lines. 
The trigger rate is a few Hz dominated by atmospheric 
muons. A method using a multiple layered neural 
network as a classifier was developed to estimate 
the relative contribution of proton and iron showers 
from the energy and multipUcity distribution of the 
muon tracks reaching the ANTARES detector. The 
performance of the method estimated from simulation 
will be discussed. 

I. Introduction 

Although the main goal of ANTARES telescope is to 
look for high energy neutrinos coming from the deep 
space, it also provides us opportunities to study cosmic ray 
physics. One of the most important topics is to distinguish 
the different chemical compositions around knee region 
of its spectrum. The cosmic ray spectrum is known as 
a power law with power index about -2.7 up to few 
PeV. Then the slope changes to -3.1 until the energy 
around 4x10^® eV [IJ. The knee was first observed by 
the MSU group in 1970s, then has been confirmed by 
many groups afterwards. The origin of the knee could be 
generally summarized into either astrophysical origin or 
particle physics origin. Nevertheless it is still a puzzle and 
generally believed to be a key issue to the problem of the 
origin of galactic cosmic rays. 

Most people attributed the knee to the sudden reduction in 



Galactic trapping efficiency. A popular explanation is that 
the knee is associated with an upper limit of acceleration 
energy by galactic supernovae. Another popular scenario 
is the leakage of particles from the Galaxy, since the 
Larmor radius of a proton in the galactic magnetic field 
increases with its energy and finally exceeds the thickness 
of the galactic disk. Additionally, there is a minority of 
theorists who proposed that the knee is due to a single, 
recent and local supernova remnant (SNR) or a rapidly 
rotating pulsar interacting with radiation from its parent 
SNR|2|. 

If the knee is caused by the maximum energy attained 
during the acceleration process or it is due to leakage from 
the Galaxy, the energy spectra for individual elements with 
charge Z would exhibit a cut-off at an energy = Z x 
E^, where E'p is the cut-off energy of protons. The sum 
of the flux of all elements with their individual cut-off 
makes up the all-particle spectrum. In this picture the knee 
is related to the proton cut-off and the steeper spectrum 
above the knee is a consequence of the subsequent cut- 
off of heavier elements, resulting in a relatively smooth 
spectrum above the knee [3|. 

II. Composition Model 

In the so-called polygonato model, the general form for 
the flux of primary nuclei of charge Z and energy Eq is 



- A7 



where the transition energy Etrans could be determined 
according to three different scenarios. The parameter Ec 
determines the smoothness of the transition, and 7^ is the 
hypothetical slope beyond the knee. Etrans is the cut- 
off energy. Under the polygonato model, three different 
scenarios are proposed : 

• Rigidity dependent: from the astrophysical point of 
view, a rigidity dependent cut-off Etrans = Ep ■ Z 
is the most likely description if we take into account 
the acceleration and propagation of the cosmic rays. 

• Mass dependent: this model predicts that the change 
in the power law index depends on the mass Etrans 
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= Ep ■ A instead of the charge. This scenario 
leads to a steeper energy spectrum after the cut-off. 
However, the sharp cut-off would be hard to explain 
on astrophysical reasons. Maybe a nearby source or 
a new type of interaction in the atmosphere could 
yield such cut-off |3|. 
• Constant composition: Etrans = Ep. The knee is 
explained by a common steepening in the energy 
spectrum; it occurs for all the particles at the same 
energy. 

So far, the best measurement of the compositions around 
the knee region was done by KASCADE f4]. The mea- 
sured primary energy spectra show that the knee in the 
all particle spectrum is due to a steepening of the light 
elements spectra. 

III. Analysis Strategies 

The ANTARES detector is located at 40 km off the 
coast of Toulon, France, at a depth of 2475 m in the 
Mediterranean Sea. It consists of 12 flexible strings, each 
with a total height of 450 m, separated by about 60 m. 
They are anchored to the sea bed and kept near vertical by 
buoys at the top of the strings. Each string carries a total of 
75 10-inch Hamamatsu photo multipliers (PMTs) housed 
in glass spheres, the so-called optical modules (OMs) [5J. 
The OMs are arranged in 25 storeys (three optical modules 
per storey) separated by 14.5 m. The detector started 
taking data in 2007 and was fully completed in 2008. 
Since ANTARES is deeply buried under the sea, only 
the muon components from the air shower will survive at 
detector level. The muons will emit Cherenkov radiation 
only when passing through the sea water, which can be 
detected using photomultiplier tubes. To be registered by 
the ANTARES detector, muons have to travel at least 
2.5 km of sea water and still be energetic enough to 
trigger the detector. At large zenith angles, the threshold 
increases because of the increasing depth of the sea water. 
The muon bundles properties (such as multiplicity) are 
strongly related to the primary energy and species of the 
nuclei. However, the ANTARES detector cannot resolve 
individual muons from a muon bundle. Hopefully, the 
topology of the hit distributions in space and time could 
give pertinent information about the properties of the 
muons in the bundles. Two useful analysis methods are 
combined and then applied on MC samples. The first one 
is a cluster finding algorithm, and the second one is an 
electromagnetic shower searching algorithm ||6l. 
In both analysis methods, we rotate our software coordi- 
nate system such that the z axis is along the reconstructed 
muon track. We define the plane which is passing through 
the detector center and perpendicular to the reconstructed 
track as detection plane. Fig. [T| shows the snapshot of 
muon bundle patterns on the detector plane from the same 
primary energy of proton and iron nuclei. The projections 
of all the hits positions on the detection plane are also 
calculated, including the time correction information. 




Fig. 1: Typical space distributions of muon bundles on 
the detector plane. The muon bunldes are from primary 
protons (left) and iron nuclei (right) with the same energy. 



In the cluster finding algorithm, the clusters are formed 
if the hits within the cluster fulfill the following three 
conditions: 

• \T, - Tj \ <\ Ti - rj\ ICIUg where and are 
the times for any two hit pairs in the hit cluster 
and Tj are the associated positions of these two hits. 
Ug is the index of refraction in the sea water and C 
is the speed of light. 

• Any two hits within the cluster should be in the same 
or neighboring strings or floors. 

• \T,-T,\< \(z, - Zj)\IC + d X tanO, Ic + T,,t, 
where Tg^t are maximum extra time, here we set 20 
ns; Ti and Tj are again time informations of two hits. 
Zi and Zj are the rotated Z positions of the two hits. 
9c is the Cherenkov angle in the sea water 

In order to quantify the cluster patterns from cosmic 
showers initiated by different groups of elements, we try 
to parametrize the hit patterns. There are two kinds of 
parameters in our analysis: 

• Cluster-wise parameters: the parameters which are 
related to each cluster. For instances, the "size" (total 
number of photoelectrons in one cluster). 

• Event-wise parameters: the parameters which are 
related to each event, such as the total number of 

clusters Nduster- 

In total, we have about 25 event- wise and 12 cluster-wise 
parameters. One example of pattern after parametrization 
superimposed with reconstruction position of the muon 
tracks (bundles) is presented in Fig |2] 

On the other hand, the electromagnetic shower searching 
method gives us additional information about the muon 
bundles. The high energy muons suffer from catastrophic 
energy losses. Once it happens, electromagnetic showers 
are initiated either by 7 or e+ and pairs. Instead 
of projecting hits information on the detector plane, we 
project all the hits on the reconstructed axis of the event 
and search for the peaks using the TSpectrum function 
in the ROOT package. The details of the algorithm can 
be found in |6|. Two more parameters are obtained by 
this method. They are the number of showers (Nshower) 
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Fig. 2: The hit patterns projection on the detector plane 
and then parametrized by cluster algorithm. The black 
points are projection of hits. The red circle is the center of 
gravity (COG) of the four clusters. Green points are the 
positions of the muon in the bundles. Empty circle is the 
COG of the muons. The pink triangle is the reconstruced 
muon position. The blue diamond is COG of all the 
clusters. 



and the amplitude of the baseline {Nhaseiine) from the 
algorithm. 

Combining the cluster-finding and shower-finding algo- 
rithms, we have in total 37 parameters. Each parameter has 
a different discrimination power for the chemical species. 
To achieve multi-dimensional comparisons, we use the 
existing package "TMVA" (Toolkit for Multivariate Data 
Analysis with ROOT ) for this analysis 18j. TMVA is 
a toolkit which hosts a large variety of multivariate 
classification algorithms. Training, testing, performance 
evaluation and application of all available classifiers is 
carried out simultaneously. 

Several methods are implemented inside TMVA pack- 
ages. In order to cross check the results, four different 
methods were chosen for this analysis. They are Multi- 
layer Perception (MLP), (MLPBNN), (TMlpANN) and k- 
Nearest Neighbour (k-NN). The performances of multi- 
variate analysis are sensitive to the correlation of the 
training parameters. Thus, it is necessary to reduce the 
un-important parameters. We input each parameter into 
the analysis individually and then we kept only the five 
with the largest discriminating power These five param- 
eters are Nhu, Nduster and Nnpe, representing number 
of hits, clusters, and photoelectrons from cluster-finding 
algorithm. Nghower and Nbase are the numbers of showers 
and the number of photoelectrons of the baseline from EM 
shower-finding algorithm. 

IV. Analysis on MC samples 

A full MC simulation was adopted in this analysis. The air 
showers induced by the primary nuclei with energy rang- 



ing from 1 to 10^ TeV /nucleon and zenith angle between 
0° and 85° are simulated using the CORSIKA (Version 
6.2) [7J and the hadronic interaction model QGSJET.Olc. 
All muons reaching the sea level, with energies larger than 
the threshold energy, are propagated through sea water 
to the detector. At last, muons are transported through 
the ANTARES sensitive volume, Cherenkov light is pro- 
duced and the detector response is simulated. Background 
noises were added afterwards. The standard ANTARES 
trigger is simulated. The muon direction and position are 
reconstructed using a multi-stage fitting procedure, which 
basically maximizes the likelihood of the observed hit 
times as a function of the muon direction and position||9l. 
We further divide the MC samples into 3 independent 
subsets for training, testing and evaluating the TMVA 
analysis. The relative sizes of the three subsets are 40%, 
30% and 30%. Each event is weighted according to 
different models of CR composition. The iron component 
is tagged as signal whereas the protons are tagged as 
background. The "pure" signal distributions are obtained, 
if we assume that all the coming cosmic rays are iron. 
On the other hand, the background samples come from 
assuming all cosmic rays are protons. The distributions of 
Nshower and Nduster assuming the all particle spectrum 
are, respectively, proton and iron, are shown in Fig [3] 
We carefully checked the output of the neural network 
from test samples in order to avoid the over-training 
effects. 



V. Results and Conclusions 

The training and test events fed into the neural network 
were subjected to a series of cuts. The main quality cut 
is the so-called A cut |91. Its purpose is to keep good 
quality on reconstructed events. For testing the method, we 
mixed proton and iron components from our pseudo-data 
set with w 98.2%-1.8%, individually. The output of the 
neural network training on the proton and iron components 
is shown on the Fig |4] The green curve corresponds to 
the true proton events, while the red curve corresponds 
to the true iron events. The blue curve in the plot shows 
the psedo-data events with the 98.2%-1.8% configuration. 
We adjust proton and iron distibutions to fit the psedo- 
data distribution ("Template fitting" in roofit). We found 
that the x^/ndf is 18.07/21. The fitted numbers from the 
two-component model are 98.5±0.7% for proton. With 
this method, further calculating the efficiencies of triggers, 
different cuts and effective areas are necessary in order to 
find the ratio in the original spectrum. Another approach 
is to put different weight in each event according to 
differnt physics models and check the output of different 
physics models from neural work. An example is shown 
on the Fig [5] From the plot, the differences between 
three extreme models are clearly seen. In the future, with 
this method, we are able to compute the compatibilities 
between different realistic physics models and data based 
on the neural work. 
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Fig. 3: The relative distribution of Nshower (left plot) and Nduster (right plot) parameters after the reconstruction 
quality cuts assuming that all the particles from cosmic rays are proton (red points) and iron (blue points), respectively. 
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Fig. 4: Output computed by the neural network in the 
presence of a mixture of proton and iron components 
from the pseudo-data sets. The green curve region (left) 
corresponds to the true proton events, while the red curve 
region (right) corresponds to the true iron events. The 
artificial mixture of pseudo-data is in blue. 



In summary, we have developed a method to estimate the 
ratio of the heavy elements in the triggered cosmic ray 
events based on the information from the muon tracks 
and electromagnetic showers in ANTARES detector. To 
combine all the discrimination powers from established 
multi-parameters needs the help from multi-variate analy- 
sis (neural network). The analysis of ANTARES real data 
with the goal of deriving ratio between different groups 
of elements in the cosmic ray spectra is ongoing and will 
yield results in the near future. 
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Abstract. Amongst the atmospheric muons recorded 
by neutrino telescopes are muons produced by the 
interaction of cosmic gamma rays with the Earth's 
atmosphere. Although they are not numerous, it has 
been suggested that such muons could be distinguished 
by neutrino telescopes from the isotropic background 
by correlating their direction with known sources of 
gamma rays. 

The Antares neutrino telescope is taking data at the 
bottom of the Mediterranean Sea in its full configura- 
tion since May 2008. Its expected sensitivity to steady 
gamma ray sources is discussed, as well as the gamma 
ray induced neutrino contamination of cosmic neutrino 
signals. It is shown that the expected signal from 
steady gamma ray sources is well below the Antares 
detection abiUty, and that gamma rays are a negUgible 
source of atmospheric neutrinos background. 

I. Introduction 

A large number of astrophysical sources are known to 
emit photons above the TeV range. When interacting 
with the Earth's atmosphere, such photons will produce 
electromagnetic air showers which may contain high en- 
ergy muons. Their detection by neutrino telescopes has 
long been proposed and investigated theoretically (see e.g. 
references in Ref. 1 1 1). Although its location more than 
2000 m depth in the Mediterranean Sea is supposed to 
suppress most of atmospheric muons, which are consid- 
ered as background for cosmic neutrino searches, some of 
these gamma ray induced muons may reach the Antares 
neutrino telescope. Its configuration being optimized for 
the detection of upgoing muons |2|, their detection is 
however challenging. 

Even though it may not be competitive with atmospheric 
Cerenkov telescopes and extended air shower arrays, the 
ability of Antares to detect gamma rays and correlate 
them to their source could be used as a calibration tool 
and would be an independant way to check the absolute 
pointing of the detector. 



The following is an update of preliminary results pre- 
sented earlier 11]. In particular, more accurate software 
is used. Most notably, the treatment of muon pair pro- 
duction from photons has been corrected in the CORSIKA 
program jSj, and the treatment of light dispersion in water 
as well as the trigger simulation have been improved in 
the Antares simulation tools. 

II. Monte-Carlo simulation 

Three Monte-Carlo productions of gamma rays have been 
generated in three contiguous energy ranges, following an 
E-'^ flux : 2.5 X 10^ photons between 1 and 100 TeV, 
2.5 X 10^ photons between 100 and 1000 TeV and 5 x 10^ 
photons between 1 and lOPeV. Interactions in the atmo- 
sphere are processed with the CORSIKA program (version 
6.960) L4J, which makes use of the Egs4 code system 
for the electromagnetic interactions 13, and the selected 
hadronic model is Qgsjet |6|. Sources are considered as 
fixed in the sky, with an azimut angle of 0° (North) and 
a zenith angle of 20 or 60° (downward vertical particles 
having a zenith angle of 0°). 

The propagation of muons in water is performed by 
a dedicated program which makes use of the Geant3 
program for the emission of secondary particles |7| and 
Music tables for the propagation of particles itself |8|, 
and which generate Cerenkov photons using photon tables 
previously generated. 

Optical background is added using standard medium qual- 
ity real data, and hits on the detector are selected using a 
standard trigger strategy ll). 

Muon tracks are reconstructed using a standard likelihood- 
based reconstruction algorithm relying on the minimiza- 
tion of hit time residuals 0. A cut is performed on the 
reconstruction quality estimator so that most events with 
an angular resolution better than 2° are selected while 
most of badly reconstructed events are rejected. 

III. Antares effective area 

Dividing the number of surviving events by the incident 
flux, one can compute the effective area of the detector as 
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Fig. 1: ANTARES gamma ray effective area at trigger and 
reconstruction level, for two zenith angles. 



a function of the incident particle energy, which expresses 
the detector efficiency. Figure [T] shows the effective area of 
the ANTARES telescope to gamma rays at standard trigger 
level, which represents all remaining events after applying 
the standard trigger strategy, and at the reconstruction 
level, representing all reconstructed events surviving the 
reconstruction quality cut, for both simulated zenith an- 
gles. 

Although almost three times lower at trigger level, the 
effective area for a 60° zenith angle exceeds the effective 
area for a 20° zenith angle at reconstruction level. The 
reason for this is comes from the granularity of the 
detector, which is more important in the horizontal plane 
than in the vertical dimension : lines are separated by 
about 60 m, while storeys on a line are separated from 
14.5 m. Consequently, although muons are less numerous 
to reach the detector at 60° because of the increased 
distance to travel in water, their Cerenkov light is detected 
by more lines than at 20°. Their direction can thus be 
reconstructed more efficiently. This affects particularly the 
determination of the azimut angle. 

IV. ANTARES SENSITIVITY 

The expected number of events for a few sources of inter- 
est, assuming a 100 % visibility over one year and a fixed 
position in the sky, is presented in Table [l] The sources are 
selected both for their visibility from ANTARES (Table |ll| 
and for their strong flux. The case of a fictional source 
with a very powerful flux (dN/dE = IGOOE:-^ ni-2s-i) 
is also considered. The higher bound extrapolates gamma 
ray fluxes up to lOPeV and is thus quite unlikely since the 
validity of flux parametrizations is limited to a few tens 
of TeV at best and since interactions of ultra-high energy 
photons with photons from the extragalactic background 
light (EBL) limit their range to galactical distances. The 
lower bound is more realistic and does not extrapolate 
above 100 TeV. 



source 


Ntrir 
'trig 




Crab 


4-15 


0.1-0.3 




0.6-2.6 


0.1-0.4 


lES 1959+650 


0.5-25 


0.01-0.5 




0.1-6 


0.01-0.8 


Mkn 421 


1-25 


0.01-1 




0.1-2 


0.01-1 


Mkn 501 


1-100 


0.03-10 




1-50 


0.03-15 


fictional 


820-1860 


22-68 




226-660 


29-92 



TABLE I: Expected number of gammar ay events seen by 
ANTARES at trigger level (Ntrig) and after reconstruction 
(Nrccoi with a cut on the reconstruction quality and recon- 
structed within 2° from source), for one year assuming a 
100% visibility, for a 20° zenith angle (straight font) and 
a 60° zenith angle (italic). The lower bound assumes a 
cutoff at 100 TeV while the upper bound assumes events 
as energetic as lOPeV. 



source 


visibility 


{e) 


Crab 


62% 


51.7° 


lES 1959-^650 


100% 


49.7° 


Mkn 501 


78% 


49.4° 


Mkn 421 


76% 


49.2° 



TABLE II: Actual visibility and mean zenith angle of 
selected sources from the ANTARES down-going muons 
field of view. 



These results show that it is unlikely for ANTARES to 
detect steady gamma ray sources : in the same optical 
background conditions, the number of reconstructed tracks 
surviving the selection cuts within 2° from source in a 
sample from actual data extrapolated to one year is about 
53800 (17900) for a zenith angle of 20° (60°), which is 
still far above the expected number of events detected from 
the fictional powerful source. Reducing the observation 
cone does not improves the signal over background ratio 
enough to compensate the loss of signal it induces. 

Results are reported in Figure |2] as the power law flux 
normalization needed to obtain a 3 or 5 standard devia- 
tions sensitivity (after selection of events) as a function 
of the spectral index. Also shown is the area covering 
Crab flux parametrizations and the position of the fictional 
powerful source. It follows that ANTARES would need a 
source more than three orders of magnitude larger than the 
Crab or alternatively a powerful source with an extremely 
hard spectrum to obtain a reasonable sensitivity, again 
assuming a 100 % visibility and a fixed position in the sky. 
Sources providing such characteristics are very unlikely to 
exist since similar features have yet never been measured 
even for short strong outbursts |10|. 

V. UNCERTAINTIES 

This study suffers many sources of systematic uncertain- 
ties. Given the low sensitivity of ANTARES to gamma 
ray induced muons, the various contributions are not 
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Fig. 2: Antares one year sensitivity to steady gamma ray 
sources, as a function of flux normalization and spectral 
index assuming a power law spectrum. 



computed in a quantitative way. Many of the systematics 
have a negligible effect on the simulation. This is for 
instance the case of the Earth magnetic field, the LPM 
effect 1 11] or the preshower effect fT2). 

Important systematics are mainly due to the simulation 
of interactions in the atmosphere. Charm photoproduction 
cross-section remaining unknown at very high energy, it 
is not included in CORSIKA. It is small but increases with 
the energy and may not be negligible above a few tens of 
TeV 1 13]. Furthermore, the charm carries a large fraction 
of the incident photon energy and has thus a high proba- 
bility to reach the detector, which makes this contribution 
difficult to estimate. The lack of charm production is hence 
an important flaw of this simulation. The hadronic model 
also introduces important uncertainties. Switching the 
hadronic model to VENUS [14] can increase the number 
of muons with an energy higher than 500 GeV at sea 
level by about 10% and the number of events producing 
such muons by 7 %. Finally, the photoproduction cross- 
section used in CORSIKA is quite conservative, and the 
muon yield might be increased by about 10 % by other 
realistic scenarios 1 15J . According to these effects, the 
present simulation may be thought as pessimistic. 

On the other hand, the description of the sea water optical 
properties used for the light simulation in water is not the 
most accurate and overestimates the number of detected 
Cerenkov photons and hence the number of reconstructed 
muons. The number of detected muons could be increased 
by as much as 20 % with regard to reality. The geometry 
of the detector also affects the detection of events : the 
azimuth angles for which several detector lines are aligned 
are favoured, and the number of reconstructed events may 
be increased by a factor of two in the extreme cases. 

Furthermore, a neutrino telescope is sensitive to variable 
sources of uncertainties. The muon yield itself depends 
on atmospheric variations : the higher the pressure and 



temperature, the lower the muon yield, by a few per- 
cents. This affects both signal and muon background. The 
sea water optical properties and most of all the optical 
background due to radioactive decays (mostly "^^K) and 
bioluminescence are also strong sources of variability. 
The present simulation is valid for standard conditions of 
optical background qualified for data analysis, but small 
degradation or improvement of these conditions may have 
an effect as large as 20 % on the number of reconstructed 
tracks. Finally, the accuracy of the software used for the 
simulation of light propagation and detection has been 
greatly improved since versions used in this simulation. 

Last but not least, in order to ensure statistical consistency 
of the simulation, the source zenith angle is fixed. In 
reality ANTARES sensitivity to gamma ray sources is 
obviously affected by their motion in the sky. 

VI. Discussion 

Although this simulation suffers large uncertainties, it is 
unlikely that the results would vary by an order of mag- 
nitude or more. It is thus clear from this study that steady 
gamma ray sources are out of reach of the ANTARES 
neutrino telescope under the conditions discussed here. 

This analysis makes use of standard triggers and recon- 
struction algorithms, optimized for detection of upgoing 
neutrinos. Dedicated strategies and background discrimi- 
nation methods would have to be developped to improve 
significantly the sensitivity : directional trigger, use of the 
muon-poor electromagnetic showers characteristic, muon 
pair discrimination. . . 

In addition, the present study assumes that the gamma 
ray source flux at very high energy can be parametrized 
by a power law (possibly with an exponential cutoff). A 
significant deviation of this assumption at higher energies 
would have a strong impact on the presented results. 

VII. Gamma-induced neutrino background 

Gamma rays interacting with the atmosphere may also 
produce neutrinos. These so-called atmospheric neutrinos 
form an irreducible background for neutrino telescopes, 
since there is no way to distinguish them from cosmic 
neutrinos. Furthermore, such atmospheric neutrinos are 
localized from the direction of gamma ray sources, which 
are cosmic neutrino emitter candidates. 

The present study is the opportunity to estimate the 
number of upgoing atmospheric neutrinos polluting the 
signal from the direction of a gamma ray source. This is 
simply done by multiplying ANTARES neutrino effective 
area |16| by the gamma ray induced neutrino flux at sea 
level. 

It is found that the expected number of reconstructed 
events is smaller by several orders of magnitude than the 
expected cosmic neutrino signal for all simulated sources. 
The contamination induced by steady gamma ray sources 
to a potential source of neutrino can thus be considered 
as negligible. 



VIII. Conclusions 

It has been demonstrated through a complete Monte-Carlo 
simulation that the observation and identification of very 
high energy gamma rays from steady sources is out of 
reach of the Antares neutrino telescope with its standard 
trigger and reconstruction strategies. It is also clear from 
this study that the gamma ray flux of such somces will 
not interfere with the identification of a possible neutrino 
sources. 
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Abstract. The primary goal of the ANTARES tele- 
scope is the detection of high energy cosmic muon 
neutrinos. The neutrinos are identified by the upward 
going muons that are produced in charged current 
neutrino interactions in the vicinity of the detector. 
The Cherenkov light produced by the muons in the 
detection volume is measured accurately by an array 
of photosensors. Muons that are going downward are 
background for neutrino searches. These muons are 
the decay products of cosmic-ray collisions in the 
Earth's atmosphere. The energy loss in water of a 
muon with an energy above a TeV is characterized by 
discrete bursts of Cherenkov light originating mostly 
from pair production and bremsstrahlung (electro- 
magnetic showers). This paper presents a method to 
identify and count electromagnetic showers produced 
by the muons. The method can be used to select a 
sample of highest energy muons with the ANTARES 
detector. 

I. Introduction 

The ANTARES neutrino telescope is located on the bot- 
tom of the Mediterranean Sea, roughly 40 km off the 
French coast. The main objective is the observation of 
extraterrestrial neutrinos. Relativistic charged leptons pro- 
duced by neutrino interactions in and around the detector 
produce Cherenkov light in the sea water This light is 
detected by an array of photomultiplier tubes (PMTs), 
allowing the muon direction to be reconstructed. The 
muon energy loss can be estimated from the sum of the 
measured number of photoelectrons. 

The detector is installed at a depth of 2475 m and 
consists of twelve vertical lines approximately 450 m long 
equipped with a total of 885 PMTs. The lines are separated 
from each other by about 65 m, anchored to the sea bed 
by a dead weight and held taut by a buoy located at the 
top. The instrumented part of the line starts 100 m above 
the sea bed, with 25 storeys every 14.5 m along the line. 
A storey consists of three PMTs pointing downward at 
an angle of 45° with respect to the vertical direction. 



in order to maximise the detection efficiency of upward 
going tracks. 

ANTARES is operated in the so called all-data-to-shore 
mode: all pulses above a threshold (typically 0.3 photo- 
electrons) are digitized off-shore and sent to shore to be 
processed in a computer farm. This computer farm applies 
a set of trigger criteria in order to separate muon-induced 
Cherenkov light from background light. The main sources 
of the background light are the decay of ""^K nuclei and the 
bioluminescence of organisms in the sea water A detailed 
description of the ANTARES detector is given in [1]. 

Although ANTARES is optimised for upward going par- 
ticle detection, the most abundant signal is due to the 
atmospheric downward going muons produced by the 
interaction of primary cosmic-rays in the atmosphere. 
Being the most penetrating particles in such air showers, 
muons with enough energy can reach the detector and are 
reconstructed by the detection of the Cherenkov light they 
emit when propagating through water The ANTARES 
detector has the capability to follow highly energetic 
muons over a few hundred metres. 

The processes contributing to the energy loss of a 
muon in water include ionisation, e+e~ pair produc- 
tion, bremsstrahlung, and photonuclear interactions. Be- 
low 1 TeV, the muon energy loss is dominated by the con- 
tinuous ionisation process. Above 1 TeV, the muon energy 
loss is dominated by pair production and bremsstrahlung 
||2l, which are radiative processes classified as electro- 
magnetic showers. They are characterized by large energy 
fluctuations and discrete bursts along the muon track. The 
average muon energy loss per unit track length due to 
these electromagnetic showers increases linearly with the 
energy of the muon allowing its energy to be determined. 
Counting electromagnetic showers along muon tracks 
gives an estimate of the muon energy |3| and can help 
in designing a better energy reconstruction algorithm. A 
similar measurement technique as the one presented in 
this article has been published recently by the Super- 
Kamiokande experiment (3). 
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Fig. 1: Schematic view of the Cherenkov light detection. 
The thick line represents the muon trajectory, the thin line 
the path of the Cherenkov light and the thin dashed line 
the path of the shower light. The muon goes through a 
reference point (r, z, t). The Cherenkov light is emitted at 
an angle 0cK with respect to the muon track at point (^p^ 
and is detected by a PMT as a hit at point (r^, z^, ti). The 
shower light is emitted at point Q and is detected by the 
same PMT at different time. 



11. Algorithm for shower identification 

The purpose of the electromagnetic shower identification 
is to distinguish Cherenkov photons emitted continuously 
along the muon track, from the Cherenkov photons in- 
duced by electromagnetic showers. Because of the radia- 
tion length in water (Xq — 35 cm), these showers never 
extend more than a few metres and can be considered 
point-like light sources in the ANTARES detector. The 
showers can be identified from a localised increase of the 
number of emitted photons above the continuous baseline 
of Cherenkov photons emitted by a minimum ionizing 
muon. 

In what follows, a hit is a photomultiplier signal exceeding 
a threshold of 0.3 photoelectrons. The shower algorithm 
consists of two steps. The first step is to identify and 
reconstruct downward going muon track candidates. In 
the second step, a distinct shower candidate is identified 
by looking for an accumulation of hits on a point along 
the muon path. The criteria to isolate the accumulation 
are defined from a simulation code based on Corsika |5|. 

A. Simulation 

Cosmic -ray interactions in the atmosphere including 
shower development were simulated for primary energies 
between 1 TeV and 10^ TeV, and incident angles between 
zero (vertical) and 85 degrees with Corsika. The primary 
cosmic-ray composition and flux model considered was 
a simplified versiorQ of the Horandel model fSJ. The 
hadronic interaction model chosen was QGSJET fT\. The 
result of the Corsika simulation consists of muons with 
their positions and times and kinematic vectors on the 
surface of the sea. The muon propagation through water, 
the discrete energy losses at high energies, the Cherenkov 
light production and propagation, including scattering, and 
the response of the detector was simulated using the 
ANTARES KM3 code |8|. KM3 uses tables generated 
from a simulation with GEANT 3 which parametrise the 



arrival time and the amount of light detected by indi- 
vidual PMTs. These tables take into account the angular 
dependence of the acceptance of the PMT as well as 
the measured properties of the water at the ANTARES 
site. The muon propagation uses the MUSIC code l9l 
and is done in short steps (1 m). If the energy loss of 
the muon over the step exceeds a threshold (1 GeV), an 
electromagnetic shower is initiated and shower photons 
are emitted, otherwise if the energy loss of the muon over 
the step is below the threshold muon Cherenkov photons 
are emitted. The optical background was assumed flat at 
a rate of 60 kHz on each photomultiplier 

B. Algorithm 

Muon events are reconstructed by using an existing al- 
gorithm ifTOl . ifm which provides an estimate of the 
direction and position of the muon at a given time. Mea- 
sured hit times are compared to the expected arrival time 
of direct Cherenkov photons. The expected Cherenkov 
photon arrival time for each hit i is calculated as 
(see Figure [TJ: 



f. 



CK 



t + 



tan^CK^ c sin^cK 



, (25) 



where t is the time where the muon passes point (r, z), c 
is the speed of light in vacuum, n is the refraction index 
of water {n is about 1.38), 6'ck is the Cherenkov angle 
for a relativistic muon in water (f?cK = 42°) and Vi is 
the perpendicular distance between the muon trajectory 



and the PMT. Equation ( 25 1 separates the direction along 



the track and the direction perpendicular to the track. The 
direction along the track (z-coordinate) is given by the 
muon momentum vector The direction perpendicular to 
the track (r-coordinate) is given by the photon momentum 
vector in water Hits too far in time from the expected 
muon hit time —tmin <U — tf^ < t,nax, are assumed to 
be background hits and are rejected, whereas direct hits 
have a roughly Gaussian distribution at zero with a width 
of 20 ns. The value for tmin — 20 ns is given by mainly 
the dispersion of light in water and the timing resolution 
of the PMT, whereas t,nax — 200 ns is defined by the 
value where the number of signal hits approaches the level 
of background hits. Furthermore, the above defined time 
interval is subdivided into two intervals. The early interval 
contains mostly muon Cherenkov photons and is given 
by \ti — tp^l < train- The Cherenkov photon emission 
position (Cf'^) along the muon trajectory is given by: 



CK 



tan 6'ck 



(26) 



The late interval is defined by tmin < ti—tf^ < t^ax and 
contains mostly electromagnetic shower photons. These 
shower photons may not necessarily be emitted at the 
Cherenkov angle from the muon track. The emission angle 



1 . The primary composition of the flux is subdivided into only 
five mass groups, namely proton, helium, nitrogen, magnesium 
and iron. 
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Fig. 2: The number of detected photons projected along 
the muon trajectory for an atmospheric muon event 
with two electromagnetic showers. The peaks indicate 
the shower positions on the muon trajectory, with the 
reconstructed shower positions indicated by the triangles. 



is left as a free parameter and, with the photon emission 
taking place at Q (see Figure [T}, the hit time is given by 



U=t 



(27) 



Equation (27i can be solved for C^f, yielding: 



2A 



where 



A = 

Bi = 2{r?Zi~ z- c{ti-t)), 

C, = c^[U-tf + 2cz{U-t)+z'^ -v?{r1 + z'^^). 

All calculated Cp^^ Ct ^^'^ C[ positions along the muon 
track are collected in a one-dimensional histogram. As an 
example of such a histogram. Figure [2] shows all the calcu- 
lated photon emission positions along the muon trajectory. 
The electromagnetic showers are identified by an excess of 
photons above the continuous baseline of photons emitted 
by a minimum ionizing muon. Two excesses are visible 
that can be attributed to the two reconstructed showers. 
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Fig. 3: Muon event rate as a function of the shower 
multiplicity for reconstructed showers and for generated 
showers which have shower photons detected on at least 
five different storeys for the Corsika simulation. 



> 

Z 5- 

0) 
Q. 

"> 

o 

in 3: 



Detectable showers 
Reconstructed showers 



0) 

E 
3 



2r 



1.5 2 



2.5 3 3.5 
log10(E 



4 4.5 
, [GeV]) 



Fig. 4: Average number of detectable showers which have 
shower photons detected on at least five different storeys 
per atmospheric muon event as a function of the muon 
energy. Also shown is the average number of reconstructed 
showers. 



III. Selection and results 
The selection and performance of the shower identification 
algorithm has been studied and validated with a sample 
of simulated atmospheric multi-muons with constant back- 



ground light as described in section II-A 



A. Selection 

The shower algorithm makes use of tracks fitted with the 
muon reconstruction algorithm described in ifTOl with two 
additional criteria. These criteria require the tracks to be 
traced for at least 125 m and to have a minimum of 
twelve hits used in the track reconstruction. These selected 
tracks account for around 65% of all reconstructed tracks. 
The advantage of the selected tracks is not only that the 
direction of the tracks is better reconstructed, but also that 



the tracks are long enough to have a high probability to 
emit showers. 

The parameters in the shower reconstruction algorithm are 
the width and the height of the peak. The analysis has 
been tuned to select showers with a high level of purity, 
possibly at the expense of efficiency. For each selected 
peak, the number of hits is integrated in a ±5 m interval 
around the peak center. Only peaks having at least 10 hits 
over the muon-track Cherenkov photon baseline in this 
interval of 10 m are selected. The number of baseline hits 
is defined as the average density of hits along the track 
times the interval width of 10 m. In addition, in order to 
suppress fake identified showers, hits from at least five 
different storeys are required in each peak. 
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B. Results 

The main result of the shower reconstruction algorithm is 
the shower multiplicity per atmospheric muon event. The 
atmospheric muon events are usually muons in a bundle 
with an average multiplicity of around 3.3. Figure [3]shows 
the event rate as a function of the number of generated 
showers with shower photons detected on at least five 
different storeys. Also shown is the number of recon- 
structed showers selected with the cuts mentioned in sec- 
tion |III-A| The average shower reconstruction efficiency 
over all shower energies is around 4%. The reconstruction 
algorithm starts to be efficient for shower energies above 
300 GeV. 

Figure |4] shows the number of generated showers with 
shower photons detected on at least five different storeys 
per atmospheric muon event as a function of the muon 
energy. The muon energy refers to the muon with the 
largest energy in the bundle. The number of generated 
shower and reconstructed shower increases as a function 
of the muon energy. Figure |5] shows the energy distribu- 
tion of all muons as well as the one with at least one 
reconstructed shower The muons have an average energy 
of 1.2 TeV, whereas muons with at least one identified 
shower have on average 2.5 times higher energy. 

IV. Conclusion 
A method to identify showers emitted by atmospheric 
muons has been applied to simulated data of the 
ANTARES detector The main differences between the 
shower light and the muon Cherenkov light are that the 
shower light is produced on a point along the muon path. 
The essential element of the algorithm is that the selection 
of shower photons is reduced to a one-dimensional prob- 
lem. The performance of the identification algorithm has 
been validated using a sample of simulated atmospheric 
muon events. 

The aim of this proceeding has been more to demonstrate 
the capability to detect electromagnetic showers than to 
make precise measurements. A more elaborated analysis 
is needed, in order to use the number of electromagnetic 
showers as a robust estimate of the muon energy. More- 
over, the method discussed here for selecting showers 
emitted by downward going muons could be used also for 
upward going muons with the main purpose of selecting 
the highest energy upward going muons. 
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Fig. 5; The solid line shows the energy distribution of 
all reconstructed atmospheric muons. The dotted lines 
represent the muon energy distribution for the muons with 
at least one reconstructed shower. 
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Abstract. The ANTARES neutrino telescope is located 
at a depth of 2475 m in the Mediterranean Sea. Its 
main objective is the observation of extraterrestrial 
neutrinos. Relativistic charged leptons produced by 
neutrino interactions in and around the detector pro- 
duce Cherenkov Ught in the sea water detected by a 
three dimensional grid of photomultiplier tubes. The 
propagation of Cherenkov Ught depends on the optical 
properties of the sea water and their understanding 
is crucial to reach the optimal performance of the 
detector. This paper presents the measurements made 
between 2008 and 2010 of the light velocity and atten- 
uation length at the ANTARES site with a system of 
Ught sources (LEDs and laser) at different wavelengths 
(between 400 nm and 532 nm). The time variability of 
the optical properties are presented and the derived 
values are compared with theoretical predictions. 

I. Introduction 

The ANTARES neutrino telescope |[T1 is located at the 
bottom of the Mediterranean Sea at a depth of 2475 m, 
roughly 40 km offshore from Toulon in France. Sea water 
is used as the detection medium for the Cherenkov light 
emitted by relativistic charged particles resulting from 
interaction of neutrinos around or inside the detector The 
particle direction is reconstructed from the arrival time 
of detected photons through the array of photomultiplier 
tubes (PMTs). 

The measurement of the refractive index is performed 
with a pulsed light source shining through water and the 
time of light distributions of photons detected by PMTs at 
different distances from the source. The attenuation length 
is measured by the amount of light detected by these 
PMTs. The optical properties are measured at wavelengths 
between 400 nm and 532 nm and are compared with 
theoretical predictions. 

A precise measurement of the optical properties minimizes 
the uncertainty on many physical results as seen in 12)1 . 
Moreover, the optical properties may change in time due 



to sea current. Several measurements of those presented 
in this paper have been performed in the past 13], ID, ||5l. 

II. Experimental setup 

The ANTARES detector consists of a three dimensional 
array of 885 PMTs arranged in twelve vertical approxi- 
mately 450 m long lines. Along each line with a vertical 
separation of 14.5 m, PMTs are grouped in triplets and 
oriented with their axis pointing downward at an angle 
of 45° with respect to the the vertical line direction. The 
horizontal separation between lines is about 70 m. 

The PMTs are sensitive to single photons in the wave- 
length range between 350 nm and 600 nm. They have 
a peak quantum efficiency of about 25% between 350 
nm and 450 nm. The PMT measures the arrival time and 
charge amplitude of the detected photons. 

The Optical Beacon system consists of a series of pulsed 
light sources distributed through the detector The primary 
aim of the Optical Beacon system is the time calibration 
between the PMTs to reach the best angular resolution of 
the detector. In addition the Optical Beacon system can 
be also used to determine the optical properties of water. 

There are four LED Beacons per line and one laser 
Beacon at the bottom of the central line. One LED Beacon 
contains 36 individual LEDs distributed over six vertical 
faces shaping an hexagonal cylinder. On each face, five 
LEDs point radially and one upwards. All the LEDs emit 
light at a nominal wavelength of 470 nm except two LEDs 
located on the lowest Beacon of line 12 which emit light 
at nominal wavelength of 400 nm. 

The LEDs emit light with a maximum intensity of ^160 
pj and a pulse width of ^4 ns (FWHM). The laser is a 
more powerful device and emits shorter pulses than the 
LEDs. The laser emits pulses of light with a maximum 
intensity of (^1 pJ) and pulse width of ^0.8 ns (FWHM) 
at a nominal wavelength of 532 nm. The LEDs and laser 
flash at a frequency of 330 Hz. Further details about the 
Optical Beacon system can be found elsewhere L6J, UJ, 

m. 
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Fig. 1: Arrival time distribution of LED light detected by a 
PMT at around 100 m. The inset shows a zoom around the 
signal region. Superimposed is the result of the fit to the 
convolution of a gaussian distribution and an exponential 
distribution (see text). 



The light spectrum of the three sources (two LEDs with 
nominal wavelengths of 470 run and 400 nm and one 
laser with nominal wavelength of 532 nm) were mea- 
sured using a high resolution calibrated spectrometer from 
Ocean Optics HR4000CG-UV-NIR. The spectrometer was 
crosschecked with the Green Nd-YAG laser (532 nm). 
The measured peak wavelengths of the LEDs in pulsed 
mode operation are 468.5 ± 14.4 nm and 403.1 ± 6.9 
nm respectively. 

III. Data acquisition and data analysis 

In the special calibration data taken with Optical Beacon 
system, the emission time and the position of the isotropic 
light flash, as well as the arrival time and the position 
when the light reaches the PMTs are known. From the 
time and position difference between the PMTs the re- 
fractive index is measured, whereas the amount of light 
collected by these PMTs gives an information about the 
attenuation length. 

A. Data acquisition 

The various measurements of the water optical properties 
were taken from May 2008. Only data with stable back- 
ground rates and below 100 kHz have been analyzed. The 
special calibration data duration takes only few minutes. 
One single upward looking LED of the lowest Optical 
Beacon in the line emits more than 10^ light flashes 
towards the PMTs above. The photons are collected by 
the PMTs in the line. Figure [T] shows the arrival time 
distribution of photons for a LED Beacon illuminating a 
PMT located at 100 m. The small tail of delayed photons 
is due to the light scattering between LED and PMT. The 
zero time is defined by the illumination time of the LED. 
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Fig. 2: The figure shows the arrival time as a function 
of the distance between the LED and the different PMTs. 
The slope of a linear fit to the arrival time versus distance 
gives the inverse of the light velocity. 

The flat distribution before and after the peak shows the 
detected optical background. 

B. Data analysis 

Since the light path used in the analysis (few 100 m) 
has the same order as the scattering length the scattering 
has in some way to be taken in account. The arrival 
time distributions are fitted to a function which is the 
convolution of a gaussian distribution and an exponential 
distribution |9|, |10|. The gaussian distribution models the 
transit time spread of the PMTs, the time width of the 
optical sources and the effect of the chromatic dispersion 
in water, while the exponential distribution models the 
scattering of photons in water. The fit function is given 
by 

/W = & + /^-e-^.Erfc(^i=(^-^)), (28) 

where t is the arrival time of the photons and the fit 
parameters are the background (6), the height of the 
fit function (/i), the mean and width of the gaussian 
distribution {pL, a) and the exponential constant (t). The 
Erfc(t) is the complementary error function distribution 
Erfc(t) — ^ e* dt' . An example of such a fit is 
shown in the zoom of Figure [T] The fit is made in the 
range from -100 ns to the time of the most populated 
bin plus 20 ns. The arrival time at each PMT is given 
by the fitted mean value of the gaussian distribution. The 
fit is stable with respect to changes in the fit range and 
histogram binning. 

The distances between Optical Beacon and PMT versus 
the measured arrival times are shown in Figure |2] (note 
there can be up to three PMTs at each distance). Due 
to the small wavelength range of the individual LEDs the 
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Fig. 3: Amount of light collected by the PMTs as a 
function of the distance. The fit result of an exponential 
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Fig. 4: The measured refractive index for a total of 42 
runs with an emission wavelength of 470 nm. The result 
to a fit to a gaussian distribution is shown. 



light dispersion is minimized. By selecting only the PMTs 
of the same line as the emission of light the uncertainty on 
the PMT angular acceptance and the PMT positions are re- 
duced. In this analysis, no correction of the PMT positions 
due to line movements have been applied. The minimal 
distance is set to eliminate the PMT which receives too 
much light and have an erroneous time estimation caused 
through the early photon effect (for explanation see |8|) 
. The minimal distance for the fit range is defined as 
the distance where the average collected charge per hit 
(usually referred to as the amplitude) in the signal region 
is below 1.5 photo-electrons. The maximum distance is 
introduced due to avoid noise fluctuations. The signal has 
to be significantly larger than the average background 
(above seven sigmas). The slope of a linear fit through 
the measured points gives the inverse of the measured 
velocity of light in water (wm)- The measured refractive 
index is defined as 



c/v„ 



(29) 



with c — 3 ■ IQ^m/s. The error of the refractive index 
given in Figure |2] is the error estimated by the linear fit. 

As seen in Figure [T] the PMTs perform a time resolved 
measurement of the collected light. The amount of light 
detected depends on the attenuation length, whereas the 
shape of the arrival time distribution of the detected light 
is related to the photon path length distribution of the 
scattered photons. For the attenuation length measurement 
a similar selection criteria as in the refraction index 
measurement is used. An exponential fit to the collected 
charge as a function of the distance is shown in Figure 
[3] for two runs with sources with wavelength of 470 nm 
(tagged as Blue) and wavelength of 400 nm (tagged as 
UV). The attenuation length L is obtained by 



IR^ 



(30) 



where /q is the intensity at the source and / the intensity 
detected by a PMT at a distance R. 



IV. Monte Carlo measurement 

The Monte Carlo simulation takes into account the ge- 
ometry of the detector and the optical properties of the 
water (refraction index, scattering length and absorption 
length). Monte Carlo generated time distributions with 
a large variety of optical properties has been used to 
analyze the time distributions and calculate the refraction 
index and also to check the stability of the analysis 
method. The analysis method was first validated with a 
Monte Carlo sample without scattering at three different 
refractive index. The variation of the absorption length 
between 30 m and 120 m has nearly no influence on 
the refraction index, with variation of the refraction index 
of less than 0.1%. A variation in the scattering length 
between 20 m and 50 m produces an uncertainty of 0.3% 
in the measurement of the refractive index. 

V. Data measurement 

Between May 2008 and March 2010 a total of 42 runs 
were taken with a nominal wavelength of 470 nm, 14 
runs with a nominal wavelength of 400 nm and 13 runs 
with a nominal wavelength of 532 nm and have been 
analyzed according to the methods explained in section 
IIII-BI The measured refraction of the 42 runs are shown in 
the Figure |4]and the mean value is evaluated by fitting the 
distribution with a gaussian. The fitted attenuation lengths 
for some of these runs as function of the data period are 
shown in Figure [5] The fitted values are reasonably stable 
with time. 

Since November 2010 runs with a modified Optical 
Beacon have been collected with light emission at six 
additional wavelengths with nominal peak values of 385, 
400, 440, 460, 505 and 518 nm. First results are expected 
soon. 

These measured refractive index with their systematic 
errors estimated in section IV are shown in Figure [6] 
Also shown is the parametric formula of the refractive 
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Fig. 5: Time evolution of the attenuation length parameter 
for one and a half year of data taking for two different 
wavelengths. 



index. The refractive index at the ANTARES site depends 
on the wavelength, on the temperature, on the salinity of 
the water and the pressure at the depth of the detector. 
The parametric formula for the phase refractive index of 
Quan and Fry [11], based on data from Austin and Halikas 
[121, is modified with appropriate pressure corrections as 
suggested in Q. The phase refractive index for sea water 
as a function of wavelength, temperature (T), salinity (5) 
and pressure (p) is given by 

np(340 < \{nm) < 560, = 3.844, 
T{°C) = 13.2, 200 < piatm) < 240) = 

1.32292 + (1.32394- 1.32292) ^ ~ ■ 

16.2561 4382 



240 - 200 
1.1455 • 10^ 



A 



A2 



A3 



(31) 



where A is the wavelength of light. 

The group refractive index (n) [4J is related to its phase 
refractive index (rip) through 



(32) 



-[^ I ^ drip 

7lr, dX 



The parametrization of the refractive indexes n and rip is 
shown in Figure |6] for the given values of temperature, 
salinity and for a pressure between 200 atm and 240 atm. 
The measurements are in agreement with the parametriza- 
tion of the group refractive index. 



VI. Conclusion 

Pulsed light sources with wavelengths between 400 nm 
and 532 nm shining through sea water and the time of light 
distribution detected by the PMTs at distances between 
few tens and few hundred meters from the source have 
been used to measure the refraction index. Dedicated 
Monte Carlo simulation has been used to validate the 
analysis method and to evaluate the systematics. The data 
results are compatible with the parametrization of the 
group refractive index. 
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Fig. 6: Group and phase refractive index as a function of 
the wavelength for a given temperature and salinity with 
measured data points and its systematical error bars. The 
upper curves (solid) show the group refractive index, the 
lower curves (dotted) give the phase refractive index. 
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Abstract. The ANTARES neutrino telescope is operat- 
ing in the Mediterranean sea in its full configuration 
since May 2008. While designed to observe up-going 
neutrinos, it also records many signals of down-going 
muons produced by the interaction of high energy 
cosmic rays in the atmosphere. The shadowing of 
cosmic rays by the Moon produces a deficit of muons 
coming from this direction at the ANTARES level. The 
observation of this deficit of events can be used to 
check the correct pointing of the detector with respect 
to a known object. 

The strategies to observe this Moon shadow are dis- 
cussed ant the current results are presented. 

I. Introduction 

The ANTARES neutrino telescope is operating in the 
Mediterranean Sea in full 12 lines configuration since May 
2008 [1]. Thanks to the good optical properties of the wa- 
ter of the deep sea, the median angular resolution obtained 
with the current reconstruction algorithm is estimated to 
be 0.4° for cosmic neutrinos with an E^^ flux Q. 

The pointing performance of the detector relies on the 
knowledge of parameters, such as the relative delays of 
the optical sensors within each detector line as well as 
between the lines, the instantaneous shape of the detector 
lines which is influenced by the sea current or the absolute 
orientation of the detector. 

The values of these parameters have been calibrated 
before immersion of the detector lines and are regularly 
measured in situ during operation. The relative time delays 
are measured using a laser beacon and LED beacons, 
and also with physical quantities such as the ^"K de- 
cay or the atmospheric muon tracks [3 |. The alignment 
measurements are performed continuously with acoustic 
emitters and transceivers, tiltmetres and compasses |4|. 
The absolute orientation of the telescope is obtained by 
the triangulation of acoustic signals between anchors of 
the lines and the deployment vessel at the sea surface, 
positioned by GPS [Sj. The uncertainties on these values 
are small enough to ensure the optimal performances of 




Fig. 1 : Illustration of the Moon blocking the cosmic rays, 
thus producing a deficit of muons originating from this 
direction. 



the telescope, but it remains very important to verify the 
correct operation of the detector using a calibration source. 
While no such source exists, some possibilities remain: 

• The observation of a coincidence between an air 
shower at sea level with a surface array and the 
resulting muons at ANTARES level f6\. 

• The observation of the Moon shadow. This technique, 
first proposed in 1957 JS), relies on the absorption of 
cosmic rays by the Moon, as illustrated on figure [T] 
As the cosmic rays are the source of the down-going 
muons at the telescope level, this absorption induces 
a deficit of the number of observed muons in the 
Moon direction. 

II. Data and Monte Carlo 
A. Data set and reconstruction 

The apparent Moon radius is on average Rmoou = 0.259°, 
which is smaller than the resolution of ANTARES, in 
particular for down-going atmospheric muons; thus the 
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Fig. 2: Top: Distribution of the reconstructed zenith of the 
events of the data set used in this analysis: data (black), 
Monte Carlo (red=yeUow+blue). Bottom: Data/MC ratio. 



0.8 



0.6 



0.4 



0.2 + 



/ \ 



0.5 



Moon 
PSF 

PDF Moon 
PDF no Moon 



i"rH--i-.t.j- 



1.5 2 2.5 3 
radial distance (deg) 



Fig. 3: Black: Parametrization of the PSF of the events at 
less than 10° from the Moon position. Blue: Constant PDF 
in the Hq hypothesis. Red: PDF in the Hi hypothesis. 



observation of the Moon shadow requires large statistics. 
For the work reported here, 884 days of data taken 
between 2007 and 2010 (including periods with 5, 10 and 
12 lines configurations) have been analysed. 

These data have been processed using the standard track 
reconstruction algorithm of ANTARES. As this algorithm 
is optimised for up-going tracks of cosmic origin, an 
algorithm dedicated to the reconstruction of down-going 
atmospheric muons would improve the results presented 
in this paper. 

The event selection criteria used here are the same as 
presented in fTl, i.e. an error estimate /? < 1° and a quality 
of reconstruction of the track A > —5.5. 

B. Monde Carlo 

The point spread function (PSF) used to generate pseudo 
experiments and to compute our search strategy is ob- 
tained from Monte Carlo simulation (MC). Each data 
run is simulated as close as possible to the experimental 
condition (same detector condition, bioluminescence rate, 
etc.), generating both atmospheric muons and atmospheric 
neutrinos. An example of data/MC comparison is repre- 
sented in figure |2] 

The PSF is extracted from this MC using all the events 
reconstructed in a region of 10° around the Moon position 
(figure [3]l. With the chosen selection criteria, the median 
angular resolution is 0.75°. 

Concerning the shape of the event density in the Moon 
region, we assume the event probability density function 
(PDF) is uniform in the abscence of MoorQ we call this 
hypothesis Hq. The PSF in the Moon hypothesis (Hi) is 



obtained by computing the 2D convolution product of the 
PSF with the Moon shape: 

PDF{x\Hq) cx 1 
PDF{x\Hi) oc (1 - PSF ^ 6Moon)ix) 

III. Search strategy 

Different search strategies have been studied, a simple 
on/off search, a test of histograms [7| and the one 
presented here, a likelihood search. The results of these 
different approaches are compatible. 

A. Likelihood method 

For the likelihood search we build a test function based 
on the likelihood ratio using the expected distributions 
corresponding to the Moon hypothesis Hi and no Moon 
hypothesis Hq. The test function used is the logarithm of 
the likelihood ratio: 

^ ^_PDF{n\Hi) 



t 



log 



PDF{n\Ho) 



B. Pseudo-experiments 



In order to test the search method, to estimate its power 
and to compute the significance, pseudo-experiments are 
performed. After generating uniform event distributions in 
an area around the Moon position, a smearing is applied 



1. If we consider a 10° radius region around the Moon 
direction, there is actually a counting rate modulation of ±20%, 
but this is a dipolar modulation which cancels out during the 
likelihood computation, so it is neglected for simplicity. Only 
second or higher order modulations could change the likelihood. 
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Fig. 4: Distributions of the likelihood test function values, for the Hq hypothesis in blue and for Hi in red. 
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Fig. 5: Event density around the Moon location (0, 0). 



ra 300 
■a 

>280 
'm 

g 260 



g 240 
> 

LLI 

220 



200 
180 
160 
140 
120 



0.5 



1.5 



2 2.5 3 
Moon distance (deg) 



Fig. 6: Event density as a function of the distance to the 
Moon. 



using the PSF, followed by a shadowing which is or is not 
applied for the events originating from the Moon direction. 
The pseudo experiments presented here are generated with 
the same event statistics as available in the 2007-2010 data 
sample. 

The distributions of the values of the test function obtained 
with the pseudo experiments are presented in figure |4] For 
a given value of the test function, its p-value (probability 
to obtain a value at least as important from a background 
fluctuation) is computed using the Hq distribution (blue). 
The significance of this p-value can then equivalently be 
expressed in term of sigma. 

The Hi distribution (red) is used to estimate the power of 
the test. With the current event statistics, we could expect 
2.llQ'g(T. With 2, 5 and 10 times this statistics, we obtain 
respectively 2.8 ± 0.9, 4.2 ± 0.9 and 5.9 ± 0.9ct. 



IV. Preliminary results 



The experimental event density in the Moon region is 
represented in figure [5] and the density as a function of 
the distance to the Moon position is shown on figure |6] 

The value of the test function obtained experimentally for 
the Moon position is Q = 1.17, which has a p-value 
of 0.007. This corresponds to 2.7a, which is within the 
expected range. This is however not significant enough to 
unambiguously reject the Hq assumption. 

Scanning over directions around the Moon position (see 
figure [7]i, we find similar hot spots; which are compatible 
with statistical fluctuations. 
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Fig. 7: Significance of the search around the Moon loca- 
tion (0,0). 



V. Conclusion and outlook 

A search for the observation of the Moon shadow with 
the ANTARES neutrino telescope has been performed 
on a dataset corresponding to 884 days of livetime. The 
resulting Moon shadow significance, 2.7(t, is within the 
expected range and compatible with other analysis |7|. At 
this point, this result is not strong enough to put interesting 
constrains on ANTARES pointing capabilities. 

The increase of statistics or the development of a track 
reconstruction algorithm optimized for down-going atmo- 
spheric muons should improve this results. 

ANTARES will also perform campaigns of measurement 
of atmospheric showers by mean of a surface array in- 
stalled on a boat above the detector, allowing to compare 
the direction of muons in coincidence with ANTARES. 
This will provide an additional check of the ANTARES 
pointing performance. 
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Abstract. The AMADEUS system is an integral part 
of the ANTARES neutrino telescope in the Mediter- 
ranean Sea. The project aims at the investigation 
of techniques for acoustic neutrino detection in the 
deep sea. Installed at a depth of more than 2000 m, 
the acoustic sensors of AMADEUS are based on 
piezo-ceramics elements for the broad-band recording 
of signals with frequencies ranging up to 125kHz. 
AMADEUS was completed in May 2008 and comprises 
six "acoustic clusters", each one holding six acoustic 
sensors that are arranged at distances of roughly Im 
from each other. The clusters are installed with inter- 
spacings ranging from 15 m to 340 m. Acoustic data are 
continuously acquired and processed at a computer 
cluster where onUne filter algorithms are applied to 
select a high-purity sample of neutrino-like signals. In 
order to assess the background of neutrino-like signals 
in the deep sea, the characteristics of ambient noise 
and transient signals have been investigated. In this 
article, the AMADEUS system will be described and 
recent results will be presented. 

I. Introduction 

Measuring acoustic pressure pulses in huge underwater 
acoustic arrays is a promising approach for the detection 
of cosmic neutrinos with energies exceeding 100 PeV. The 
pressure signals are produced by the particle showers that 
evolve when neutrinos interact with nuclei in water The 
resulting energy deposition in a cylindrical volume of a 
few centimetres in radius and several metres in length 
leads to a local heating of the medium which is instan- 
taneous with respect to the hydrodynamic time scales. 
This temperature change induces an expansion or contrac- 
tion of the medium depending on its volume expansion 
coefficient. According to the thermo-acoustic model [1|, 
lIJl, the accelerated expansion of the heated volume — a 
micro-explosion — forms a pressure pulse of bipolar shape 
which propagates in the surrounding medium. Coherent 
superposition of the elementary sound waves, produced 
over the volume of the energy deposition, leads to a 
propagation within a flat disk-like volume (often referred 
to as pancake) in the direction perpendicular to the axis 



of the particle shower After propagating several hundreds 
of metres in sea water, the pulse has a characteristic 
frequency spectrum that is expected to peak around 
10 kHz 111, H, im. As die attenuation length in sea 
water in the relevant frequency range is about one to two 
orders of magnitude larger than that for visible light, a 
potential acoustic neutrino detector would require a less 
dense instrumentation of a given volume than an optical 
neutrino telescope. 

The AMADEUS project |6| was conceived to perform a 
feasibility study for a potential future large-scale acoustic 
neutrino detector. For this purpose, a dedicated array 
of acoustic sensors was integrated into the ANTARES 
neutrino telescope |7|. In the following, the AMADEUS 
device will be described and recent results will be pre- 
sented. 

II. The ANTARES Detector 

The ANTARES neutrino telescope was designed to detect 
neutrinos by measuring the Cherenkov light emitted along 
the tracks of relativistic secondary muons generated in 
neutrino interactions. A sketch of the detector, with the 
AMADEUS modules highlighted, is shown in Figure [T] 
The detector is located in the Mediterranean Sea at a 
water depth of about 2500 m, roughly 40 km south of 
the town of Toulon at the French coast at the geographic 
position of 42°48'N, 6°10'E. ANTARES was completed 
in May 2008 and comprises 12 vertical structures, the 
detection lines. Each detection line holds up to 25 storeys 
that are arranged at equal distances of 14.5 m along 
the line, starting at about 100 m above the sea bed and 
interlinked by electro-optical cables. A standard storey 
consists of a titanium support structure, holding three 
Optical Modules (each one consisting of a photomultiplier 
tube inside a water-tight pressure-resistant glass sphere) 
and one cylindrical electronics container 

A 13th line, called Instrumentation Line (IL), is equipped 
with instruments for monitoring the environment. It holds 
six storeys. For two pairs of consecutive storeys in the IL, 
the vertical distance is increased to 80 m. Each line is fixed 
on the sea floor by an anchor equipped with electronics 
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Fig. 1: A sketch of the ANTARES detector. The six 
acoustic storeys are highUghted and a photograph of a 
storey in standard configuration is shown. L12 and IL 
denote the 12th detection line and the Instrumentation 
Line, respectively. 

and held taut by an immersed buoy. An interlink cable 
connects each line to the Junction Box from where the 
main electro-optical cable provides the connection to the 
shore station. 

III. The AMADEUS System 

Within the AMADEUS system |6|, acoustic sensing is 
integrated in the form of acoustic storeys that are mod- 
ified versions of standard ANTARES storeys, in which 
the Optical Modules are replaced by custom-designed 
acoustic sensors. Dedicated electronics is used for the 
amplification, digitisation and pre-processing of the ana- 
logue signals. Figure |2] shows the design of a standard 
acoustic storey with hydrophones. Six acoustic sensors per 
storey were implemented, arranged at distances of roughly 
1 m from each other. The data are digitised with 16 bit 
resolution and 250 k samples per second. 




Fig. 2: Drawing of a standard acoustic storey, or acoustic 
cluster, with hydrophones. 



The AMADEUS system comprises a total of six acoustic 
storeys: three on the IL, which started data taking in 
December 2007, and three on the 12th detection line 
(Line 12), which was connected to shore in May 2008. 
AMADEUS is now fully functional and routinely taking 
data. 

Two types of sensing devices are used in AMADEUS: 
hydrophones Acoustic Modules |6 |. The acoustic sen- 
sors employ in both cases piezo-electric elements for the 
broad-band recording of signals with frequencies ranging 
up to 125 kHz. For the hydrophones, the piezo elements 
are coated in polyurethane, whereas for the Acoustic 
Modules they are glued to the inside of standard glass 
spheres which are normally used for Optical Modules. 

The measurements presented in this article were done with 
the hydrophones. Their calibration will be discussed in 
Sec.lWl 

The AMADEUS on-shore trigge^^ searches the data by 
an adjustable software filter; the events thus selected are 
stored to disk. This way the raw data rate of about 
1.5TB/day is reduced to about lOGB/day for storage. 
Currently, three trigger schemes are in operation |6|: A 
minimum bias trigger which records data continuously 
for about 10 s every 60min, a threshold trigger which is 
activated when the signal exceeds a predefined amplitude, 
and a pulse shape recognition trigger For the latter, a 
cross-correlation of the signal with a predefined bipolar 
signal, as expected for a neutrino-induced shower, is 
performed. The trigger condition is met if the output of the 
cross-correlation operation exceeds a predefined threshold. 
For the latter two triggers, the thresholds are automatically 
adjusted to the prevailing ambient noise and the condition 
must be met in at least four sensors of a storey. 

IV. Ambient Noise 

Ambient noise, which can be described by its charac- 
teristic power spectral density (PSD), is caused by en- 
vironmental processes and determines the minimum pulse 
height that can be measured, if a given signal-to-noise 
ratio (SNR) can be achieved with a search algorithm. To 
measure the ambient background at the ANTARES site, 
data from one sensor on the IL07 taken from the beginning 
of 2008 until the end of 2010 were evaluated. After 
quality cuts, 27905 minimum bias samples (79.9% of 
the total number recorded in that period) were remaining 
for evaluation, each sample containing data continuously 
recorded over a time-span of ^lOs. For each of these 
samples, the noise PSD (units of V^/Hz) was integrated in 
the frequency range / = 10 — 50 kHz, yielding the square 
of the ambient noise for that sample, as quantified by 
the output voltage of the hydrophone. Preliminary studies 
using the shower parametrisation and algorithms from f4l 
indicate that this range optimises the SNR for the expected 
neutrino signals. 

1. While this functionality might be more commonly denoted 
as filtering, it is ANTARES convention to refer to the "on-shore 
trigger". 
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The frequency of occurrence distribution of the resulting 
noise values, relative to the mean noise over all samples, 
is shown in Fig. [3] Also shown is the corresponding 
cumulative distribution. For 95% of the samples, the noise 
level is below 2(crnoiso), demonstrating that the ambient 
noise conditions are stable. 
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Fig. 3: Frequency of occurrence distribution for the am- 
bient noise in the range 10 — 50 kHz, relative to the mean 
ambient noise recorded over the complete period of three 
years that was used for the analysis (left scale, filled 
histogram). Also shown is the cumulative distribution, nor- 
malised to the total number of entries of the distribution 
(right scale, dotted line). 

All sensors have been calibrated in the laboratory prior to 
deployment. The absolute noise level can be estimated 
by assuming a constant sensor sensitivitjj^] of —145 ± 
2dBre lV/|J.Pa. With this value, the mean noise level is 
(fnoisc) — 10.1^2 '^P^ ^i'^h the median of the distribution 
at 8.1 mPa. 

Currently, the detection threshold for bipolar signals corre- 
sponds to a SNR of about 2 for an individual hydrophone. 
For this SNR, the median of the noise distribution corre- 
sponds to a signal amplitude of ^15mPa, equivalent to a 
neutrino energy of 1.5 EeV at a distance of 200 m [3]. By 
applying pattern recognition methods that are more closely 
tuned to the expected neutrino signal, this threshold is 
expected to be further reduced. 

V. Transient Sources 

Transient sources, e.g. from sea mammals, may create 
signals containing the characteristic bipolar pulse shape 
that is expected from neutrino-induced showers. Further- 
more, the pulse shape recognition trigger (see Sec. IIIi 
selects events with a wide range of shapes. Therefore, 
a classification scheme is being developed that selects 
neutrino-like events and suppresses background events 
with high efficiency. After selecting neutrino candidates on 
the level of a storey, measurements from multiple storeys 
can be combined to search for patterns that are compatible 
with the characteristic "pancake" pressure field resulting 
from a neutrino interaction. 



A. Source Position Reconstruction 

The sensors within a cluster allow for efficient triggering 
of transient signals and for direction reconstruction. The 
combination of the direction information from different 
acoustic storeys yields the position of an acoustic source. 
Figure |4] shows the reconstructed directions of all sources 




Fig. 4: Map of directions of sources as reconstructed with 
an acoustic storey on Line 12. Zero degrees in azimuth 
correspond to the north direction, the polar angle of zero 
corresponds to the horizon of an observer on the acoustic 
storey. At the bottom, the signals of the emitters of the 
ANTARES positioning system are visible. 

that were triggered during a period of one month. The 
dark bands of increased acoustic activity can be associated 
with shipping routes and points of high activity with the 
directions of local sea ports. It is obvious from Fig. |4]that 
a fiducial volume for the determination of the background 
rate of bipolar events must exclude the sea surface. 

B. Signal Classification 



The pulse shape recognition trigger described in Sec. Ill 
selects a wide range of events each of which can be 
allocated to one of four classes: Genuine bipolar events 
that are compatible with signals expected from neutrinos 
("neutrino-like events"), multipolar events, reflections of 
signals from the acoustic emitters of the ANTARES 
positioning system and random events, where the latter 
class contains all events that do not fit into any of the other 
classes. For the classification, simulated signals represent- 
ing the four classes in equal proportions were produced 
and a set of features extracted which are highly discrim- 
inant between the classes. This feature vector is then fed 
into a machine learning algorithm |8|. Classification is 
performed for the signals from individual hydrophones. 
Subsequently, the results from individual hydrophones are 
combined to derive a classification for a given acoustic 



2. The ambient noise is originating mainly from the sea 
surface and hence displays a directivity which has to be folded 
with the variations of the sensitivity over the polar angle to obtain 
an effective average sensitivity. For the results presented here, 
the noise has been assumed to be isotropic. 
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storey. Several algorithms were investigated , the best of 
which yielded a failure rate (i.e. wrong decision w.r.t. 
simulation truth) at the 1% -level when applied to the two 
signal classes "neutrino-like" and "not neutrino-like". 

VI. Monte Carlo Simulations 

Monte Carlo simulations based on fS), PI are currently 
being implemented for the AMADEUS detector setup. 
Figure |5] shows the simulated density of the energy depo- 
sition of a 10^° GeV hadronic shower, projected into the 
xz-plane. The z- and x-coordinates denote the directions 
along the shower axis and a direction orthogonal to the 
shower axis, respectively. It is mostly the radial energy 
distribution within the shower which is responsible for the 
shape and amplitude of the acoustic pulse that is observed 
in the far field. The resulting pulse is shown in Fig. |6] The 
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Fig. 5: Density of the energy deposition of a 10^° GeV 
hadronic shower resulting from a neutrino interaction, 
projected from a three-dimensional distribution upon the 
xz-plane. Bin sizes are 0.01 m in x and 0.1 m in z. 

corresponding neutrino interaction was generated such 
that the centre of the hadronic shower for a vertically 
downgoing neutrino lies within the same horizontal plane 
as a storey denoted "Storey 2", at a distance of 200 m. This 
way, the storey lies within the "pancake" of the pressure 
field. On a storey 14.5 m below that storey, denoted 
"Storey 1", no signal is observed. This configuration 
corresponds to two adjacent acoustic storeys on L12 or 
the two lowermost storeys on the IL07, see Fig. [T] This 
simulation illustrates the characteristic three-dimensional 
pattern expected from neutrino-generated pressure waves. 
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Fig. 6: Simulated acoustic signals as recorded with hy- 
drophones in two acoustic storeys with a vertical spacing 
of 14.5 m. See text for details. For Storey 2, signals from 
two different hydrophones are shown. 



is on the classification of transient bipolar events to min- 
imise the irreducible background for neutrino searches. In 
addition, Monte Carlo Simulations are under development. 
AMADEUS is excellently suited to assess the background 
conditions for the measurement of bipolar pulses expected 
to originate from neutrino interactions. 
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VII. Summary and Conclusions 

Recent results from the acoustic neutrino detection test 
system AMADEUS, an integral part of the ANTARES 
detector in the Mediterranean Sea, have been presented. 
Measurements of the ambient noise at the ANTARES site 
show that the noise level is very stable and at the expected 
level, allowing for measurements of neutrino energies 
down to ^ 1 EeV. The current focus of the analysis work 
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